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INFRARED  OPTICAL  FIBERS 

Gavan  Rosman  of  Telecom  Australia 
Research  Laboratory,  Clayton,  spoke  of  the 
Telecom  Monash  efforts  of  the  past  decade 
to  produce  a  better  optical  fiber.  It  is  diffi¬ 
cult  to  surpass  silica-based  glass.  One  wants 
a  glass  with  low  loss  over  a  wide  frequency 
range,  but  light  loss  is  boxed  in  between  two 
mechanism-scattering  of  light  at  short 
wavelengths  and  multiphonon  absorption 
on  the  long  wavelength  side.  To  broaden 
the  window,  researchers  look  for  materials 
with  a  lower  temperature  of  glass  formation 
and  with  heavier  elements  less  stiffly  bonded. 
Compared  to  silica  the  resulting  glasses  suffer 
in  mechanical  properties.  Nonoxide  glasses 
are  candidates.  Telecom  has  worked  exten¬ 
sively  on  heavy  metal-fluoride  fibers. 

UQUIDS-PERSPECTIVES, 
ADVANCES,  AND  FRONTIERS 

I.K.  Snook,  Royal  Melbourne  Insti¬ 
tute  of  Technology  (Ref  2).  By  experimen¬ 
tal  studies,  particularly  scattering  experi¬ 
ments,  by  large-scale  computer  simulation, 
and  by  analytic  methods  such  as  hard  sphere 
perturbation  methods,  the  static  and  dynamic 
structure  of  liquids  are  now  measured  with 
high  precision  and  are  increasingly  under¬ 
stood.  Theory  and  numerical  calculations 
agree  with  experimental  findings  on  equa¬ 
tions  of  state  and  on  transport  properties  of 
liquids  and  mixtures  such  as  liquid  inert 
gases,  ionic  melts,  liquid  metals,  and  liquid 
water.  Fluids  at  interfaces,  complex  fluids, 
suspensions  of  particles  in  fluids,  and  melt¬ 
ing  at  interfaces  are  among  areas  of  current 
research. 


THE  STRUCTURE  OF  UQUID 
SURFACES  AND  WETTING 
PHENOMENA 

David  Beagelhole,  Victoria  University, 
Wellington,  New  Zealand.  On  a  melting 
solid  in  equilibrium  with  its  vapor  there  is  a 
liquid  surface  layer.  How  the  liquid  density 
varies  through  this  transition  region  or  even 
how  thick  this  layer  is  are  not  easy  issues  to 
resolve  experimentally.  For  thicknesses  small 
compared  to  the  wavelength  of  light  (not  the 
case  in  the  interesting  critical  region),  ellip- 
sometry  measures  the  spatial  integral  of 
[(dielectric  constant  of  the  liquid  minus  that 
of  the  solid)  times  (dielectric  constant  of  the 
liquid  minus  that  of  the  vapor)].  Character¬ 
ized  by  surface  layer  thickness,  liquids  appear 
to  fall  systematically  into  certain  definite 
empirical  groups.  Measurements  have  also 
been  made  on  CC14  at  the  critical  tempera¬ 
ture,  where  the  density  of  the  liquid  equals 
the  density  of  the  vapor.  In  wetting,  a  thin 
layer  of  a  liquid  spreads  over  another  liquid 
or  solid.  By  ellipsometry  the  microscopic 
thickness  of  such  layers,  the  spread  of  drops, 
and  capillary  phenomena  can  be  studied. 
Water  wetting  layers  appear  to  be  present 
on  ice  below  0  °C-the  surface  melts  first. 
Inert  impurities  have  a  large  effect  on  sur¬ 
face  melting  and  on  the  (temperature  depen¬ 
dent)  thickness  of  the  premelting  layer. 

EXPERIMENTAL  VERIFICATION 
OF  THE  AHARONOV-CASHER 
EFFECT  FOR  NEUTRONS  WITH  A 
CRYSTAL  INTERFEROMETER 

Goeffrey  I.  Opat,  School  of  Physics, 
University  of  Melbourne.  The  Aharonov- 
Bohm  effect,  a  startling  confirmation  of  the 
fundamental  role  in  quantum  physics  of 
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it  starts  to  transform  to  a  paired  superfluid 
state.  By  0.1  mK  most  of  the  fluid  is  con* 
densed  in  the  paired  ground  state  and  the 
fraction  of  normal  fluid  is  low.  The  normal 
fluid  then  constitutes  excitations,  quasipar¬ 
ticles,  whose  mean  free  path  is  long  because 
of  their  low  density.  These  ballistic  quasi¬ 
particles  then  collide  with  the  containing 
vessel.  Vibrating  wires  can  be  used  as  ther¬ 
mometers,  heaters,  quasiparticle  generators, 
and  quasiparticle  detectors.  Properties  of 
the  elementary  excitations,  including  their 
intrinsically  nonlinear  damping,  are  being 
elucidated  by  vibrating  wire  experiments. 

HYSTERETIC  LOSSES  IN  HIGH 
T  SUPERCONDUCTORS 

c 

Karl-Heinz  Muller,  CSIRO  Division 
of  Applied  Physics,  Lindfield,  New  South 
Wales  (Ref  6).  In  applications  such  as  ac 
power  transmission,  ac  loss  in  the  supercon¬ 
ducting  state  is  an  important  material  engi¬ 
neering  parameter.  So  far,  unfavorably  high 
ac  losses  have  been  found  in  ceramic  bulk 
high  Tc  superconductors.  In  ordinary  type  II 
superconductors  low  frequency  hysteresis 
loss  is  the  result  of  irreversible  dissipation 
accompanying  motion  of  vortices  of  mag¬ 
netic  flux  due  to  the  Lorentz  force.  This 
does  not  appear  to  be  the  origin  of  the  high 
losses  in  the  ceramics.  Muller  and  coworkers 
show  that  at  low  magnetic  fields  the  high 
losses  in  polycrystalline  YBCO  are  mainly 
pinning  losses  originating  from  Josephson 
vortices  moving  between  grains,  while  at  high 
fields  the  losses  are  due  to  moving  Abrikosov 
vortices  inside  the  grains.  The  minimum  in 
the  ac  loss  versus  dc  field  strength  curve  is 
found  to  be  due  to  the  weak  field  depen¬ 
dence  of  the  intergranular  pinning  force 
density,  not  to  surface  barriers  to  flux  entry 
or  exit  as  in  conventional  superconductors. 


THE  UNIVERSAL  ABSENCE  OF 
MACROSCOPIC  LORENTZ  FORCE 
CONTRIBUTIONS  TO  HIGH  T 
SUPERCONDUCTOR 
ENERGY  DISSIPATION 

G.B.  Smith,  Department  of  Applied 
Physics,  University  of  Technology,  Sydney 
(Ref  7).  Smith,  with  groups  at  the  University 
of  Technology  and  at  CSIRO,  Lindfield,  has 
performed  experiments  that  complement 
the  investigations  of  loss  mechanisms  in  high 
Tc  superconductors  described  in  the  pre¬ 
vious  abstract.  In  conventional  supercon¬ 
ductors  the  principal  dissipation  mechanism 
is  flux  jumps  between  pinning  centers.  The 
high  T  ceramics  are  different.  Lorentz  forces, 
at  least  on  a  macroscopic  scale,  appear  to 
play  no  part  in  dissipation.  Loss  curves  in 
high  Tc  superconductors  are  independent  of 
the  direction  of  current  relative  to  ambient 
magnetic  field  (after  correcting  for  demag¬ 
netization).  If  loss  depended  upon  the 
Lorentz  force  JxB,  it  is  implausible  that  it 
would  be  the  same  with  J _LB  as  with  J  ||  B;  in 
the  latter  case  there  should  be  no  loss. 
Experiments  have  been  performed  upon 
highly  oriented  thin  films,  randomly  oriented 
bulk  samples,  with  current  densities  ranging 
over  more  than  two  orders  of  magnitude,  in 
fields  from  0.01  to  7  tesla,  and  at  tempera¬ 
tures  ranging  down  to  below  the  irreversible 
temperature,  at  which  the  flux  creep  mech¬ 
anism  is  operative  in  conventional  super¬ 
conductors.  Loss  results  from  some  activa¬ 
tion  process,  and  the  activation  energies 
depend  upon  the  magnitudes  of  both  the 
current  and  the  magnetic  field,  but  not  upon 
their  relative  direction.  Localized  fluctua¬ 
tions  at  Josephson  junctions  seem  to  be 
implicated. 


SCANNING  TUNNEUNG  of  the  intermetallic  compound.  By  con- 

MICROSCOPY  OF  POLYMERS  tinued  ball  milling  beyond  the  superplastic 

domain  the  mixture  can  be  made  so  com- 

J.W.  White,  Research  School  of  Chem-  plete,  within  each  particle,  that  the  atoms 

istry,  Australian  National  University,  are  within  a  diffusion  length  of  their  ordered 

Canberra.  Scanning  tunneling  microscopy  location,  but  the  mixture  remains  disordered, 

is  performed  only  with  difficulty  on  weakly  A  somewhat  inhomogeneous  alloy  is  formed, 

conducting  substances.  This  talk  constitutes  can  be  compacted,  and  is  stable.  Its  Gibbs 

the  first  report  of  imaging  by  scanning  tun-  free  energy  is  not  far  from  that  of  a  linear 

neling  microscopy  of  a  polymer  made  con-  combination  (at  that  composition)  of  the 

ducting  by  dissolving  an  amphophilic  sub-  free  energies  of  the  two  terminal  compo- 

stance.  Polyacetylene  can  be  made  soluble  nents.  Since  the  composition  is  within  the 

in  organic  solvents  by  copolymerizing  with  a  range  that  would  have  been  a  liquid  on  the 

soluble  species  such  as  polyisoprene.  Scan-  equilibrium  diagram  except  for  the  exis- 

ning  tunneling  microscopy  reveals  that  the  tence  of  the  intermetallic  compound,  and 

copolymer  molecules  thus  formed  aggre-  since  room  temperature  is  below  the  glass 

gate  in  the  same  way  as  soap  bubbles  but  transition,  continued  mechanical  alloying 

form  lamellar  structures.  A  reasonable  will  often  form  an  amorphous  alloy.  But  the 

explanation  is  that  aggregation  happens  intermetallic  has  a  lower  free  energy.  So 

because  the  molecules  are  amphophilic  (can  when  the  solid  metal  pellet  is  heated  to  a 

have  either  positive  or  negative  valence).  temperature  at  which  diffusion  and  struc- 

Conductivity  increases  markedly  at  about  tural  transformation  can  occur,  a  self- 

20  percent  polyacetylene  component,  pre-  sustaining  solid  state  reaction  takes  place, 

sumably  a  percolation  threshold.  heat  is  evolved-sometimes  enough  to  melt 

the  pellet-and  the  ordered  intermetallic 
METASTABLE  PHASES  OF  compound  forms. 

METALS  AND  ALLOYS 

NUCLEAR  ORIENTATION  STUDY 

Earl  Callen,  Office  of  Naval  Research  OF  THE  QUASI-1 D  ANTIFERRO- 

Far  East  (Ref  8).  Equilibrium  phase  dia-  MAGNET  CsMnCI3*2H20 
grams  can  be  used  to  advantage  to  predict 

and  understand  metastable  behavior.  As  G.J.  Bowden,  School  of  Physics,  Uni- 

one  example  (of  many  surveyed  in  this  review  versity  of  New  South  Wales  (NSW), 

talk)  consider  two  substances  whose  phase  Kensington  (Ref  9).  Bowden,  with  groups  at 

diagram  would  have  had  a  eutectic  but  for  NSW,  Kensington,  and  at  the  Department 

the  existence  of  an  intermetallic  compound.  of  Physics,  University  College,  Australian 

Intermetallic  compounds,  despite  low  Gibbs  Defense  Force  Academy,  Campbell,  has  been 

free  energy,  can  often  be  suppressed  because  using  the  dilution  refrigerator  of  the 

they  are  ordered  and  exist  only  in  a  narrow  Australian  Defense  Force  Academy  to  cool 

composition  range.  Suppose  one  then  CsMnCl3*2H20  to  less  than  15  mK.  The 

mechanically  alloys  at  a  temperature  (room  results  support  the  picture  of  a  very  large 

temperature)  too  low  for  atomic  diffusion  (30  percent)  spin  reduction  in  this  quasi-ID 

powders  of  the  two  components  in  the  ratio  antiferromagnet.  This  is  reinforced  by  a 


random  phase  approximation  calculation. 
(But  of  course  the  calculation  is  based  upon 
die  initial  assumption  of  a  Neel  state.  Though 
it  is  too  far  from  self-consistent  to  provide  a 
quantitatively  reliable  estimate  of  the  spin 
reduction,  the  calculation  does  suggest  that 
it  is  large.)  Bowden  also  sees  evidence  of 
solitons.  In  the  ferromagnet  a  soliton  is  a 
moving  domain  wall;  the  domain  wall  lowers 
magnetostatic  energy.  In  the  antiferromag- 
net  the  soliton  is  a  propagating  boundary  of 
an  antiphase  domain.  There  is  no  magneto¬ 
static  driving  force  for  antiphase  domains, 
only  the  luck  of  the  draw  at  nucleation.  But 
solitons  should  be  present  in  thermodynamic 
equilibrium,  and  should  propagate,  as  to  do 
so  requires  no  energy  in  the  perfect  crystal. 
In  the  antiferromagnet  they  are  not  sup¬ 
pressed  or  otherwise  much  affected  by  a 
magnetic  field,  even  a  spin  flop  field.  This  is 
observed. 
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(1989). 

8.  E.  Callen,  F.  Pettit,  and  P.  Shingu,  Scientific 
Information  Bulletin  14(2),  1  (April/June  1989). 

9.  GJ.  Bowden  and  J.P.D.  Martin,/  Phys.  Condens. 
Matter  1,  179  (1989);  J.P.D.  Martin,  Ph.D.  thesis, 
University  of  New  South  Wales  (1988). 


EariCaflen  is  a  member  of  the  staff  of  the  Office  of  Naval  Research  Far  East.  He  is  a  Professor  Emeritus 
of  The  American  University.  He  has  been  active  in  the  physics  of  magnetoelastic  phenomena  and 
amorphous  magnetism. 
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10TH  AUSTRALASIAN  FLUID 
MECHANICS  CONFERENCE 


HideoYoshihara 


Selected  papers  presented  at  the  10th  A  ustralasian  Fluid  Mechanics  Conference  at 
the  University  of  Melbourne  are  reviewed.  Followingthis  discussion,  conversations 
with  Professors  Graeme  Bird  (direct  simulation  Monte  Carlo  method)  and  Clive 
Fletcher  (computational  fluid  dynamics  and  supercomputers),  both  of  the  University 
of  Sydney,  are  described. 


1 0TH  AUSTRALASIAN  FLUID 
MECHANICS  CONFERENCE 

This  conference,  which  was  held  from 
1 1  to  15  December  1989  at  the  University  of 
Melbourne,  emphasized  fundamental  fluid 
dynamics  much  like  the  American  Physical 
Society  fluid  dynamics  meetings.  The  empha¬ 
sis  of  this  meeting  was  computational  and 
experimental  turbulence.  This  series  of 
meetings,  held  every  3  years,  originally  cov¬ 
ered  more  empirical  engineering  topics  such 
as  hydraulics  and  was  primarily  a  national 
meeting  with  New  Zealand.  During  the  last 
few  conferences  emphasis  was  shifted  to 
more  fundamental  aspects  of  fluid  dynamics 
and  the  meeting  was  expanded  to  include 
international  participation.  The  conference 
was  chaired  by  Professor  Anthony  Perry  of 
the  Mechanical  Engineering  Department, 
University  of  Melbourne,  and  honored 
Professor  P.  Joubert  of  the  University  of 
Melbourne,  a  pioneer  of  Australian  fluid 
dynamics,  upon  his  retirement.  Total  atten¬ 
dance  was  229  with  the  breakdown  as  fol¬ 
lows:  Australia,  154;  United  States,  34; 


Germany,  11;  United  Kingdom,  7;  Japan,  5; 
New  Zealand,  5;  and  India,  Singapore, 
Canada,  Switzer1  ard,  and  Israel,  13. 

Presentations  included  10  invited  talks 
of  60  minutes  duration  given  at  the  begin¬ 
ning  of  the  morning  and  afternoon  sessions 
each  day.  The  contributed  talks  were  of 
20  minutes  duration  presented  in  four  con¬ 
current  sessions  with  six  papers  per  day  in 
each  session.  The  highlight  of  the  conference 
was  the  10  invited  talks  by  outstanding  fluid 
dynamicists.  These  included  George 
Batchelor  (Cambridge  U.),  William  Reynolds 
(Stanford  U.),  Hans  Homung  (Cal  Tech), 
John  Lumley  (Cornell  U.),  Anthony  Leonard 
(Cal  Tech),  Roddam  Narasimha  (National 
Aeronautical  Lab  (NAL)-Bangalore),  Brian 
Cantwell  (Stanford  U.),  Morton  Denn  (U. 
of  Calif.-Berkeley),  Ray  Stalker  (U.  of 
Queensland),  and  J.F.  Clarke  (Cranfield  Inst 
of  Tech.). 

Invited  Papers* 

George  Batchelor  (Cambridge  U.), 
“The  Dynamics  of  a  Mixture  of  Two  Com¬ 
ponents”:  Professor  George  Batchelor  is  an 


•  It  was  fortunate  for  this  report  that  all  of  the  written  versions  of  the  invited  talks  were  in  the  Proceedings,  and  this 
infrequent  accomplishment  is  to  be  credited  to  Professor  Tony  Perry  and  Dr.  M.  Chong,  who  edited  the 
Proceedings. 
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internationally  recognized  researcher  of  Newtonian  flows.  Construction  of  the  stream- 
fundamental  turbulence  and  is  an  alumnus  lines  and  vorticity  lines,  given  the  Navier/ 
of  Melbourne  University.  He  reviewed  the  Stokes  solution,  has  been  the  subject  of 
dynamics  of  fluids  containing  macroscopic  other  interesting  contributions  by  Univer- 

particles  by  dividing  the  flows  into  three  sity  of  Melbourne  researchers  Tony  Perry 

classes.  The  first  is  the  well-developed  case  (chairman  of  the  conference)  and  Dr.  Chong 

where  there  is  no-slip  between  the  gas  and  and  was  also  the  subject  of  several  other 

the  particles.  More  difficult  and  less  estab-  papers  at  the  conference, 
lished  are  the  cases  with  slip.  He  first  reviewed  John  Lumley  (Cornell  U.),  “Low 

the  special  case  of  the  dusty  gas  where  there  Dimensional  Models  of  the  Wall  Region  of 

are  a  number  of  results  available.  The  final  a  Turbulent  Boundary  Layer,  and  the  Pol¬ 
and  not-yet  satisfactorily  formulated  case  bility  of  Control”:  Coherent  features  of  the 

was  the  slip-case  with  fluid  dynamic  interac-  wall  region  of  the  boundary  layer  such  as 

tion  between  the  particulates.  Difficulty  bursts  (generators  of  turbulence)  are 

arises  here  since  the  particulate  velocities  modeled  by  low-order  eigenfunctions  gen- 

must  be  treated  statistically.  Use  of  exam-  erated  by  periodic  boundary  conditions  in 

pies  made  clearer  the  physics  of  these  flows  the  homogeneous  directions.  This  low  wave 

and  their  theoretical  difficulties  and  limita-  number  modeling  then  was  used  to  suggest 

tions.  means  to  reduce  turbulence,  for  example, 

William  Reynolds  (Stanford  U.),  using  polymer  additives. 

“Effects  of  Rotation  on  Homogeneous  Anthony  Leonard  (Cal  Tech),  “New 

Turbulence”:  Professor  Reynolds  is  a  long-  Approaches  to  Problems  in  Fluid  Mechan- 

time  expert  on  fundamental  turbulence  and  ical  Mixing”:  Knowing  a  chaotic  or  turbu- 

is  a  key  member  of  the  NASA-Ames/Stanford  lent  solution,  there  is  then  the  task  of  inter- 

Center  for  Turbulence  Research.  The  pres-  preting  the  mass  of  data,  for  example,  by 

ent  paper  was  motivated  by  a  direct  Simula-  tracking  particles  (transport  modeling)  or 

tion  calculation  showing  that  non-isotropy  interfaces  (for  reacting  flows.)  In  this  pre- 

of  turbulence  was  reduced  by  rotation.  He  sentation  a  first  step  is  taken  considering  a 

used  the  rapid  distortion  theory  for  homo-  time-periodic  planar  case  of  a  vortex  pair  in 

geneous  turbulence  to  evolve  the  mecha-  an  external  oscillating  strain-rate  field.  The 

nism  for  this  significant  effect  of  rotation.  results  are  analyzed  by  studying  lobe-like 

These  results  are  the  stepping  stones  to  structures  relative  to  invariant  curves  in  the 

turbulence  models  applicable  in  vortices  to  Poincare  velocity  map.  There  is  interest  in 

enable  prediction  of  bursting.  extensions  of  this  technique  to  the  three- 

Hans  Homung  (Cal  Tech),  “Vorticity  dimensional  case;  investigations  are  planned 

Generation  and  Transport”:  Professor  Hans  for  the  future. 

Homung  is  also  an  alumnus  of  Melbourne  Roddam  Narasimha  (NAL- 

University  and  is  presently  the  director  of  Bangalore),  “Building  Bridges  Between 

the  Graduate  Aeronautical  Laboratory  Dynamic  Chaos  and  Open-Flow  Turbu- 

(GAL),  replacing  Hans  Liepmann.  Hispaper  ler^t.”:  A  dynamic  system  defined  by  a  set  of 

assumes  the  solution  to  be  known  at  the  no-  three  nonlinear,  first  order  ordinary  differ- 

slip  wall,  and  from  this  he  constructs  the  ential  equations  was  evolved,  producing  chaos 

topology  of  the  vorticity  lines  for  planar 
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with  many  key  features  of  open-flow  turbu¬ 
lence  such  as  the  proper  cascading  of  “tur¬ 
bulent  energy.”  Previous  systems  did  not 
possess  these  important  features.  Professor 
Narasimha  expressed  optimism  in  the  rele¬ 
vance  of  dynamic  system  chaos  to  open-flow 
turbulence. 

Brian  Cantwell,  G.  Lewis,  and  J.  Chen 
(Stanford  U.),  “Topology  of  Three- 
Dimensional  Variable  Density  Flows”:  Pro¬ 
cedures  to  interpret  unsteady  compressible 
flow  fields  obtained  either  experimentally 
or  computationally  are  presented  that  fol¬ 
low  earlier  approaches  by  the  Melbourne 
University  researchers.  Specifically,  the 
topology  of  the  flow  is  obtained  by  analyzing 
singular  points.  Examples  of  the  unsteady 
diffusion  flame  and  a  wake  were  given.  The 
recommended  quantity  to  track  was  the 
Galilean-invariant  pressure  gradient.  Efforts 
such  as  this  are  essential  to  interpret  com¬ 
plex  turbulent  flow  fields  for  evolving 
improved  turbulence  models. 

Morton  Denn  (U.  of  Calif.-Berkeley), 
“Mechanics  of  Anisotropic  Fluids”:  The 
fluid  dynamics  of  an  anisotropic  liquid  flow, 
where  the  anisotropy  is  caused  by  suspen¬ 
sions  of  high-aspect  ratio  particles  or  by 
rigid  high-aspect  ratio  molecules,  is  reviewed. 
These  effects  are  important  for  the  perfor¬ 
mance  of  fiber-reinforced  composite  mate¬ 
rials  or  in  advanced  optical  devices. 

J.F.  Clarke  (Cranfield  Inst,  of  Tech.), 
“Developments  in  Detonation  Theory”:  The 
evolution  of  detonation  waves  is  studied 
using  analytical  and  computational  methods. 
It  was  found  that  the  unsteady  induction 
region  plays  an  important  role  in  the  subse¬ 
quent  detonation  process. 

Ray  Stalker  (U.  of  Queensland), 
“Dissociating  Flows  in  Hypervelocity  Aero¬ 
dynamics”:  Simple  pre-design  level  tools 
for  computing  detached  shock  flows  over 
blunt  noses  were  suggested,  assuming  direct 


dissociation  without  recombination.  Stalker’s 
results,  however,  are  not  applicable  in  prac¬ 
tical  cases  where  the  flow  is  not  inviscid  and 
heat  transfer  is  significant. 

Professor  Stalker,  before  and  after  his 
keynote  talk,  passed  out  brochures  on  his 
private  company,  WBM-STALKER,  which 
designs  and  fabricates  hypervelocity  shock 
tunnels  for  testing  aerospaceplane-type 
configurations  and  scramjet  engines.  It 
presently  has  three  contracts  in  the  United 
States  with  the  Rockwell  Rocketdyne 
Company,  General  Applied  Science  Labo¬ 
ratory  (Long  Island),  and  the  California 
Institute  of  Technology.  Professor  Stalker  is 
highly  experienced  in  this  area,  having  oper¬ 
ated  high  temperature  shock  tunnels  at  the 
University  of  Queensland  for  many  years. 
He  was  invited  by  the  AGARD  Fluid  Dynam¬ 
ics  Panel  in  1987  to  present  an  invited  paper 
on  his  hypersonic  test  experiences.  Most 
impressive  was  the  Rocketdyne  shock  tun¬ 
nel,  which  is  piston-driven  with  an  expan¬ 
sion  tube  following  the  shock  tube.  The 
brochure  gave  the  overall  length  as  over 
100  meters  with  a  weight  of  more  than 
700  tons.  This  facility  is  sized  to  enable 
testing,  for  example,  of  a  full-scale  scramjet 
engine.  The  Cal  Tech  shock  tunnel  is  smaller, 
about  50  meters  long,  and  it  will  be  installed 
on  the  roof  of  the  GAL  building.  The  Mach 
number  and  total  temperature  of  these  facil¬ 
ities  were  not  given,  but  it  was  stated  that 
testing  can  be  carried  out  for  near-orbital 
speeds.  It  was  mentioned  by  Professor 
Homung  that  support  for  the  Cal  Tech  facil¬ 
ity  was  furnished  in  part  by  Rockwell. 

Contributed  Papers 

In  general,  the  quality  of  these  papers 
was  high,  comparing  favorably  in  many  cases 
with  the  quality  of  some  of  the  invited  talks. 
The  subjects  covered  and  the  number  of 
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papers  on  each  subject  are  as  follows:  tur¬ 
bulent  flow  including  transition  (30),  com¬ 
pressible  flow  (15),  acoustics  (12),  stratified 
flows  (9),  computational  fluid  dynamics 
(CFD)  (9),  geophysics  (9),  experimental 
techniques  (9),  multi-phase  flows  (6),  envi¬ 
ronmental  fluid  dynamics  (6),  and  wind 
engineering  (6).  Other  sessions  included 
reacting  flows  and  heat  transfer/convection. 
A  discussion  of  the  presentations  attended 
follows. 

D.  Driver  and  J.  Johnson  (Stanford 
U.),  “Three-Dimensional  Boundary  Layer 
Flow  with  Streamwise  Adverse  Pressure 
Gradient”:  Laser  velocimeter  measurements 
of  the  Reynolds  stresses  were  made  for 
boundary  layers  with  large  crossflow.  Reduc¬ 
tion  of  the  streamwise  Reynolds  stress  by 
the  crossflow  was  quantitatively  determined, 
and  this  effect  will  be  particularly  important 
in  the  shock/boundary  layer  interaction  for 
transonic  swept  wings. 

J.  Flores  (NASA-Ames),  “The  Simu¬ 
lation  of  Hypersonic  Flow  Using  a  Zonal 
Approach”;  D.  Chaussee  (NASA-Ames), 
“Hypersonic  3-D  Flow  Past  Winged  Bodies”; 
N.  Chadeijian  and  J.  Flores  (NASA-Ames), 
“Transonic  Viscous  Flow  Computations 
About  a  Complete  Aircraft  Using  the  Navier/ 
Stokes  Equations”:  These  three  papers,  as 
well  as  that  of  Dr.  T.  Edwards  below,  demon¬ 
strate  the  significant  progress  in  the  calcula¬ 
tion  of  complex  configurations  using  Navier/ 
Stokes  methods.  Computing  costs  must  be 
significantly  reduced  before  such  procedures 
can  become  viable  design  tools.  Unfortu¬ 
nately,  the  authors  were  unable  to  present 
their  papers  which  were,  however,  ably 
presented  by  Dr.  Edwards. 

A.  Hilgenstock,  H.  Vollmers,  and 
U.  Dallmann  (DLR-Gottingen),  “Open 
Separation  and  Vortex  Formation-Numerical 
Simulation  and  Analysis  of  a  Delta  Wing 


Flow”:  Further  building  blocks  for  the 
understanding  of  the  puzzle  of  vortex  burst¬ 
ing  were  given. 

T.  Lim,  T.  Nickels,  and  M.  Chong  (U. 
of  Melbourne),  “A  Note  on  the  Cause  of 
Rebound  in  Vortex  Ring/Wall  Interactions”: 
Water  tank  tests  were  skillfully  carried  out, 
determining  the  wall  boundary  layer  to  be 
the  cause  for  the  rebound. 

K.  Gersten  (U.  of  Bochum,  West 
Germany),  “Some  Open  Questions  in  Tur¬ 
bulence  Modeling  from  Viewpoint  of 
Asymptotic  Theory”:  Large  Reynolds 
number  limit  was  used  to  weed  out  inconsis¬ 
tent  algebraic  turbulence  models. 

T.  Edwards  (NASA-Ames),  “Numer¬ 
ical  Investigation  of  Hypersonic  Exhaust 
Plume/Afterbody  Flow  Fields”:  Timely  and 
important  propulsion/airframe  interactions 
were  calculated  with  the  Navier/Stokes  code. 

VISIT  TO  THE  UNIVERSITY 
OF  SYDNEY 

Professor  Graeme  Bird  (Aeronauti¬ 
cal  Engineering  Department) 

Graeme  Bird,  by  consensus,  is  the 
“premier  expert”  on  the  direct  simulation 
Monte  Carlo  (DSMC)  method.  I  last  talked 
with  Professor  Bird  around  1960  in  San  Diego, 
and  he  has  been  working  on  DSMC  since 
then.  He  will  retire  in  February  1990  from 
the  University  of  Sydney  but  will  continue 
DSMC  research  in  his  own  company,  which 
operates  out  of  a  small  room  in  the  base¬ 
ment  of  his  home,  where  we  met  during  my 
visit  to  Sydney.  He  has  been  supported  by 
NASA-Langley  during  the  past  5  years, 
received  an  American  Institute  for  Aero¬ 
nautics  and  Astronautics  (ALAA)  Award  in 
1989  for  his  DSMC  research,  and  was  recently 
elected  Fellow  of  the  AIAA. 
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In  the  outer  fringes  of  the  earth’s 
atmosphere  (80  to  100  km),  the  flow  about 
a  space  vehicle  cannot  be  treated  as  a  con¬ 
tinuum  but  must  be  treated  by  a  method  that 
takes  into  account  the  molecular  makeup  of 
the  air.  The  DSMC  method  is  just  such  a 
method;  it  traces  the  trajectories,  not  of 
individual  molecules,  but  of  representative 
particles  made  up  of  many  molecules.  The 
key  is  to  model  the  representative  particles 
to  reflect  the  statistical  average  behavior  of 
the  molecules  they  contain.  (Needless  to 
say,  there  is  some  controversy  as  to  the 
proper  formulation  of  the  representative 
particles.)  Graeme  Bird’s  DSMC  method 
has  yielded  impressive  results,  for  example, 
matching  closely  with  Navier/Stokes  results 
for  weak  normal  shock  profiles,  which  is  a 
severe  test  case.  This  method  has  received 
wide  support  in  the  United  States,  in  partic¬ 
ular  by  NASA-Langley  and  aerospace 
companies  such  as  Boeing  and  General 
Dynamics,  who  have  purchased  his  produc¬ 
tion  code.  Though  the  method  entails 
repeated  independent  calculations  for  each 
representative  particle,  extensive  scalar 
preparatory  processing  in  each  determina¬ 
tion  renders  large  vector  supercomputers 
ineffective.  Accordingly,  Professor  Bird 
carried  out  his  calculations  on  the  U.S.  space 
shuttle  at  120  km  from  his  basement  office 
with  a  PC  with  8  MB  memory.  An  essential 
element  in  his  production  code  is  the  adap¬ 
tive  cartesian  mesh,  which  is  used  both  for 
sampling  the  macroscopic  properties  as  well 
as  determining  collision  pairs  to  be  used  in 
his  model.  His  current  NASA  grant  calls  for 
adapting  his  code  to  a  large  supercomputer, 
utilizing  the  faster  speed  and  larger  memory 


to  carry  out  expeditiously  more  extensive 
calculations  such  as  extending  the  altitude 
down  to  100  km  for  the  space  shuttle.  In  the 
near  future  calculations  can  be  expected 
down  to  Knudsen  numbers  overlapping 
Navier/Stokes  solutions,  and  this  would  be  a 
landmark  accomplishment.  (See  Refer¬ 
ences  1  and  2  for  discussions  relating  to 
DSMC  methods.*) 

Professor  Bird’s  calculations  are 
impressive,  but  even  more  impressive  will 
be  the  results  that  can  be  expected  from  him 
during  his  “retirement”  years. 

Dr.  Clive  Fletcher  (Mechanical 
Engineering  Department) 

Dr.  Fletcher  holds  the  rank  of  Reader, 
which  falls  between  Professor  and  Asso¬ 
ciate  Professor.  He  is  on  many  international 
CFD  committees  and  is  a  leading  CFD 
researcher  in  Australia.  I  visited  Dr.  Fletcher 
to  obtain  an  overview  of  the  supercomputer 
environment  in  Australia. 

Fluid  dynamics  research  at  the  Uni¬ 
versities  of  Sydney  and  Melbourne  is  mostly 
of  a  fundamental  nature,  treating  simple 
basic  flows,  so  that  onsite  minicomputers 
largely  serve  their  needs.  There  is,  however, 
a  need  for  supercomputers,  for  example,  for 
turbulence  research,  and  a  Cray  XMP 
(probably  a  low  end  model)  is  available  at 
the  Leading  Edge  Technology  Company  in 
Melbourne  that  can  be  accessed  via  tele¬ 
phone.  This  is  a  private  firm  that  operates  in 
close  collaboration  with  Cray  Research- 
Australia.  Leading  Edge  Technology  is  also 
scheduled  shortly  to  receive  a  YMP. 


*  1.  G  .A.  Bird,  “Perception  of  numerical  methods  in  rarefied  gasdynamics,”  in  Rarefied  Gas  Dynamics:  Theoretical 
and  Computational  Techniques,  edited  by  E.P.  Muntz,  D.P.  Weaver,  and  D.H.  Campbell,  Vol.  18  of  Progress  in 
Astronautics  and  Aeronautics  (AIAA,  Washington,  DC,  1989). 

2.  G.S.  Bird,  “Direct  simulation  of  gas  flows  at  the  molecular  level,”  Communications  in  Applied  Numerical 
Methods  4, 165-172  (1988). 
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The  primary  supercomputer  user  in  Keynote  addresses  were  given  by  represen- 

Australia  is  the  nuclear  industry,  as  is  the  tatives  of  supercomputer  firms:  Blaskovich 

case  in  many  countries.  A  Fujitsu  VP- 100  (Cray),  “Supercomputers:  One  Third  of  the 

(267  MFLOP)  is  installed  at  the  Australian  Environmental  Formula”;  dementi  (IBM), 

Nuclear  Science  and  Technology  Labora-  “Modern  Techniques  in  Computational 

tory.  This  organization  has  sponsored  inter-  Chemistry”;  Tahara  (Fujitsu),  “Trend  of 

national  supercomputer  conferences;  the  Supercomputer  Usage  in  Japan”;  and 

second  of  the  series  was  held  at  Wollongong  T annenbaum  (NNSX,  U.S.  representative 

University  (just  south  of  Sydney)  during  12-  of  NEC),  “Technological  Advances  in  the 

14  December  1989,  concurrent  with  the  Next  Generation  of  Supercomputing.”  The 

Australasian  Fluid  Dynamics  Conference.  Fujitsu  speaker,  N.  Tahara,  participated  in 

Dr.  Fletcher  gave  me  a  copy  of  the  program.  the  Office  of  Naval  Research  Navier/Stokes 

benchmark. 


Hideo  Yoshihara  arrived  in  T okyo  in  April  1 988  for  a  2-year  assignment  as  a  liaison  scientist  for  the  Office 
of  Naval  Research.  His  assignment  is  to  follow  the  progress  of  advanced  supercomputers  and  to  review  and 
assess  the  viscous  flow  simulation  research  in  the  Far  East.  Dr.  Yoshihara  formerfy  was  with  the  Boeing 
Company,  where  he  was  Engineering  Manager  for  Applied  Computational  Aerodynamics.  He  was  also  an 
affiliate  professor  in  the  Department  of  Aeronautics  and  Astronautics  of  the  University  of  Washington,  an  AIAA 
Fellow,  and  a  former  member  of  the  Fluid  Dynamics  Panel  of  AGARD/NATO. 
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JAPANESE  ADVANCED 
TECHNOLOGY  FOR  NOVEL 
SHIP  PROPULSION  CONCEPTS 


Justin  H.  McCarthy 


Four  ongoing  Japanese  development  efforts  on  ship  propulsors  and  powering  are 
described.  They  deal  with  superconducting  magnetohydrodynamic  propulsion, 
propeller  boss  cap  fins  and  contrarotating  propellers  to  increase  propulsive  efficiency, 
and  a  newly  established  research  and  development  program  on  super  high-speed 
vessels. 


INTRODUCTION 

I  returned  to  Japan  in  September  1989 
to  revisit  a  number  of  research  institutions 
that  I  had  visited  3  years  earlier  during  a  14- 
month  tour  of  duty  with  the  Office  of  Naval 
Research  Far  East  in  Tokyo.  Despite  the 
lingering  depression  in  Japanese  as  well  as 
worldwide  shipbuilding,  Japan  continues  to 
allocate  significant  resources  to  ship  research 
and  development  (R&D)  with  emphasis  on 
far-term,  moderate-  and  high-risk  concepts. 
This  article  describes  four  novel  technology 
developments  related  to  ship  propulsion 
which  may  be  of  general  interest  to  readers 
of  this  publication.  As  in  previous  visits,  I 
continue  to  be  impressed  by  the  compe¬ 
tence,  versatility,  and  vision  of  the  Japanese 
professionals  conducting  these  and  other 
investigations. 

The  first  investigation  concerns  R&D 
on  magnetohydrodynamic  propulsion  and 
its  demonstration  on  a  30  m  vehicle  designed 
for  this  purpose.  This  work  is  being  led  by 
the  Japan  Foundation  for  Shipbuilding 
Advancement  under  Ryoichi  Sasakawa.  The 


second  and  third  investigations  are  concerned 
with  increased  propulsive  efficiency  achieved 
primarily  by  reduction  of  kinetic  energy  swirl 
losses  downstream  of  conventional  propel¬ 
lers.  One  concept,  developed  and  exten¬ 
sively  employed  by  Mitsui  O.S.K.  Lines,  uses 
evenly  spaced  short-span  fins  on  the  bossing 
immediately  aft  of  a  propeller.  The  second 
concept,  involving  contrarotating  propel¬ 
lers  and  advanced  machinery  developments, 
was  recently  demonstrated  by  Ishikawajima- 
Harima  Heavy  Industries  Co.,  Ltd.  on  a  bulk 
carrier.  Finally,  a  new  and  ambitious  5-year 
Japanese  national  effort  to  create  a  super 
high-speed  vessel,  “Techno-Superliner  ’93,” 
is  described. 

SUPERCONDUCTING 
MAGNETOHYDRODYNAMIC 
(MHD)  SHIP  PROPULSION 

In  June  1985  the  Japan  Foundation 
for  Shipbuilding  Advancement  (JAFSA) 
formed  a  committee  to  lead  a  pioneering 
R&D  effort  on  superconducting  MHD 
propulsion  for  ships.  The  objectives  of  the 
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R&D  effort  were  to  design,  develop,  manu¬ 
facture,  and  conduct  at-sea  trials  on  a  super¬ 
conducting  coil  for  a  shipboard  MHD  pro¬ 
pulsion  device.  References  1-3  summarize 
progress  through  mid-1989.  Successful 
execution  of  the  effort  has  entailed  the 
construction  of  new  facilities  for  evaluation 
of  MHD  propulsion  at  model  scale.  Three 
facilities  made  of  nonmagnetic,  noncorrod¬ 
ing  materials,  largely  stainless  steel,  were 
constructed  for  saltwater  testing  at  J  AFS  A’s 
Tsukuba  Institute:  (1)  a  tank  23.5  m  long, 
3  m  wide,  and  2  m  in  depth  for  model  MHD 
self-propulsion  experiments;  (2)  a  circulat¬ 
ing  water  channel  with  a  test  section  5  m 
long,  1.6  m  wide,  and  1.85  m  deep  having  a 
maximum  flow  speed  of  1.5  m/s;  and  (3)  a 
1  m  long,  0.1  by  0.2  m  cross-section  electrode 
testing  tube  having  flow  and  speeds  up  to 

3  m/s.  These  facilities  have  proved  indis¬ 
pensable  to  the  R&D  efforts  on  MHD  pro¬ 
pulsion. 

Work  has  proceeded  on  schedule  with 
construction  of  a  MHD  propulsion  vehicle 
to  be  completed  by  Mitsubishi  Heavy 
Industry’s  Kobe  Shipyard  in  April  1990.  The 
at-sea  evaluations  of  the  device  are  planned 
for  November  1990.  The  vehicle  will  have 
an  aluminum  hull  30  m  in  length  with  a  185- 
ton  displacement  and  is  predicted  to  have  a 
top  speed  of  8  knots.  Figure  1  shows  an 
artist’s  picture  of  the  vehicle  and  its  internal 
arrangements.  Propulsion  is  via  twin  MHD 
internal  ducts  having  an  estimated  propul¬ 
sive  efficiency  of  less  than  4  percent.  The 
actual  value  will  not  be  known  until  the  trials 
are  completed  and  many  well  be  less  than 

4  percent. 

MHD  propulsion  is  based  on  Fleming’s 
left-hand  rule  of  electromagnetism  applied 
to  the  magnetic  field  created  in  seawater  by 
a  magnet  attached  to  a  hull.  When  current 
flows  through  a  coil  perpendicular  to  the 
magnetic  field  an  electromagnetic  force 


(EMF)  is  created  in  the  (weakly  conduct¬ 
ing)  seawater  by  the  interaction  between 
the  electric  and  magnetic  fields.  Propulsive 
thrust  results  from  the  seawater’s  reaction 
to  the  electromagnetic  interaction.  Two 
types  of  MHD  propulsion  were  considered, 
external  and  internal.  In  the  external  ver¬ 
sion,  the  EMF  acts  on  the  seawater  exterior 
to  the  vehicle.  In  the  internal  version,  the 
EMF  acts  on  seawater  passing  through  a 
longitudinal  duct  running  through  the  vehi¬ 
cle.  Model  experiments  were  conducted  on 
both  types  of  propulsors  in  1987  at  JAFSA’s 
Tsukuba  laboratory.  Overall  propulsive 
efficiencies  of  no  more  than  1  percent  were 
achieved  in  these  tests,  which  employed  liquid 
helium  to  achieve  superconductivity  in  the 
coils.  A  decision  was  made  to  adopt  the 
internal  type  of  MHD  propulsor  because  it 
allows  minimization  of  the  amount  of  radia¬ 
tion  leaked  to  the  environment. 

The  twin  internal  MHD  propulsor  units 
adopted  in  the  final  design  each  consist  of  six 
dipole  coils  around  six  ducts  arranged  in  a 
ring  so  as  to  reduce  magnetic  radiation  and 
eliminate  magnetic  shielding  except  in  way 
of  the  vehicle’s  control  room.  Electric  power 
comes  from  two  generators  driven  by  two 
diesel  engines.  Reference  4  describes  the 
experimental  evaluation  of  the  hydrodynamic 
performance  of  the  propulsor  units.  A  draw¬ 
ing  of  the  propulsor  cross  section  is  shown  in 
Figure  2  and  the  principal  design  perfor¬ 
mance  measures  are  approximately  as  fol¬ 
lows  (Ref  4): 

Magnetic  Flux  Density  at  Duct  Center  *  tesla  (T) 
Normal  Electrode  Current  2,000  A 

Total  Thrust  (tno  units)  8,000  N 

Total  Capacity  of  Main  Generators  4,060  kW 
Total  Height  of  Propulsion  Equipment  130  tons 

In  estimating  the  total  thrust  a  50  percent 
efficiency  was  assumed  for  converting  Lorentz 
force  into  thrust. 
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Figure  1.  Experimental  MHD  vessel  ( courtesy  of  Japan  Foundation  for  Shipbuilding  Advancement). 
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Figure  2.  Six-duct  MHD  propulsor  arrangement. 


The  low  values  of  efficiency  achieved 
to  date  result  mainly  from  the  Joule  heating 
loss  due  to  the  low  conductivity  of  seawater. 
Efficiency,  which  is  directly  proportional  to 
ship  speed  (U)  and  the  square  of  the  mag¬ 
netic  flux  density  (B),  can  be  significantly 
increased  if  much  more  powerful  magnetic 
fields  can  be  generated.  A  twofold  increase 
in  B  from  4  to  8  T  and  in  speed  from  8  to 
16  knots  would  result  in  a  maximum  overall 
propulsive  efficiency  of  about  30  percent  in 
comparison  with  the  estimated  maximum 
overall  propulsive  efficiency  of  4  percent  for 
the  unit  to  be  tested  at  sea.  Thus,  the  devel¬ 
opment  of  more  powerful,  lightweight 
magnets  is  one  of  the  keys  to  efficient  and 
commercially  economical  MHD  ship  pro¬ 
pulsion  in  the  future.  The  availability  of  high 
temperature  superconducting  materials 
would  also  provide  a  major  advance. 

Other  important  factors  affecting 
future  viability  include  the  choice  of  nonpol¬ 
luting  electrode  materials  for  long  life,  a 
subject  proprietary  to  JAFSA,  and  the  con¬ 
trol  of  bubbles  and  large-scale  turbulence 
associated  with  the  MHD  process.  The 


presence  of  copious  amounts  of  bubbles 
and  associated  noise  would  be  unacceptable 
in  naval  applications.  Experimental  research 
investigations  are  currently  underway  at 
JAFSA’s  Tsukuba  laboratory  to  evaluate 
electrode  materials,  alternative  electrode 
geometries,  and  bubble  generation.  The 
initial  choice  of  electrode  material  is  titanium 
coated  with  an  oxide,  but  other  materials 
with  reduced  chlorine  generation  are  under 
consideration.  Many  difficulties  have  obvi¬ 
ously  presented  themselves,  but  Japan’s 
JAFSA  leadership  is  continuing  its  signifi¬ 
cant  commitment  to  the  high-risk  develop¬ 
ment  of  MHD  propulsion. 

PROPELLER  BOSS  CAP  FIN 
(PBCF) 

The  PBCF  device  to  increase  propul¬ 
sive  efficiency  has  been  developed  during 
3  years  of  collaboration  between  Kazuyuki 
Ouchi  of  the  Mitsui  OSK  Lines,  Ltd.  and 
Masahiro  Tamashima  of  the  West  Japan 
Fluid  Engineering  Company  in  Sasebo,  Japan 
(Ref  5  and  6).  The  PBCF  concept,  depicted 
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in  Figure  3,  consists  of  small  fins  placed  on 
the  hub  fairwater  (boss  cap)  aft  of  the  pro¬ 
peller.  These  fins  are  designed  to  reduce  the 
swirl  downstream  of  the  propeller,  thereby 
reducing  wake  kinetic  energy  losses  and  at 
the  same  time  reducing  the  tendency  for 
hub  vortex  cavitation.  During  the  past  2  years 
59  vessels,  ranging  in  size  from  a  300-ton 
fishing  boat  to  a  98,000-ton  bulk  carrier, 
have  been  fitted  with  PBCF  devices.  Depend¬ 
ing  on  the  ship,  efficiency  gains  of  between 
2  and  8  percent  have  been  realized  in  sea 
trials,  with  an  average  gain  of  about  4  per¬ 
cent.  Model  experiments  have  shown  signif¬ 
icant  delays  in  the  onset  of  hub  vortex  cavi¬ 
tation. 

One  of  the  surprises  in  the  Japanese 
work  has  been  the  discrepancy  between 
efficiency  gains  realized  in  model  self¬ 
propulsion  experiments  and  at  full  scale,  the 
latter  being  up  to  6  percentage  points  higher 
than  the  former.  No  satisfactory  explana¬ 
tion  for  the  difference  has  been  found,  but  it 
is  believed  to  be  caused  by  Reynolds  num¬ 
ber  “scale  effects”  between  model  and  full- 
scale  test  conditions.  Because  of  a  factor  of 
100  difference  in  model-  and  full-scale 


Figure  3.  Propeller  boss  cap 


Reynolds  numbers,  the  hull  boundary  layer 
flow  into  the  propeller  will  differ  in  the  hub 
region  and,  in  the  case  of  the  model,  laminar 
flow  and  separation  may  occur  on  the  back 
side  of  the  model  propeller  blade  near  the 
hub.  Both  effects  could  degrade  the  perfor¬ 
mance  of  the  PBCF  device  at  model  scale, 
but  quantitative  resolution  of  the  scaling 
problem  can  be  achieved  only  by  model  and 
ship  velocity  measurements  in  the  region  of 
the  PBCF.  Full-scale  velocity  measurements 
are  difficult  and  have  been  conducted  rarely. 
As  a  result  Ouchi  and  Tamashima  will  rely 
simply  on  the  empirical  finding  based  on 
extensive  data,  that  propulsive  efficiency 
gains  with  a  PBCF  device  at  full  scale  will  be 
about  double  those  measured  at  model  scale. 

In  order  to  better  understand  the 
powering  performance  of  PBCF  devices, 
Gearhart  and  McBride  of  the  Applied 
Research  Laboratory  at  Pennsylvania  State 
University  have  conducted  detailed  model 
flow  and  powering  measurements  that  they 
have  theoretically  analyzed  (Ref  7).  For  the 
particular  case  examined  the  following  ledger 
of  PBCF  effects  was  constructed: 


fin  (PBCP)  arrangement. 
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Efficiency 

Change 


Factor  (X) 


Increase  in  thrust  due  to  pressure 

rise  on  hub  +2 

Decrease  in  torque  due  to  swirl 

reduction  +4 

Added  frictional  drag  of  FBCF 

vanes  - 1 

Net  Efficiency  Gain  5 


As  pointed  out  by  Gearhart  and  McBride, 
the  magnitude  of  the  efficiency  gain  is  depen¬ 
dent  on  the  amount  of  swirl  that  is  imparted 
by  the  propeller.  Propeller  designs  that 
consider  turning  in  the  blade  root  sections 
induced  by  secondaiy  flows  will  not  have  a 
strong  root  vortex  and  therefore  will  not 
benefit  from  use  of  a  PBCF  device. 

CONTRAROTATING  PROPELLERS 

In  January  1989  sea  trials  were  suc¬ 
cessfully  completed  for  a  contrarotating 
propeller  (CRP)  system  fitted  to  the  37,000- 
ton  bulk  carrier,  JUNO  (Ref  8).  Develop¬ 
ment  of  this  system,  which  entailed  machin¬ 
ery  innovations  on  epicyclic  gears,  contraro¬ 
tating  bearings,  and  seals,  was  initiated  by 
Ishikawajima-Harima  Heavy  Industries  Co., 
Ltd.  (IHI)  in  1984  in  order  to  reduce  ship 
fuel  costs.  The  recent  trials  demonstrated  a 

15  percent  power  saving  at  a  ship  speed  of 

16  knots  in  either  ballast  or  fully  loaded  con¬ 
ditions.  As  in  the  case  of  the  PBCF  device, 
the  at-sea  power  saving  was  higher  than  had 
been  predicted  from  model  experiments. 
Additional  model  tests  are  planned  to  resolve 
the  discrepancy. 

A  CRP  system  consists  of  two  propel¬ 
lers  on  coaxial  shafts  that  rotate  in  opposite 
directions  (see  Figure  4).  The  theoretical 
lifting-line  methods  used  by  IHI  to  design 
the  propellers  are  described  in  References  9 
and  10.  The  CRP  concept,  pioneered  by 


John  Ericsson  in  the  mid-1800s,  has  long 
been  known  to  provide  significant  efficiency 
gains  for  high-powered,  single-screw  ships. 
An  American  study  (Ref  11)  published  in 
1964  showed  an  efficiency  gain  of  8  percent 
fora  CRP  application  to  a  tanker.  The  CRP 
concept  can  also  result  in  reduced  cavitation 
and  propeller-induced  hull  vibration,  and 
because  of  its  smaller  diameters  it  can  be 
applicable  to  shallow  draft  vessels.  How¬ 
ever,  because  of  the  increased  complexity, 
reliability  problems,  and  cost  of  contraro¬ 
tating  shafting  and  machinery,  there  have 
been  few  applications  of  CRP  systems  to 
ships.  Recent  new  gear,  shafting,  and  seal 
technology  advanced  at  IHI  has  removed 
many  of  the  reliability  concerns  and  it  is 
expected  that  the  higher  initial  cost  of  the 
CRP  system  may  one  day  be  outweighed  by 
the  fuel  savings  realized  over  the  lifetime  of 
a  ship. 

The  features  of  the  contrarotating  gear 
system  developed  by  IHI  are  shown  in 
Figure  5.  The  system,  driven  by  a  single 
engine  and  designed  to  produce  equal  rpms 
and  torques  on  the  forward  and  aft 
contrarotating  propellers,  consists  of  star- 
compound  epicyclic  gears,  oil  lubricated  roller 
bearings,  seals,  and  a  Geislinger  elastic 
coupling  for  torsional  vibration  isolation. 
The  epicyclic  gear  system  consists  of  a  sun 
gear  on  the  inner  propeller  shaft  that  drives 
two  planet  gears  in  a  direction  opposite  to 
the  sun  gear.  The  planet  gears  drive  a  ring 
gear  on  the  outer  shaft  in  a  direction  opposite 
to  the  inner  shaft  to  produce  contrarotation. 
The  shaft-mounted  gear  arrangement  was 
adopted  for  protection  of  the  gears  during 
hull  deformation.  IHI  believes  that  one  of 
the  most  important  technical  achievements 
is  the  design  of  the  contrarotating  radial 
roller  bearing  between  the  contrarotating 
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tested  for  500  hours  to  confirm  its  expected 
satisfactory  performance.  No  problems  have 
been  encountered  during  at-sea  operations. 


shafts.  Tapered  and  static-pressure-type 
bearings  were  considered  but  rejected  for  a 
number  of  reasons.  Prior  to  installation  on 
JUNO ,  the  entire  CRP  system  was  shop 


Figure  4.  Stem  arrangement  of  contrarotating  propellers  (CRP). 
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SUPER  HIGH-SPEED  VESSEL  - 
‘Techno-Superliner  ’93” 

An  ambitious  new  5-year  program, 
initiated  in  April  1989 by  Japan’s  Ministry  of 
Transport  and  involving  an  association  of 
Japan’s  principal  shipbuilders  and  govern¬ 
ment  laboratories,  is  aimed  at  developing  a 
super  high-speed  vessel  for  transporting  high 
value  cargo  (Ref  12).  The  vessel  is  called 
“Techno-Superliner  ’93”  (TSL)  and  is  to  be 
demonstrated  in  1993  on  a  large-scale  model, 
from  20  to  30  m  in  length.  The  prototype 
TSL  is  projected  to  be  put  into  service  before 
the  end  of  the  century.  The  goals  call  for  a 
50-knot  speed,  1,000-ton  payload,  and  500- 
nautical  mile  range.  TSL  is  to  be  operable  in 
the  ocean  up  to  sea  state  6  (significant  crest- 
to-trough  wave  heights  up  to  6  m).  It  is 
viewed  as  providing  an  alternative  to  air¬ 
plane  and  merchant  ship  transport  for  carry¬ 
ing  more  lightweight  cargo  than  planes  in 
shorter  times  than  ships.  To  set  the  stage  for 


TSL  a  high-speed  craft  symposium,  orga¬ 
nized  by  the  Society  of  Naval  Architects  of 
Japan,  was  held  in  June  1989  (Ref  13). 

A  key  consideration  in  TSL’s  design  is 
the  minimization  of  wave  making  resistance 
by  having  its  main  hull  rise  above  the  sea 
surface  at  high  operational  speeds  of  up  to 
50  knots.  Two  types  of  craft  are  being 
explored.  (1)  In  the  familiar  surface  effect 
ship  (SES)  the  hull  rides  on  an  air  cushion  on 
the  water  surface.  The  air  is  trapped  between 
shallow-draft  side  walls  and  flexible  skirts  at 
the  bow  and  stem.  (2)  A  small  waterplane 
area  twin  hull  (SWATH)  configuration  has 
submerged  foils  that  lift  the  main  hull  above 
the  ocean  surface.  In  this  concept  relatively 
thin  vertical  support  struts,  two  on  each  side, 
pass  through  the  ocean  surface  connecting 
the  main  hull  to  twin  underwater  hulls  that 
contain  the  propulsion  machinery  and  sup¬ 
port  the  lifting  foils.  Figure  6  shows  an 
artist’s  picture  of  this  concept. 


Figure  6.  Techno-Superliner  ’93  (TSL). 
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The  SES  concept  is  being  evaluated  at 
Mitsubishi  and  Mitsui  Heavy  Industries;  the 
SWATH  concept  is  being  evaluated  by  an 
association  of  Kawasaki,  IHI,  NKK, 
Sumitomo,  and  Hitachi  Heavy  Industries. 
For  both  hull  concepts  water-jet  propulsion 
with  gas  turbine  and/or  superconducting 
motor  drives  is  being  considered.  The  first 
3  years  of  the  national  effort  will  deal  with 
preliminary  design,  hydrodynamic  perfor¬ 
mance  prediction,  structures,  (new)  mate¬ 
rials,  and  the  propulsion  system.  The  last 
2  years  will  be  devoted  to  the  ride  control 
system  and  a  large-scale  model  demonstra¬ 
tion. 

If  successful,  the  high-risk  TSL  effort 
will  provide  transit  times  of  under  12  hours 
within  Japan,  between  Tokyo  and  Kyushu 
or  Hokkaido,  and  from  Nagasaki  to  Shanghai, 
China.  The  TSL  effort  is  sure  to  be  followed 
closely  by  all  maritime  nations. 

POSTSCRIPT 

Three  years  ago  I  reported  on  two 
other  Japanese  propulsor  concepts  (Ref  14 
and  15)  that  were  undergoing  evaluation  at 
model  scale  at  that  time:  the  off-center 
placement  of  propellers  on  single-screw  ships 
and  the  “wave-devouring  propulsor”  (WDP). 
Both  concepts  have  been  recently  evaluated 
at  sea  with  success,  the  former  by  NKK 
Corporation  on  a  233,000-ton  ore  carrier 
and  the  latter  by  Hitachi  Zosen  Corporation 
on  a  16  m  fishing  vessel. 

Trials  for  the  off-center  propeller, 
placed  1.3  m  to  starboard  of  the  ship’s  cen¬ 
terline  and  rotating  clockwise  so  as  to  oppose 
boundary  layer  vorticity  rotating  in  a  counter¬ 
clockwise  direction,  increased  propulsive 
efficiency  to  the  extent  that  annual  fuel  costs 
will  be  reduced  by  about  $75,000  (Ref  16). 
The  second,  more  unorthodox,  concept 


consists  of  a  bow-mounted,  spring-loaded 
horizontal  foil  that  is  tuned  to  oscillate  when 
ocean  waves  cause  the  ship  to  pitch,  gener¬ 
ating  a  net  thrust  force  on  the  foil  and  reduc¬ 
ing  the  pitch  motions  of  the  ship.  It  is 
intended  as  an  energy-saving,  auxiliary  pro¬ 
pulsion  device.  In  the  recent  sea  trials  in 
head  seas  the  WDP  produced  a  speed  of 
about  2  knots  on  its  own  and  reduced  pitch 
motions  by  about  40  percent  (Ref  17). 
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increasing  as  the  level  of  COz  drops 
(S.  Miyachi,  Univ.  of  Tokyo).  The  dual  activ¬ 
ities  of  the  enzyme  rubisco  also  vary  with  the 
concentration  of  COz.  At  high  levels  of  C02 
reduction  is  favored  over  photo-oxidation 
reaction.  The  dual  action  and  responses  of 
rubisco  in  marine  algae  are  yet  to  be  care¬ 
fully  studied.  It  is  thus  too  early  to  evaluate 
the  possible  effects  on  photosynthetic  effi¬ 
ciency  of  these  algae  by  the  increase  in  C02 
content  of  the  atmosphere  (N.E.  Tolbert, 
Michigan  State  Univ.). 

The  activities  within  the  food  web 
determine  the  fate  of  the  reduced  carbon; 
however,  also  significant  is  the  deposition  of 
intra  and  extracellular  carbon  as  calcium 
carbonate.  These  calcite  and  aragonite 
particles  rain  down  from  the  surface  waters 
to  be  deposited  onto  the  seafloor 
(M.  Okazaki,  Gakugei  Univ.,  Tokyo; 
P.  Westbroek,  Leiden  Univ.,  The 
Netherlands).  Whether  the  rate  of  diffusion 
of  C02  into  the  ocean  would  compensate  for 
the  rapid  rise  in  atmospheric  C02  is  still 
being  argued.  In  this  context  I  was  surprised 
that  the  role  of  the  specialized  living  com¬ 
munity  that  occupies  the  upper  surface  of 
the  ocean,  the  Neuston,  was  not  mentioned. 
There  was  a  flurry  of  interest  in  the  Neuston 
about  20  years  ago  and  I  think  the  time  is 
ripe  for  its  reawakening.  It  is  likely  that  the 
Neuston  organisms  contain  a  high  concen¬ 
tration  of  carbonic  anhydrase,  which  will 
accelerate  the  solubilization  of  COz  across 
this  interphase. 

I  did  not  find  among  the  papers  that 
were  presented  at  the  sessions  that  I  attended 
many  exciting  or  novel  ideas,  but  I  think  it  is 
worthwhile  to  summarize  the  type  of  work 
that  is  presently  being  undertaken  under  the 
title  of  marine  biotechnology. 


DISCUSSION  OF 
CONFERENCE  PAPERS 

Microorganisms 

Examples  of  the  potentials  for  bio¬ 
technology  of  the  microbes  from  extreme 
environments  were  presented.  W.  Reichardt 
(FRG)  reported  on  the  use  of  bacteria  iso¬ 
lated  from  cold  ocean  waters  as  sources  of 
enzymes  for  efficient  decomposition  of 
organic  substances  at  low  temperatures.  At 
the  other  extreme  the  use  of  microbes  from 
the  deep  sea  hydrothermal  vents  as  sources 
of  enzymes  was  discussed  by  D.  Prieur 
(CNRS,  France)  and  Y.  Sako  et  al.  (Kyoto 
Univ.).  Sako  reported  on  the  isolation  of  the 
enzyme  PEP-carboxylase  from  a  bacteria 
that  was  derived  from  a  submarine  thermal 
vent.  This  enzyme  was  optimally  active  at 
80  °C.  R.  Colwell  reported  on  the  isolation 
of  a  large  number  of  compounds  from  marine 
sponges  and  tunicates;  these  are  being 
screened  for  their  antitumor  activity.  Other 
reports  dealt  with  finding  of  antibacterial, 
antiviral,  or  toxins  produced  by  isolated 
marine  microbes. 

An  interesting  paper  by  K.  Kogure  et 
al.  (Univ.  Tokyo)  addressed  the  question  of 
the  origin  of  the  toxin  tetrodotoxin  (TTX), 
which  is  found  in  various  marine  animals 
and  is  particularly  infamous  in  making  the 
“puffer”  fish  highly  toxic.  Kogure  et  al., 
using  a  sensitive  in  vitro  bioassay,  showed 
that  a  large  number  of  bacterial  strains  iso¬ 
lated  from  seawater,  sea  sediments,  and 
marine  animals  are  synthesizing  TTX  or 
related  sodium-channel  blocking  agents. 
These  toxins  accumulate  in  various  marine 
animals  via  the  complex  food  web  of  the  sea. 
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Other  studies  pointed  to  the  value  of 
the  bacterial  flora  for  the  development  and 
growth  of  fish  in  the  larval  stages  (M.  Maeda, 
Univ.  of  Tokyo)  or  adults  (H.  Sugita  et  al., 
Nihon  Univ.,  Tokyo).  Thus,  growth  was 
improved  when  bacterial  strains  that  synthe¬ 
size  vitamin  B12  were  present. 

Microalgae 

Mass  cultivation  of  microalgae  for  the 
production  of  specific  high  value  substances 
is  already  a  commercial  enterprise  in 
Australia  and  Israel.  The  halophilic  alga 
Dunaliella  is  being  cultivated  primarily  for 
its  high  content  of  beta-carotenes.  The  natural 
9-cis  isomer  found  in  this  organism  is  reported 
to  be  much  more  biologically  active  than  the 
all-trans  synthetic  isomer.  Conditions  for 
maximizing  pigment  contents  of  the  cells 
are  being  studied  (A.  Ben  Amotz,  Israel; 
M.A.  &  L.J.  Borowitzka,  Australia). 

During  the  oil  crisis  there  was  a  flurry 
of  interest  in  several  species  of  microalgae 
that  produce  and  store  hydrocarbons  in  their 
cells.  Two  Japanese  laboratories  are  appar¬ 
ently  still  hard  at  work  on  cultivation  of  the 
green  alga  Botryococcus  for  the  production 
of  high  quality  hydrocarbons  (S.  Okada  and 
K.  Yamaguchi,  Univ.  of  Tokyo;  H.  Iwamoto 
et  al.,  Meiji  Univ.,  Tokyo).  Extracellular 
hydrocarbons  were  produced  when  the  cul¬ 
tures  grew  on  solid  surfaces  as  colonies,  but 
no  such  production  occurred  when  they  were 
grown  suspended  in  medium.  Although 
“high  quality”  gasoline  could  be  produced 
by  catalytic-cracking  of  these  hydrocarbons, 
the  cost  of  the  product  is  still  several  orders 
of  magnitude  too  high. 

A  novel  approach  to  cultivation  of 
carotenoids  producing  microalgae  was  pre¬ 
sented  by  Y.  Cohen  (Israel).  He  cultivates 
the  benthic  biofilms  of  hypersaline  ponds. 


These  films  are  made  up  primarily  of  the 
cyanobacteria  Microcoleus  and  a  commu¬ 
nity  of  other  microbes.  The  harvested  bot¬ 
tom  mat  from  these  ponds  was  made  up  of 
over  25  percent  polysaccharides  and 
2.5  percent  carotenoids.  Cohen  works  in  a 
desert  environment,  where  he  constructed 
small  ponds  of  evaporating  seawater.  The 
hypersaline  brine  that  accumulates  is  over- 
layed  with  light  seawater.  Nutrients  are 
supplied  to  the  lower  brine  layer  by  a  drip- 
irrigation  system. 

Cultivation  of  a  red,  single  cell  alga, 
Porpheridium ,  for  the  production  of  sulfated 
polysaccharides,  pigments  and  polyunsatu¬ 
rated  lipids  was  reported  by  S.  Arad  (Israel). 
She  introduced  a  novel  culture  system  of 
vertical  polyethylene  sleeves.  Contamina¬ 
tion  by  protozoa  was  still  a  major  problem  in 
her  cultures. 

M.  Murakami  and  K.  Yamaguchi 
(Univ.  of  Tokyo)  found  a  potent  antifungal 
agent  produced  by  the  dinoflagellate 
Alexandrium  hiranoi.  The  substance  is  a 
polyether  macrolide;  it  stops  fungal  growth 
at  a  concentration  of  0.5  //g/ml  and  blocks 
the  development  of  sea  urchin  embryos  (the 
bioassay  used  for  searching  for  bioactive 
substances)  at  0.05  /tg/ml. 

T.  Yasumoto  et  al.  (Tohoku  Univ.) 
isolated  a  number  of  polyethers  from  differ¬ 
ent  dinoflagellates  and  measured  their  anti¬ 
fungal  activity.  Seven  out  of  eight  species  of 
dinoflagellates  produced  significant  antifun¬ 
gal  substances.  Y.  Miura  (Osaka  Univ.)  and 
T.  Matsunaga  (Tokyo  Univ.  of  Agric.  & 
Tech.)  are  also  engaged  in  screening  many 
strains  of  microalgae  for  production  of  anti¬ 
biotics.  Two  of  their  strains  produced  yeast 
inhibiting  substances. 

Genetic  manipulation  techniques  for 
microalgae  were  the  subject  of  two  papers. 
The  marine  unicellular  blue-green  alga 


ONRFE  SCI  INFO  BUL 15  (2)  90 


25 


Synechococcus  was  successfully  transformed 
by  the  use  of  the  electroporation  technique 
(H.  Takeyama  et  al.,  Tokyo  Univ.  of  Agric. 
&  Tech.).  Another  paper  by  T.  Kudo  and 
N.  Saga  (Sapporo  Univ.)  suggests  the  use  of 
antibiotic  paper  discs  on  agar  plates  in  search¬ 
ing  for  transformed  or  genetically  modified 
cells. 

Solar  energy-bioconversion  by  algae 
is  also  a  viable  topic  of  research. 

Mass  cultures  of  single  cell  algae  as  a 
renewable  source  of  organic  fuel  were  studied 
for  a  long  time.  A  more  recent  approach 
was  the  possibility  of  producing  fuels  such  as 
oxygen  and  hydrogen  gas  by  photolysis  of 
water.  Conditions  for  the  photosynthetic 
production  of  hydrogen  and  oxygen  were 
investigated  by  E.  Greenbaum  (Oak  Ridge 
National  Lab).  He  screened  a  number  of 
marine  algae  and  found  several  that  could 
produce  hydrogen  after  a  period  of  adapta¬ 
tion  under  dark,  anaerobic,  and  free  C02 
incubation.  G.  Subramanian  and 
D.  Prabaharan  (India)  found  pure  cultures 
of  marine  blue-green  algae  that  produce 
hydrogen  from  seawater  upon  illumination. 

Macroalgae 

Work  on  propagation  of  seaweeds  and 
their  genetic  improvement  by  the  disasso- 
ciation  of  the  tissue  of  these  plants  to  cells 
and  protoplasts  is  proceeding  in  many  labo¬ 
ratories.  The  regeneration  of  plants  from 
the  isolated  protoplasts  occurs  readily  in 
some  cases  ( Porphyra )  but  is  still  a  problem 
for  other  algae  (. Ecklonia ). 

Zhang  et  al.  reported  on  formation  of 
hybrids  of  different  varieties  of  Gracilaria 
with  improved  properties  for  cultivation  at 
different  locations  in  China.  A.  Miura  and 
J.-A.  Shin  (Tokyo  Univ.  of  Fisheries)  reported 


on  hybridization  of  Porphyra  for  improved 
nori  production.  J.P.  Van  der  Meer  (Canada) 
reported  on  genetic  crosses  in  Gelidium. 

Experimental  modifications  in  the 
composition  of  cultivated  algae  were  studied. 
Polysaccharides  of  agarophytes  were  studied 
by  A.P.  KoUist  et  al.  (U.S.S.R.).  M.  Nishizawa 
(Defense  Academy)  and  T.  Hirano  (Tokyo 
Univ.  of  Fisheries)  attempted  to  cause  cross 
linking  of  the  algal  macromolecules  (agar) 
by  ionizing  radiation.  I  believe  that  they  are 
attempting  to  improve  the  quality  of  their 
biocolloids  by  the  cross  linking  treatment. 

The  quality  of  fatty  acids  of  the  brown 
alga  Laminaria  was  studied  by  M.  Honya  et 
al.  (Ishikawajima-Harima  Heavy  Industries). 
An  interesting  product  from  the  red  alga 
Corallina,  collected  from  nature,  was  bro- 
moperoxidase.  This  enzyme  is  interesting  as 
a  potential  catalyst  for  halogenation  of  sub¬ 
stances  in  industrial  bioreactors  (Y.  Izumi  et 
al.,  Kyoto  Univ.). 

T.  Kajiwara  et  al.  (Yamaguchi  Univ.) 
reported  on  the  cultivation  of  the  green  alga 
Ulva.  They  started  with  isolated  protoplasts 
and  grew  dense  unialgal  cultures  in  flasks. 
They  analyzed  the  culture  medium  and  the 
plants  for  the  production  of  volatile  sub¬ 
stances  that  are  responsible  for  the  typical 
flavor  and  smell  of  these  plants.  They 
reported  that  their  cultured  material 
appeared  to  be  enriched  with  some  specific 
components  when  compared  to  the  plants 
collected  from  nature. 

Marine  Animals 

Interesting  reports  on  the  genetic 
manipulations  of  cultivated  fish  were  pre¬ 
sented.  B.  Chevassus  reported  on  the 
research  in  France.  They  are  injecting  fish 
eggs  with  plasmids  that  cany  detectable  genes. 
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Their  expression  was,  however,  not  yet 
conclusively  demonstrated.  Chromosome 
manipulations,  by  mitotic  blocking  at  differ¬ 
ent  stages  of  zygote  development,  fertilizing 
UV  irradiated  eggs  or  using  UV  irradiated 
sperm  for  production  of  haploid  or  poly¬ 
ploid  fish,  are  apparently  successful. 

Transgenic  carp  and  loach  were  pro¬ 
duced  that  contained  and  expressed  the  genes 
for  rainbow  trout  growth  hormones  or  human 
growth  hormone.  These  genes  were  intro¬ 
duced  in  constructed  plasmids  and  were 
injected  into  the  recipient  eggs.  Over  half  of 
the  injected  eggs  developed  and  the  pres¬ 
ence  of  the  foreign  DNA  could  be  demon¬ 
strated  in  the  fish.  The  expression  of  these 
genes  could  also  be  detected.  A  significant 
fraction  of  the  Ft  progeny  from  crosses  of 
the  transgenic  fish  also  carried  the  foreign 
DNA.  The  transgenic  fish  grew  faster  than 
their  siblings  (T.T.  Chen  et  al.,  Univ.  of 
Maryland).  Biologically  active  tuna  growth 
hormone  was  produced  in  E.  coli  and  in 
yeast  by  recombinant  DNA  techniques 
(M.  Nonaka  et  al.,  Taiyo  Fishery,  Tokyo). 
Carp  alpha-globin  genes  were  introduced 
into  rainbow  trout  (T.  Oshiro  et.  al.,  Tokyo 
Univ.  of  Fisheries).  Direct  injection  of  cloned 
genes  into  the  nuclei  of  the  Medaka  eggs 
served  as  a  model  system  for  genetic  manip¬ 
ulations  of  fish  by  several  groups  (E.  Tamiya 
et  al.,  Univ.  of  Tokyo;  K.  Inoue  et  al.,  Nippon 
Suisan  Kaisha;  and  K.  Ozato  et  al.,  Kyoto 
Univ.). 

The  uptake  by  rainbow  trout  of  intes- 
tinally  administered  salmon  growth  hormone 
was  demonstrated  (S.  Moriyama  et  al., 
Kitasato  Univ.),  and  the  use  of  hormones 
for  improving  and  inducing  breeding  behav¬ 
ior  in  cultured  fish  was  reported  (J.M.  Little 
and  J.A.  Dawson,  Canada). 


M.  Maeda  (Univ.  of  Tokyo)  and 
H.  Sugita  et  al.  (Nihon  Univ.,  Tokyo)  report 
that  the  bacterial  flora  of  fish  intestines 
contribute  to  their  well-being,  which  sug¬ 
gests  the  possibility  of  introducing  geneti¬ 
cally  engineered  bacteria  with  desired  prop¬ 
erties  to  cultivated  fish. 

Monoclonal  antibodies  were  produced 
as  vaccines  to  immunize  cultured  fish  against 
viral  infection  and  also  for  diagnostic  pur¬ 
poses  (A.  Renard  et  al.,  Belgium;  Y.  Kamei 
et  al.,  Sapporo  Breweries  Co.,  Tokyo). 

Several  papers  reported  on  marine 
invertebrate  animals  as  sources  of  cytotoxic 
antifungal  and  antibiotic  substances 
(M.  Guyot,  CNRS,  France;  R.J.  Quinn  et 
al.,  Australia;  and  C.  Del-Rio  E.  et  al., 
Mexico). 

FINAL  REMARKS 

The  important  value  of  this  confer¬ 
ence  was  in  bringing  into  focus  the  state  of 
the  art  in  marine  biotechnology.  Most 
impressive  was  the  progress  reported  by  fish 
culturists  in  utilizing  the  developments  of 
molecular  biology.  Immunological  and 
recombinant  DNA  techniques  are  rapidly 
incorporated  into  this  technology.  The 
conference  should  also  be  viewed  in  the 
broader  context  of  the  increased  interest, 
especially  in  Japan,  in  marine  biotechnol¬ 
ogy.  The  Japanese  Government,  in  cooper¬ 
ation  with  major  industries,  is  now  construct¬ 
ing  two  major  research  centers  devoted  to 
the  development  of  marine  biotechnology, 
one  at  Shimizu  in  Shizuoka  Prefecture  and 
the  other  at  Kamaishi  in  Iwate  Prefecture. 
The  focus  of  the  research  will  be:  (1)  identi¬ 
fication  and  production  of  new  products 
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from  marine  organisms,  (2)  the  propaga¬ 
tion  of  useful  marine  organisms,  and  (3)  the 
control  of  the  marine  environment  from 
man-made  pollutants.  The  obvious  depen¬ 
dence  of  the  Japanese  economy  and  the  diet 
on  the  ocean  is  the  natural  driving  force 
behind  this  increased  interest  in  marine 
biotechnology.  I  foresee  rapid  developments 
in  marine  biotechnology,  especially  in  Japan, 
in  the  near  future. 


Aharon  Gibor  began  a  i  -year  assignment  at  the  Office  of  Naval  Research  Far  East  in  August  1 989.  Dr.  Gibor 
Is  a  professor  of  biology  at  the  University  of  California,  Santa  Barbara  and  is  now  on  a  leave  of  absence  from 
that  institution.  He  received  a  B.A.  degree  in  1950,  his  M.A.  degree  in  1952  from  the  University  of  California, 
Berkeley,  and  his  Ph.D.  degree  in  1 956  from  Stanford  University.  His  thesis  research  was  done  at  the  Hopkins 
Marine  Station.  Dr.  Gibor  was  involved  in  research  on  the  genetic  autonomy  of  cytoplasmic  organelles  of 
eukaryotic  cells,  especially  chloropiasts  and  flagella.  His  present  research  is  on  the  growth  and  development 
of  algal  cells  and  tissues  and  the  role  of  celt  walls  of  these  plants  In  controlling  their  development. 
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FIFTH  INTERNATIONAL 
CONFERENCE  ON  NUMERICAL 
SHIP  HYDRODYNAMICS 


Edwin  P.  Rood 


It  was  clear  from  the  papers  presented  at  this  conference  that  there  has  been  a 
world-wide  attempt  to  apply  supercomputers  to  prediction  of  ship  hydrodynamics, 
both  in  the  numerical  prediction  of  the  flow  and  in  the  processing  of  large  volumes 
of  experimentally  acquired  data.  In  this  sense  the  conference  was  addressed  toward 
resolving  a  technology  issue  of  solving  complex  hydrodynamics  with  brute  force 
computer  power. 


This  international  conference  was  the 
fifth  in  a  series  that  has  included  meetings  in 
the  United  States  and  in  France.  The  con¬ 
ference  was  sponsored  by  the  Shipbuilding 
Research  Association  of  Japan,  the  David 
Taylor  Research  Center,  the  Office  of  Naval 
Research,  and  the  Naval  Studies  Board. 
Attendance  at  this  meeting  in  Japan  was 
excellent:  almost  200  attendees  from  approx¬ 
imately  20  countries.  There  were  43  papers, 
with  4  invited,  dealing  with  all  aspects  of 
numerical  ship  hydrodynamics.  The  confer¬ 
ence  is  scheduled  to  be  held  again  in  1993  in 
the  United  States. 

There  were  several  topics  addressed 
in  the  sessions,  including  Navier-Stokes 
solvers,  turbulence  simulation,  viscous  flows, 
free  surface  waves,  effects  of  surface  ten¬ 
sion,  and  image  processing.  It  was  clear 
from  the  presented  papers  that  there  has 
been  a  world-wide  attempt  to  apply  super¬ 
computers  to  prediction  of  ship  hydrody¬ 
namics,  both  in  the  numerical  prediction  of 
the  flow  and  in  the  processing  of  large  volumes 
of  experimentally  acquired  data.  Although 
there  was  no  formally  stated  theme  for  the 
conference,  it  is  easy  to  imagine  that  the 


theme  could  have  been  “Where  are  we  with 
supercomputers?”  In  this  sense  the  confer¬ 
ence  was  addressed  toward  resolving  a  tech¬ 
nology  issue  of  solving  complex  hydrody¬ 
namics  with  brute  force  computer  power. 

Why  is  there  a  focused  effort  on 
numerical  ship  hydrodynamics?  This  seems 
at  odds  with  the  notion  that  naval  hydrody¬ 
namics  is  a  conservative  field  with  the  image 
of  being  mature.  Perhaps  this  notion  is  true 
in  the  context  of  the  design  of  ships  to  meet 
the  current  commercial  demands  for  power¬ 
ing  and  motions.  However,  the  technology 
does  not  permit  predictions  outside  of  the 
historical  design  parameter  space,  and  surely 
does  not  render  even  a  good  guess  at  the  real 
fluid  details  of  the  flow  around  marine  vehi¬ 
cles.  This  latter  need  is  in  response  to  recent 
developments  in  remote  sensing,  which  is 
anticipated  to  permit  world-wide  monitor¬ 
ing  of  maritime  traffic  in  the  foreseeable 
future.  The  advent  of  supercomputers  has 
had  a  measurable  effect  on  hydrodynamic 
predictions  that  has  nevertheless  been  insig¬ 
nificant  for  design  usage.  More  specifically, 
supercomputers  permit  calculations  for  much 
more  complex  geometries,  but  the  lack  of 
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physics  in  the  governing  equations  renders 
those  sophisticated  predictions  not  any  more 
useful  than  simpler  predictions.  Where  will 
the  state  of  the  art  be  in  the  future?  Super¬ 
computers  have  another  function,  and  that 
is  to  permit  the  large  data  volume  analysis 
from  image  processing  of  whole  flow  field 
measurements.  This  much  neeued  informa¬ 
tion  will  permit  identification  and  under¬ 
standing  of  the  physics  required  to  produce 
accurate  numerical  predictions,  also  obtained 
with  the  supercomputer. 

Chief  among  the  needed  physics  are 
understanding  and  prediction  of  the  momen¬ 
tum  exchanges  produced  by  three- 
dimensional  flow  separations,  a  predomi¬ 
nant  feature  of  real  ship  flows  about  hulls 
and  propellers,  and  understanding  and  pre¬ 
diction  of  turbulent  flows  at  free  surfaces,  a 
topical  area  about  which  there  is  very  little 
understanding.  The  phenomena  associated 
with  the  postbreaking  wave,  for  which  there 
is  no  identification  of  the  physical  mecha¬ 
nisms  producing  bubbles  and  vorticity,  also 
merit  scientific  study.  Basic  research  is  pro¬ 
ceeding,  chiefly  in  the  United  States  for  the 
Navy,  at  a  low  level  in  all  these  areas. 

This,  then,  is  the  technological  chal¬ 
lenge  of  the  future:  to  apply  newly  emerging 
scientific  information  to  the  supercomputer 
prediction  of  ship  hydrodynamics  in  an  effi¬ 
cient  manner.  The  field  is  ripe  for  break¬ 
throughs  in  the  scientific  understanding  of 
real  flows  because  of  new  measurement  and 
prediction  techniques  made  possible  by  the 
supercomputer.  These  breakthroughs  will 
provide  information  important  not  only  to 
ship  hydrodynamics  but  also  to  such  related 
fields  as  physical  oceanography. 

Unfortunately,  there  was  no  signifi¬ 
cant  new  science  presented  at  the  meeting. 
This  is  not  because  of  a  lack  of  interest  but 
because  of  the  urgency  to  apply  supercom¬ 
puters  to  brute  force  computations  in  the 


recent  past  and  because  of  a  lack  of  long¬ 
term  financial  commitments  to  approach 
the  science  issues  with  the  emerging  new 
capabilities  in  instrumentation.  Image  pro¬ 
cessing  is  the  single  most  important  capabil¬ 
ity  to  enter  the  scene.  It  is  still  in  its  develop¬ 
mental  stage  and  will  prove  to  be  a  great 
advantage,  especially  when  used  in  conjunc¬ 
tion  with  supercomputer  numerical  predic¬ 
tions  to  guide  exploration  of  the  flow. 

The  conference  is  a  hallmark  and 
archival  in  its  confirmation  that  science  is 
required  before  supercomputer  numerical 
prediction  capability  can  be  fully  utilized. 

Exemplicative  papers  in  support  of 
the  above  conclusions  were  presented  by 
K.  Mori  et  al.  (“Development  of  a  New 
Velocity  Measurement  System  by  Using 
Computerized  Flow  Visualization  and 
Numerical  Modelling”),  T.  Hino  (“Compu¬ 
tation  of  a  Free  Surface  Flow  Around  an 
Advancing  Ship  by  the  Navier-Stokes  Equa¬ 
tions”),  W.  Saric  (“Boundary-Layer  Stabil¬ 
ity  and  Transition”),  and  S.A.  Orszag  et  al. 
(“RNG  Modeling  Techniques  of  Complex 
Turbulent  Flows”).  The  following  provide 
brief  glimpses  at  the  significant  results  from 
these  papers. 

Mori  has  completed  a  preliminary 
endeavor  to  closely  couple  physical  and 
numerical  experiments  on  line.  He  makes 
full  flow  field  measurements  with  image 
processing  techniques,  uses  the  processed 
velocities  to  validate  numerical  predictions 
of  the  measured  flow  field,  then  uses  the 
numerical  procedure  to  predict  the  flow  in 
the  regions  difficult  to  measure.  This 
approach  is  a  significant  change  in  strategy 
and  heralds  the  future  for  predictions  of 
complex  flows.  The  essential  element  is  the 
reliance  on  numerical  hydrodynamics  with 
the  presumption  that  the  reliance  on  physi¬ 
cal  data  will  decrease  with  sophistication  in 
the  numerical  experiments. 
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Hino  consumed  some  70  hours  of 
supercomputer  time  to  calculate  the  non¬ 
linear  free  surface  flow  around  an  advanc¬ 
ing  ship,  and  even  then  the  geometry  was 
relatively  simple.  His  calculations  were  sig¬ 
nificant  in  including  the  defining  grid  for  the 
free  surface  position  as  part  of  the  solution. 
The  gross  features  of  the  flow  were  pre¬ 
dicted  adequately,  but  the  details  of  the 
wake  were  poor. 

Saric  summarized  recent  advances  in 
understanding  the  physical  mechanisms  of 
boundary  layer  transition.  It  is  now  recog¬ 
nized  that  the  flow  must  be  considered  as  an 
open  system  possibly  modelled  by  dynami¬ 
cal  systems  theory.  The  implication  is  that 
secondary  stability  analysis  and  subharmonic 
dynamics  may  characterize  the  transition 
process.  Understanding  of  this  process  is 


important  for  ship  flows  around  append¬ 
ages  and  other  geometries  of  relatively  low 
dimensional  magnitude. 

Orszag  presented  the  case  for  renor¬ 
malization  group  theory  as  a  tool  for  turbu¬ 
lence  modelling.  This  theory  provides  an 
analytical  method  to  eliminate  small  scales 
from  the  Navier-Stokes  equations,  thus  lead¬ 
ing  to  a  dynamically  consistent  description 
of  the  large  scales.  The  method  has  worked 
for  channel  flows,  cavities,  and  backward¬ 
facing  steps.  The  need  in  ship  hydrodynam¬ 
ics  is  for  a  model  of  three-dimensional  sepa¬ 
ration  and  for  a  turbulence  expression  for 
the  wake  near  the  free  surface. 

There  are  numerous  other  papers  of 
high  quality  from  this  conference;  the  above 
examples  show  the  relevance  and  timeliness 
of  the  conference  as  a  whole. 


Edwin  P.  Rood  is  a  Scientific  Officer  in  the  Office  of  Naval  Research  (ONR)  Fluid  Dynamics  Program.  He 
is  responsible  for  sponsored  research  in  the  areas  of  ship  hydrodynamics,  hydroacoustics,  and  wake/ocean 
interactions.  He  received  B.S.  (1968)  and  M.S.  (1 971)  degrees  in  physics  from  Virginia  Polytechnic  Institute 
and  from  American  University.  In  1 984,  following  a  long  period  in  the  applied  research  environment  at  the 
David  Taylor  Research  Center,  he  received  his  Ph.D.  in  fluid  mechanics  from  Catholic  University.  Prior  to 
joining  ONR,  Dr.  Rood  was  head  of  the  Special  Systems  Branch  at  the  David  Taylor  Research  Center. 
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FIFTH  INTERNATIONAL 
CONFERENCE  ON  NUMERICAL 
SHIP  HYDRODYNAMICS: 
VISCOUS-FLOW  PAPERS 


Frederick  Stern 


The  conference,  held  in  Hiroshima,  Japan,  focused  on  the  latest  developments  in 
numerical  ship  hydrodynamics,  including  the  topics  of  viscous  flow,  wavy-inviscid 
flow,  image  processing,  and  cavitation.  In  addition,  an  International  Towing 
Tank  Conference  (ITTC)  workshop  on  computational  fluid  dynamics  validation 
and  group  discussions  on  Rankine-source  methods,  Navier-Stokes  solvers,  and 
boundary-integral  methods  for  radiation-diffraction  problems  were  held.  This 
report  focuses  primarily  on  the  topic  of  viscous  flow  with  some  brief  comments 
on  the  topics  of  image  processing  and  cavitation  and  on  the  ITTC  workshop  and 
Navier-Stokes  solver  group  discussion. 


INTRODUCTION 

The  Fifth  International  Conference 
on  Numerical  Ship  Hydrodynamics  (ICNSH) 
was  held  in  Hiroshima,  Japan,  on  25-28 
September  1989.  Previous  conferences  were 
held  in  1985  (Washington,  DC);  1981  (Paris, 
France);  1977  (Berkeley,  CA);  and  1975 
(Gaithersburg,  MD).  Hiroshima  had  189 
participants  representing  18  different  nations. 
Out  of  these,  39  percent  were  from  outside 
Japan.  The  United  States  had  18  represen¬ 
tatives.  There  were  50  papers,  including  5 
keynote  lectures,  held  in  11  sessions  of  which 
all  but  3  were  double  format.  The  topics 
included  viscous  flow  (Reynolds-averaged 
Navier-Stokes  solvers,  large-eddy  simula¬ 
tion,  free-surface  effects,  and  unsteady  flow), 
wavy-inviscid  flow  (Neumann-Kelvin  for¬ 
mulations,  body  motions,  radiation- 
diffraction  problems,  and  Rankine-source 


methods),  image  processing,  and  cavitation. 
In  addition,  an  International  Towing  Tank 
Conference  (ITTC)  workshop  on  computa¬ 
tional  fluid  dynamics  validation  and  group 
discussions  on  Rankine-source  methods, 
Navier-Stokes  solvers,  and  boundary-integral 
methods  for  radiation-diffraction  problems 
were  held.  The  conference  was  very  well 
organized  with  a  balanced  social  program 
that  stimulated  technical  exchanges.  How¬ 
ever,  the  double  session  format,  which  was 
introduced  for  the  first  time,  restricted 
attendee  participation  in  the  various  topics. 
This  report  focuses  primarily  on  the  topic  of 
viscous  flow  with  some  brief  comments  on 
the  topics  of  image  processing  and  cavita¬ 
tion  and  on  the  ITTC  workshop  and  Navier- 
Stokes  solver  group  discussion.  Many,  but 
not  all,  of  the  viscous-flow  papers  will  be 
discussed. 
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VISCOUS  FLOW 

The  viscous-flow  papers  were  well 
representative  of  the  present  state-of-the- 
art  both  from  the  standpoint  of  accomplish¬ 
ments  and  dilemmas.  Considerably  more 
papers  were  presented  on  this  topic  than  in 
previous  conferences  which,  no  doubt,  is  an 
indication  of  the  shift  of  interest  of  ship 
hydrodynamics  researchers  towards  viscous- 
flow  methods.  The  status  of  the  viscous-flow 
methods  is  such  that  practical  ship  geom¬ 
etries  can  be  considered,  including  com¬ 
plexities  such  as  appendages  and  propel¬ 
lers.  Comparisons  with  experimental  data 
indicate  that  many  features  of  the  flow  are 
adequately  simulated;  however,  turbulence 
modeling  and  grid  generation  appear  to  be 
pacesetting  issues  with  regard  to  future 
developments.  Only  very  limited  work  is 
being  done  with  regard  to  unsteady  flow  and 
the  inclusion  of  free-surface  effects.  An 
important  general  observation  is  that  many 
of  the  authors  presented  their  methods  and 
computer  programs  as  if  they  were  com¬ 
pletely  unique  and  original  when,  in  fact, 
they  were  not.  Authors  should  provide  infor¬ 
mation  concerning  the  history  of  the  com¬ 
puter  programs  they  are  using  (i.e.,  original 
developer  and  previous  applications)  as  well 
as  the  specific  modifications  and  extensions 
that  they  may  have  made  for  their  applica¬ 
tion.  This  will  facilitate  the  assessment  of 
their  specific  advancements  and  difficulties. 
In  most  cases,  the  contributions  are,  no  doubt, 
significant  but  difficult  to  assess  due  to  the 
lack  of  such  information. 

In  general,  the  viscous-flow  methods 
presented  at  the  conference  can  be  conve¬ 
niently  divided  into  three  groups  depending 
on  the  type  of  velocity-pressure  approach 
used:  pseudo-compressibility  approach  (PC), 


uncoupled  momentum  and  continuity  equa¬ 
tions  approach  (UMQ,  and  coupled  momen¬ 
tum  and  continuity  equations  approach 
(CMC).  Although  finer  divisions  are  possi¬ 
ble,  e.g.,  with  regard  to  discretization  tech¬ 
niques,  turbulence  models,  grid-generation 
techniques,  and  other  numerical  procedures, 
these  are,  in  most  cases,  relatively  minor  and 
the  details  of  which  are  beyond  the  scope  of 
this  report  There  was  only  one  CMC  method, 
i.e.,  Hoekstra.  In  most  cases,  the  UMC 
methods  used  either  the  SIMPLE(R)  or 
MAC  (i.e.,  Poisson  equation)  algorithms.  In 
the  former  case,  either  finite-volume  (FV) 
or  finite-analytic  (FA)  discretization  tech¬ 
niques  were  used,  whereas  in  the  latter  cases 
the  overall  numerics  closely  followed  the 
standard  MAC  methods.  The  UMC- 
SIMPLE-FV  codes  closely  followed  the 
Patanker  model,  and  the  UMC-SIMPLE- 
FA  codes  closely  followed  the  Iowa  model. 

Most  of  the  papers  reported  results 
for  bare  bodies  with  applications  for  the 
Wigley,  Series  60  CB  =  .6,  HSVA,  or  SSPA 
hull  forms:  Kodama  (PC);  Hoekstra  (CMC); 
Zhu  et  al.  (UMC,  MAC);  Oh  et  al.  (UMC, 
SIMPLE,  FV);  Larsson  et  al.  (UMC, 
SIMPLE,  FA);  Masuko  and  Ogiwara  (UMC, 
SIMPLE,  FV);  Piquet  and  Visonneau  (UMC, 
PISO,  FA);  and  Tzabiras  and  Loukakis 
(UMC,  SIMPLE,  FV).  The  differences  in 
the  results  between  these  various  methods 
are  not  profound,  which  supports  the  pre¬ 
vious  status  statement.  Interestingly,  the 
full  transformation  used  by  Larsson  et  al. 
and  the  inclusion  of  the  near-wall  turbu¬ 
lence  model  by  Piquet  and  Visonneau  do 
not  lead  to  significant  improvements,  at  least 
for  the  present  applications.  In  the  latter 
case,  the  inclusion  of  the  near-wall  model 
slowed  convergence  such  that  the  SIMPLE 
algorithm  was  abandoned  in  favor  of  the 
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PISO  algorithm.  Tzabiras  and  Loukakis 
include  results  for  full-scale  Reynolds  num¬ 
bers. 

Fujii’s  keynote  lecture  provided  a  good 
overview  of  the  status  of  viscous-flow  com¬ 
putations  for  aerodynamic  applications. 
Although  the  magnitude  of  the  effort  is 
considerably  larger,  the  status  of  this  work 
appears  to  be  the  same  as  that  stated  pre¬ 
viously  for  ship  hydrodynamics.  Yang  et  al. 
(PC)  report  results  for  axisymmetric  bodies, 
including  investigations  of  angle  of  attack 
(separated  flo^)  and  propeller-hull  interac¬ 
tion.  The  propeller-hull  interaction  work 
closely  follows  the  earlier  work  of  Stern  and 
Kim.  Hino  (UMC,  MAC)  investigated  free- 
surface  effects  using  a  method  that  closely 
follows  the  MAC  and  SUMMAC  methods* 
however,  extensions  were  made  for  body- 
fitted  coordinates  and  high  Reynolds  num¬ 
bers  and  an  approximate  method  is  used 
whereby  the  grid  is  fixed  and  does  not  con¬ 
form  to  the  free  surface.  The  results  are 
promising  but  clearly  have  difficulty  in  accu¬ 
rately  resolving  the  boundary  layer  and  wake. 
Kinoshita  et  al.  (UMC,  MAC)  investigated 
oscillatory  flow  around  a  circular  cylinder. 
The  results  are  in  close  agreement  with 
available  experimental  data.  Stem  and  Kim 
(UMC,  SIMPLE,  FA)  significantly  extended 
their  previous  work  on  propeller-hull  inter¬ 
action  to  a  “complete”  viscous-solution 
method  for  marine  propellers.  Results  are 
presented  for  an  idealized  geometry  that 
exhibits  many  of  the  distinctive  features  of 
interest. 

Lastly,  with  regard  to  the  viscous-flow 
papers,  the  keynote  lectures  of  Saric  and 
Orszag  et  al.  were  valuable  in  providing  to 
the  ship  hydrodynamics  community  current 
information  on  boundary-layer  stability  and 
transition  and  RNG  and  spectral  methods, 


respectively.  In  the  former  case,  recent 
work  indicates  that  there  is  some  hope  for  a 
useful  stability  and  transition  theory,  whereas 
in  the  latter  case,  it  appears  that  such  methods 
may  find  application  in  ship  hydrodynamics. 
Doi  presented  results  for  turbulent  channel 
flow  using  a  large-eddy  simulation,  finite- 
difference  method.  The  finite-difference 
scheme  appears  to  attenuate  the  turbulence 
intensities  in  comparison  with  the  results 
from  other  methods. 

IMAGE  PROCESSING 
AND  CAVITATION 

The  two  papers  on  image  processing, 
Mori  and  Ninomiya  and  Fu  et  al.,  were 
interesting  and  generated  considerable  dis¬ 
cussion.  The  latter  approach  appears  very 
promising  for  quantitative  measurements. 
Kubota  et  al.  presented  an  innovative 
approach  to  the  problem  of  bubble  cavita¬ 
tion  in  which  both  the  fluid  and  bubbles  are 
treated  as  a  continuum  of  a  compressible 
viscous  fluid  with  a  large  density  variation. 
The  results  reproduce  some  of  the  observed 
trends.  Kinnas  and  Fine  presented  an  “alter¬ 
native”  linearized  supercavitating  foil  method 
that  allows  for  leading  edge,  face,  and  mid¬ 
chord  detachment  as  well  as  partial  cavita¬ 
tion.  The  method  is  applied  for  a  foil  and  in 
a  quasi-steady  manner  for  unsteady  propel¬ 
ler  cavitation. 

CONCLUDING  REMARKS 

Although  vigorous  discussions  were 
held  in  the  ITTC  workshop  concerning 
computational  fluid  dynamics  validation,  it 
appears  that  a  consensus  of  opinion  is  unob¬ 
tainable  in  the  near  future.  This  is  not 
surprising  in  consideration  of  the  difficulty 
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of  validation  even  for  experimentation  and 
the  fact  that  uniqueness  theorems  do  not 
exist  for  the  solutions  to  the  governing  equa¬ 
tions  of  computational  fluid  dynamics.  The 
Navier-Stokes  solver  group  discussion  was 
directed  towards  a  developers  and  users 
forum  on  issues  related  to  computer  pro¬ 
gram  development;  however,  as  discussed 
earlier,  in  order  to  successfully  hold  such  a 
forum,  the  various  methods  being  discussed 
must  be  properly  categorized. 

In  summary,  the  Fifth  ICNSH  was  a 
very  successful  conference  and  provided  a 
useful  exchange  of  technology. 
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TURBULENCE  MODELING 
IN  JAPAN 


Hideo  Yoshihara 


Significant  research  on  turbulence  modeling  for  Reynolds-averaged  Navier/Stokes 
equations  in  Japan  is  reviewed.  Anisotropic  k,e  models  and  large  eddy  simulation 
studies  are  covered. 


INTRODUCTION 

In  aerospace  applications  the 
Reynolds-averaged  Navier/Stokes  equations 
with  eddy  diffusivity  turbulence  models  are 
used  almost  exclusively.  The  simplest  of 
these  turbulence  models,  such  as  the  Cebeci/ 
Smith  or  Baldwin/Lomax  algebraic  mixing 
length  models,  have  yielded  satisfactory 
solutions  in  many  cases.  There  are,  how¬ 
ever,  many  key  flows  for  which  the  algebraic 
turbulence  models  have  been  seriously  defi¬ 
cient.  One  such  case  of  importance  is  the 
strong  shock/boundary  layer  interaction  with 
post-shock  separation  arising,  for  example, 
on  a  swept  wing  at  transonic  speeds.  Here 
the  failing  in  part  stems  from  the  assump¬ 
tion  of  local  equilibrium  for  the  turbulent 
transport  where  the  local  transport  is  taken 
independent  of  its  prior  history.  For  swept 
shock  separations  arising  on  a  swept  wing, 
the  turbulent  transport  becomes  strongly 
dependent  on  the  crossflow  with  the  stream- 
wise  Reynolds  stress  significantly  decreased 
by  the  crossflow.  For  swept-shock/boundary 
layer  interactions,  both  prior  history  as  well 
as  crossflow  effects  must  be  incorporated 
into  the  turbulence  model.  For  flows  with 
significant  secondary  flow  as  in  the  sepa¬ 
rated  region  of  a  backward-facing  step,  in  a 
rectangular  duct,  or  in  the  wing/fuselage 


juncture  region,  assumption  of  isotropy  of 
the  turbulence  transport  ceases  to  be  valid. 
There  are  also  flows  in  which  the  turbulent 
transport  is  locally  in  equilibrium  but  the 
mixing  length  cannot  be  suitably  defined  as 
in  a  confluent  boundary  layer  on  a  flapped 
wing.  In  some  cases  a  mixing  length  can  be 
defined,  but  its  determination  becomes 
unwieldy  as  in  an  unstructured  or  irregularly 
defined  grid.  In  these  cases  one  must  turn  to 
differential  equation  turbulence  models. 

In  Japan  there  are  several  efforts  that 
address  some  of  the  above  shortcomings.  In 
the  following  sections,  recent  efforts  on 
anisotropic  k,e  models  and  large  eddy  simu¬ 
lation  calculations  are  described  and  some 
miscellaneous  results  on  other  modeling 
studies  are  presented. 

TWO-EQUATION  k,e 
TURBULENCE  MODELS 

The  k,e  model  lies  in  the  next  level  of 
the  eddy  diffusivity  hierarchy  from  the  alge¬ 
braic  models.  The  characteristic  velocity 
and  length  scales  are  defined  by  partial  dif¬ 
ferential  equations  for  the  turbulent  energy 
(k)  and  dissipation  length  scale  (e).  For 
secondary  flows  in  the  separated  region  of  a 
backward-facing  step  and  in  a  corner  of  a 
rectangular  duct  or  along  the  wing/fuselage 
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juncture,  the  standard  isotropic  k,e  method 
is  inapplicable.  An  anisotropic  model  is 
required  as  developed  by  Yoshizawa  (Ref  1), 
which  was  described  earlier  (Ref  2).  This 
model  was  derived  using  Kraichnan’s  Direct 
Interaction  Approximation  assuming  the 
turbulence  to  be  a  dichotomy  of  large  and 
small  eddies  widely  separated  in  scale.  Such 
a  model  is  valid  only  for  large  Reynolds 
numbers  and  is  inapplicable  in  the  low 
Reynolds  number  region  close  to  the  wall 
where  important  practical  effects  of  turbu¬ 
lent  transport  occur.  A  Van  Driest-type 
damping  function  was  used  as  an  expedient 
to  bridge  the  near-wall  region. 

Nisizima  and  Yoshizawa  (Ref  3) 
applied  the  anisotropic  model  to  a  Couette 
flow  at  a  Reynolds  number  of  14,000  (based 
on  channel  half-width  and  wall  velocity)  and 
a  fully  developed  planar  channel  flow  at 
Reynolds  numbers  of  12,704  and  32,637 
(based  on  channel  half-width  and  centerline 
velocity).  Test/theory  match  in  these  cases 
was  mixed.  In  both  flows  the  agreement  of 
the  mean  velocity  distributions  across  the 
channel,  in  particular  the  velocity  profiles 
near  the  wall,  was  excellent.  Test/theory 
agreement  of  the  turbulent  intensities  and 
Reynolds  stress,  hence  the  anisotropy,  was 
reasonably  good  across  most  of  the  channel, 
but  discrepancies  occurred  near  the  wall. 
Test/theory  comparison  was  found  to  be 
qualitatively  the  same  for  the  turbulent  energy 
and  the  dissipation.  In  summary,  incorpora¬ 
tion  of  the  anisotropy  greatly  improved  the 
test/theory  agreement  relative  to  the  isotro¬ 
pic  model  for  both  the  Couette  and  planar 
channel  flows. 

Eddy  diffusivity  modeling  for  the  above 
planar  channel  flow  is  “one  dimensional.” 
The  square  duct  flow  is  a  “two-dimensional” 
modeling  problem  and  is  thus  significantly 
more  difficult.  In  Reference  4  Nisizima  and 
Yoshizawa,  using  the  above  anisotropic  k,e 


model,  calculated  the  fully  developed  square 
duct  flow  for  a  Reynolds  number  of  43,379 
based  on  the  bulk  velocity  and  the  duct 
width.  For  the  damping  function  a  product 
of  the  single-wall  damping  functions  was 
assumed.  The  calculations  were  carried  out 
by  a  finite  difference  method  for  the  Navier/ 
Stokes  equations  expressed  in  terms  of  a 
stream  function.  Calculations  were  com¬ 
pared  with  several  sets  of  experimental  data. 
Data  at  the  computed  Reynolds  number 
were  taken  from  Reference  5.  Other  data  at 
a  significantly  different  Reynolds  number 
were  used  only  for  qualitative  comparison 
purposes.  Data  were  presented  as  contour 
maps  in  a  cross-sectional  plane.  The  stream- 
wise  and  crossflow  mean  velocities,  mean 
crossflow  velocity  vectors  (to  show  second¬ 
ary  vortical  flows),  turbulent  intensities  and 
Reynolds  stress  (to  show  anisotropy),  and 
the  turbulent  energy  were  plotted.  In  gen¬ 
eral,  the  computed  contour  lines  for  the 
above  quantities  showed  the  same  topolog¬ 
ical  features  as  the  experiments,  but  there 
were  significant  quantitative  differences, 
particularly  in  the  corner  region  where  the 
calculations  are  unreliable.  The  authors 
also  suggested  that  a  partial  cause  for  the 
significant  test/theory  mismatch  was  due  to 
the  experimental  duct  flow  not  being  fully 
developed.  The  location  of  the  secondary 
flow  vortices  agreed  reasonably  with  exper¬ 
iments. 

IMPROVED  k,e  METHOD  OF 
MYONG  AND  KASAGI 

The  primary  contribution  of  Myong 
and  Kasagi  (Ref  6)  was  to  extend  the  aniso¬ 
tropic  k,e  model  of  the  type  developed  by 
Nisizima  and  Yoshizawa  to  cover  the  near¬ 
wall  region.  This  was  accomplished  in  two 
steps.  First,  the  standard  isotropic  model 
was  extended  into  the  near-wall  region  by 
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introducing  a  new  near-wall  characteristic 
length  scale.  The  latter  incorporated  the 
influence  of  viscosity  and  the  proper  near¬ 
wall  behavior  of  the  velocity  fluctuation 
components.  A  composite  dissipation  length 
scale  was  then  formed  as  the  sum  of  the 
near-wall  length  scale  and  the  outer  flow 
length  scale  from  the  standard  isotropic 
model.  The  sum  was  multiplied  by  a  Van 
Driest-type  damping  function.  The  above 
composite  length  scale  was  suggested  earlier 
by  Rotta  (Ref  7).  The  second  step  was  to 
incorporate  anisotropy  by  postulating  two 
quadratic  mean  velocity  gradient  terms  to 
be  added  to  the  Reynolds  stress  tensor 
equation.  The  first  term  was  essentially  that 
derived  by  Yoshizawa,  while  the  second 
resulted  from  the  anisotropic  diffusion  term 
in  the  turbulent  energy  equation.  The  latter 
term  has  an  essential  contribution  in  the 
near-wall  region  but  is  negligible  elsewhere. 

The  improved  anisotropic  k,e  method 
was  applied  to  the  fully  developed  planar 
channel  flow  at  Reynolds  numbers  (based 
on  centerline  velocity)  of  5,580, 27,600,  and 
100,000  using  the  finite  volume  time-marching 
method  of  Patankar.  With  the  exception  of 
the  turbulent  energy,  the  test/theory  match 
of  the  mean  velocity,  turbulence  intensities, 
normal  and  Reynolds  stresses,  and  dissipa¬ 
tion  was  excellent.  In  the  case  of  the  turbu¬ 
lent  energy,  there  was  a  nonnegligible  test/ 
theory  mismatch  in  the  peak  region  near  the 
wall.  It  was  suggested  that  the  cause  was  due 
to  an  erroneous  measurement  of  the  stream- 
wise  turbulence  intensity.  The  computed 
results  also  agreed  closely  with  the  direct 
simulation  calculations  of  Kim,  Moin,  and 
Moser  (Ref  8)  at  the  lowest  Reynolds  num¬ 
ber  of  5,580. 

Myong  and  Kobayashi  (Ref  9)  more 
recently  used  the  above  improved  k,e  method 
to  calculate  the  flow  in  a  finite  square  duct  at 
a  Reynolds  number  based  on  bulk  velocity 


and  duct  width  of 25,000.  The  product  of  the 
single-wall  damping  functions  was  used.  The 
calculated  results  were  compared  to  the 
experimental  results  of  Gessner  and  Emery 
(Ref  10).  The  addition  of  the  near-wall 
improvement  led  to  a  good  overall  test/ 
theory  match.  The  streamwise  and  lateral 
distributions  of  the  mean  streamwise  veloc¬ 
ity  component  showed  excellent  agreement 
The  agreement  of  the  mean  crossflow  veloc¬ 
ity  was,  however,  less  satisfactory,  while  the 
distributions  of  the  turbulent  energy  and  the 
Reynolds  stress  across  the  duct  compared 
reasonably  with  experiments.  There  were 
discrepancies  in  the  wall  shear  in  the  comer 
region,  but  this  should  not  be  surprising  in 
view  of  the  unproven  two-dimensional 
damping  function  that  was  assumed.  In 
summary,  the  results  from  the  improved  k,e 
method  of  Myong  and  Kasagj  yielded  greatly 
improved  test/theory  agreement  with  little 
added  complexity  over  the  standard  isotro¬ 
pic  model. 

LARGE  EDDY  SIMULATION  (LES) 

In  a  direct  simulation  of  turbulent  flow 
the  unsteady  Navier/Stokes  equations  are 
solved  using  a  sufficiently  refined  mesh  and 
time  step  to  resolve  the  turbulent  fluctua¬ 
tions.  To  reduce  the  space/time  resolution 
requirements,  in  LES  a  filter  is  applied  to 
the  Navier/Stokes  equations  to  screen  out 
the  high  frequency  end  of  the  fluctuations. 
The  well-known  difficulty  arising  by  the  fil¬ 
tering  is  that  a  subgrid  scale  viscosity  appears 
in  the  momentum  equations  to  render  the 
system  of  equations  underdetermined.  The 
key  problem  in  LES  is  the  modeling  of  the 
subgrid  scale  viscosity.  In  the  following,  the 
recent  research  of  Horiuti  (Ref  11)  and 
Miyake  and  Kajishima  (Ref  12)  for  a  planar 
channel  flow  is  described. 
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LES  Calculations  of  Horiuti 

The  subgrid  scale  viscosity  is  given  by 
the  sum  of  three  contributions,  namely,  the 
Leonard  term,  cross-stress  term,  and  the 
Reynolds  stress  term.  In  an  earlier  contribu¬ 
tion  (Ref  13)  Horiuti,  in  a  LES  calculation  of 
a  fully  developed  planar  channel  flow,  con¬ 
firmed  that  the  Leonard  and  cross-stress 
terms  were  essentially  equal  but  of  opposite 
sign,  so  that  the  previous  practice  of  neglect¬ 
ing  only  the  cross-stress  term  would  lead  to 
substantial  deterioration  of  the  result.  He 
further  showed  that  the  use  of  Bardina’s 
expression  for  the  cross-stress  and  Reynolds 
stress  term  led  to  improved  LES  results.  It 
also  restored  Galilean  invariance,  though 
this  is  not  a  factor  in  a  fully  developed  chan¬ 
nel  flow.  More  recently  Horiuti  (Ref  11) 
implemented  the  Bardina  model  in  a  finite 
difference  method.  By  suitably  formulating 
the  differencing  of  the  Leonard  and  cross¬ 
stress  terms,  the  anisotropic  form  of  the 
subgrid  scale  viscosity  due  to  Yoshizawa 
(Ref  14)  was  evolved  through  the  quadratic 
truncation  error  term.  Of  special  interest 
was  the  “side-by-side”  calculation  of  the 
planar  channel  flow  with  the  direct  simula¬ 
tion  and  LES  methods.  The  direct  simula¬ 
tion  result  was  then  filtered  and  compared 
to  the  LES  solution.  Correlation  of  the  two 
solutions  in  terms  of  the  streamwise  fluctu¬ 
ation  velocity  averaged  over  the  channel 
width  was  of  the  order  of  0.75  across  the 
channel.  Computing  time  with  the  new  dif¬ 
ference  procedure  was  half  of  that  of  the 
fully  spectral  method. 

LES  Calculations  of 
Miyake  and  Kajishima 

Calculations  were  carried  out  for  a 
fully  developed  planar  channel  flow  at  a 
Reynolds  number  of  4,940  based  on  bulk 


velocity  and  channel  height  (Ref  12).  The 
mesh  dimensions  were  256  (streamwise), 
117  (normal),  128  (spanwise)  for  a  total  of 
3.83  x  106  mesh  points.  (This  compares  to  a 
mesh  dimension  192  x  129  x  160  with  a  total 
number  of  points  of  4  x  106  points  in  the 
direct  simulation  calculations  of  Kim,  Moin, 
and  Moser  (Ref  8),  who  used  a  fully  spectral 
method.)  The  Smagorinsky  model  was  used 
for  the  subgrid  viscosity.  Finite  differencing 
was  employed  in  the  normal  direction.  Three 
cases  considered  were  the  steady  channel 
flow,  a  channel  flow  rotating  about  an  axis 
normal  to  the  walls,  and  a  steady  channel 
flow  with  sinusoidal  streamwise  area  varia¬ 
tions.  In  the  basic  channel  flow  good  agree¬ 
ment  was  achieved  with  the  direct  simula¬ 
tion  results  of  Mansour,  Kim,  and  Moin 
(Ref  15)  for  the  mean  distance  between 
streaks,  location  of  the  turbulent  energy 
peaks,  and  other  statistical  quantities  as 
flatness  and  skewness  factors.  Good  agree¬ 
ment  with  the  direct  simulation  results  con¬ 
firmed  effects  of  the  Smagorinsky  viscosity, 
required  for  stability,  to  be  negligible. 

Miyake  and  Kajishima  used  their  LES 
solution  to  assess  Launder’s  model  of  the 
redistribution  terms  and  Hanjalic’s  isotro¬ 
pic  model  for  the  dissipation  term  used  in 
the  Reynolds  stress  turbulence  model.  The 
result  was  that  the  above  models  were  rea¬ 
sonable  in  the  core  portion  of  the  duct,  but 
they  failed  badly  near  the  walls.  The  authors 
suggested  a  need  to  incorporate  some  aspects 
of  the  coherent  streak  structure  into  the 
modeling. 

In  the  rotating  channel,  the  effect  of 
the  Coriolis  force  was  correctly  reproduced 
with  Reynolds  stress  amplification  occur¬ 
ring  on  the  pressure  side  and  its  suppression 
on  the  suction  side.  At  a  sufficiently  large 
rotation  rate  a  near-relaminarization 
occurred  on  the  suction  side,  but  this  flow 
was  distorted  since  the  Smagorinsky  subgrid 
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viscosity  could  not  be  conveniently  turned 
off.  In  the  wavy  wall  calculation  a  sinusoidal 
porous  wall  condition  was  prescribed  on 
one  wall  Periodic  boundary  conditions  were 
retained  in  the  streamwise  direction.  The 
results  confirmed  that  the  turbulence  inten¬ 
sities  were  increased  in  the  decelerating 
region  and  decreased  in  the  accelerating 
region. 

MISCELLANEOUS  TURBULENCE 
MODEUNG  STUDIES 

Exploratory  Turbulence  Modeling 
Studies  on  Practical  Configurations 

Takakura,  Ogawa,  and  Ishiguro 
(Ref  16)  recently  reported  on  an  explora¬ 
tory  study  of  several  current  turbulence 
models  for  strong  shock/boundary  layer 
interaction  on  the  ONERA-M6  wing  where 
shock-induced  separation  was  present.  An 
ADI  Navier/Stokes  method  in  the  delta 
formulation  was  used,  and  the  turbulence 
models  investigated  were  the  widely  used 
algebraic  Baldwin/Lomax  model,  the  stan¬ 
dard  isotropic  k,e  model,  and  the  Smagorinsky 
subgrid  model  used  in  large  eddy  simula¬ 
tion.  The  use  of  the  latter  is  surprising  and 
inappropriate  since  the  subgrid  model  is 
intended  to  model  only  the  effects  of  the 
smallest  eddies,  so  that  significantly  finer 
mesh  and  shorter  time  step  must  be  used 
than  required  for  the  Reynolds-averaged 
simulation. 

On  the  M6  wing  the  shock  configura¬ 
tion  at  the  computed  Mach  number  and 
angle  of  attack  was  a  lambda  planform  shock 
composed  of  a  swept  forward  shock  and 
essentially  unswept  rear  and  outboard  shocks. 
Shock-induced  separation  occurred  down¬ 
stream  of  the  outboard  shock.  Calculated 
results  showed  a  significant  improvement  in 


the  test/theory  match  of  the  chordwise  pres¬ 
sure  distributions  with  the  k,e  model  rela¬ 
tive  to  the  Baldwin/Lomax  model,  particu¬ 
larly  in  the  important  outboard  shock  region 
where  separation  occurred.  The  k,e  model 
was,  however,  still  inadequate  for  the  shock/ 
boundary  layer  interaction,  producing  shock 
pressure-rises  that  were  too  large,  leading  to 
a  shock  location  that  was  too  far  aft.  The 
shortcoming  of  the  standard  k,e  method  is 
primarily  due  to  the  assumption  of  isotropy, 
where  significant  crossflow  effects  on  the 
turbulent  transport  are  omitted  in  the  inter¬ 
action  region.  Crossflow  here  is  essentially 
span  wise,  reducing  substantially  the  stream- 
wise  Reynolds  stress.  This  decrease  pro¬ 
duces  a  greater  wedging  displacement  effect 
of  the  boundary  layer  at  the  shock  which, 
then,  produces  the  needed  smaller  shock 
pressure-rise  and,  more  importantly,  the 
proper  upstream  displacement  of  the  shock. 

Another  illustrative  Navier/Stokes 
calculations  is  that  of  Sawada  (Ref  17)  on 
the  Japanese  space  shuttle  HOPE  configu¬ 
ration  at  M  =  3.5  and  10°  angle  of  attack. 
The  mesh  was  generated  by  a  multi-block 
system  using  an  algebraic  method.  With  this 
method  mesh  anomalies  arose  in  the  form 
of  “fictitious  corners”  at  the  vehicle  surface 
where  the  ‘normal”  coordinate  in  a  cross- 
sectional  plane  switched  coordinate  fami¬ 
lies.  This  complicated  the  determination  of 
the  mixing  length,  necessitating  use  of  a 
differential  equation  model,  in  this  case 
Coakley’s  q,w  two-equation  model.  Here 
the  algebraic  Baldwin/Lomax  turbulence 
would  have  sufficed  if  the  mesh  anomalies 
were  absent.  Thus,  if  a  more  appropriate 
mesh  generation  scheme  were  employed 
for  the  relatively  simple  shuttle  configura¬ 
tion  considered,  the  much  simpler  Baldwin/ 
Lomax  algebraic  turbulence  model  could 
have  been  used  as  was  done,  for  example,  by 
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Yamamoto,  Suzuki,  and  Mori  (Ref  18). 
Stiffness  effects  of  the  turbulence  equations 
arising  via  the  nonhomogeneous  terms  were 
circumvented  by  taking  a  sufficiently  small 
marching  step. 

Large-Scale  “Direct  Simulation” 

In  the  direct  simulation  of  turbulence 
the  mesh  must  be  sufficiently  refined  and 
the  time  step  sufficiently  short  that  the 
important  scales  of  turbulence  are  resolved. 
The  consequence  of  an  inadequate  mesh 
and  time  step  is  that  truncation  errors  will 
appear  that  distort  the  turbulence.  In  a 
number  of  papers  (e.g.,  Ref  19)  Kuwahara 
and  his  colleagues  assumed  that  the  large- 
scale  unsteady  portion  of  the  flow  as  the 
wake  of  a  bluff  body  could  be  simulated  by 
the  Navier/Stokes  equations  without  resolv¬ 
ing  the  small-scale  turbulence.  That  this  is 
erroneous  can  be  shown  by  considering  the 
simple  case  of  a  circular  cylinder.  Here  it  is 
known  that  the  location  of  the  separation 
point  and  the  rate  of  vorticity  shed  along  the 
separation  streamline  are  directly  depen¬ 
dent  on  the  nature  of  the  boundary  layer  just 
upstream  of  the  separation  point.  If  the 
mesh  and  time  step  are  inadequate,  the 
correct  character  of  the  boundary  layer 
upstream  of  the  separation  point  would  not 
be  properly  simulated.  As  a  result  the  sep¬ 
aration  point  and  vorticity  shed  would  be 
erroneous,  thus  distorting  the  large-scale 
wake  flow. 

SUMMARY  AND  CONCLUSIONS 

Unquestionably  outstanding  research 
on  turbulence  modeling  is  being  carried  out 
in  Japan.  Professors  Yoshizawa,  Miyake, 
and  Kasagi  and  Dr.  Horiuti  would  certainly 
rank  in  the  top  levels  of  competence  in  the 
United  States.  They  are  well  supported  by 


the  Ministry  of  Education,  Culture,  and 
Science  in  terms  of  research  staff  and  com¬ 
puter  support.  It  is  also  true  that  centers  of 
turbulence  modeling  are  few  in  number  in 
Japan,  at  least  an  order  of  magnitude  fewer 
than  in  the  United  States. 

The  basic  understanding  of  turbulence 
has  advanced  significantly,  particularly 
through  direct  simulation  calculations  as  those 
described  above.  Needs  of  the  applied 
computational  fluid  dynamicist,  however, 
are  not  much  closer  to  being  met.  Thus 
consider  the  relatively  simple  case  of  the 
comer  flow  where  “overlapping”  boundary 
layers  occur.  The  anisotropic  k,e  model  of 
Myong  and  Kasagi  greatly  improved  the 
prediction  of  the  secondary  flow  in  the  square 
duct.  In  a  more  severe  comer  flow  as  in  the 
wing/fuselage  juncture  in  which  a  strong 
separation  vortex  is  embedded,  there  will 
not  be  a  tolerance  to  a  comer  damping 
function  as  in  the  square  channel  case.  It  is 
not  clear  whether  renormalization  group 
theory  of  turbulence  might  not  provide  a 
viable  transport  model  free  of  unphysical 
comer  damping  functions. 

Turbulence  modeling  for  shock- 
induced  separation  downstream  of  a  swept 
shock  is  of  a  much  greater  complexity.  Exist¬ 
ing  modeling  approaches  as  those  described 
above  are  inadequate.  Modeling  the  depen¬ 
dence  of  turbulent  transport  on  the  cross¬ 
flow  in  this  situation  can  only  proceed  phe¬ 
nomenologically. 
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RECENT  JAPANESE  RESEARCH 
ON  OCEAN  WAVE-STRUCTURE 
INTERACTIONS 


Henry  T.  Wang 


Japanese  research  on  ocean  waves  and  their  interactions  with  offshore  structures 
ranges  over  a  broad  spectrum  of  activities  and  applications.  Many  of  the  research 
activities  aim  at  introducing  new  techniques  for  wave  measurement,  generation, 
and  modification.  Recent  research  efforts  on  ocean  wave-structure  interactions 
were  surveyed  during  site  visits  to  11  Japanese  marine  engineering  research 
institutes. 


INTRODUCTION 

During  the  last  2  weeks  of  September 
1989,  both  before  and  after  the  Fifth  Inter¬ 
national  Conference  on  Numerical  Ship 
Hydrodynamics  held  in  Hiroshima,  I  made  a 
series  of  visits  to  Japanese  marine  engineer¬ 
ing  research  institutes  to  survey  the  latest 
research  activities  on  ocean  waves  and  their 
interactions  with  offshore  structures.  This 
is,  of  course,  one  of  the  major  research  areas 
for  this  island  nation,  which  is  the  world’s 
leading  shipbuilder  and  much  of  whose  land 
is  located  within  60  miles  of  a  seashore.  In 
all,  I  made  a  total  of  10  on-site  visits  and  one 
off-site  visit,  listed  below  in  the  order  that 
they  were  made: 

1.  Institute  of  Industrial  Science  (IIS), 
Tokyo  University,  Roppongi  Campus. 
Host:  Prof.  T.  Kinoshita. 

2.  Dept,  of  Naval  Architecture,  Tokyo 
University  (TNU),  Hongo  Campus. 
Hosts:  Profs.  H.  Kajitani  and  H.  Kato. 


3.  Dept,  of  Naval  Architecture  and  Ocean 
Engineering,  Yokohama  National  Uni¬ 
versity  (YNU).  Host:  Prof.  T.  Hirayama. 

4.  Ishikawajima-Harima  Heavy  Industries 
Co.,  Ltd.  Research  Institute  (IHI), 
Yokohama.  Hosts:  Drs.  S.  Yamashita 
and  S.  Ogiwara. 

5.  Ship  Research  Institute  (SRI),  Ministry 
of  Transport,  Mitaka  (near  Tokyo). 
Host:  Mr.  S.  Kato. 

6.  Akishima  Laboratory,  Mitsui 
Engineering  &  Shipbuilding  Co.,  Ltd. 
(MES),  Akishima  (near  Tokyo).  Hosts: 
Drs.  K.  Shimada  and  M.  Kobayashi. 

7.  Dept,  of  Naval  Architecture  and  Ocean 
Engineering,  Osaka  University  (ONU), 
Suita  (near  Osaka).  Host:  Prof.  K.  Saito. 

8.  Dept,  of  Naval  Architecture,  Osaka 
Prefecture  University  (OPU),  Sakai 
(near  Osaka).  Hosts:  Profs.  Y.  Ikeda 
and  Y.  Himeno. 
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9.  Akashi  Ship  Model  Basin  Co.,  Ltd. 
(ASMB),  Akashi  (near  Osaka)  (off-site 
visit).  Host:  Dr.  M.  Takagi. 

10.  Tsuyazaki  Sea  Safety  Research  Labo¬ 
ratory,  Research  Institute  for  Applied 
Mechanics  (RIAM),  Kyushu  Univer¬ 
sity,  Tsuyazaki  (40  miles  north  of 
Fukuoka).  Host:  Prof.  K.  Kawatate. 

11.  Dept,  of  Naval  Architecture  and  Ocean 
Engineering,  Hiroshima  University 
(HIU),  Higashi-Hiroshima  (20  miles  east 
of  Hiroshima).  Host:  Prof.  M.  Takaki. 

The  above  institutes  include  national 
universities  (TNU,  YNU,  ONU),  research 
centers  affiliated  with  universities  (IIS, 
RIAM),  locally  supported  universities  (OPU, 
HIU),  industrial  research  laboratories  (IHI, 
MES,  ASMB),  and  a  government  labora¬ 
tory  (SRI).  It  is  not  by  coincidence  that  all  of 
the  above  sites  lie  close  to  stops  on  the  high 
speed  Shinkansen  east-west  rail  line  con¬ 
necting  Tokyo  to  Fukuoka.  This  is  the 
country’s  single  most  important  transporta¬ 
tion  line.  In  the  large  majority  of  cases,  the 
visit  consisted  of  a  guided  tour  of  the  wave¬ 
making  and  other  experimental  facilities, 
technical  discussions  with  the  host  and  other 
researchers  of  the  laboratory,  and  the  trans¬ 
mittal  of  reports  and  papers  describing  their 
recent  research.  At  my  request,  in  order  to 
obtain  information  on  the  most  recent 
research,  a  number  of  the  documents  were 
in  Japanese.  Typically,  research  results  are 
first  published  in  Japanese  and,  after  a 
number  of  refinements,  later  appear  in 
English  language  versions. 


BRIEF  DESCRIPTION  OF 
INSTITUTES 

I  will  not  attempt  here  to  give  a  detailed 
description  of  the  administrative,  organiza¬ 
tional,  or  physical  features  of  each  institute. 
Instead,  I  will  concentrate  on  some  general 
institutional  features  which  differ  signifi¬ 
cantly  from  American  practice.  I  will  also 
mention  noteworthy  features  or  recent 
developments  of  interest  for  particular  insti¬ 
tutes. 

In  the  case  of  the  universities,  it  is  of 
interest  to  note  that  the  school  year  runs 
from  early  April  to  the  end  of  February 
instead  of  from  September  to  May.  The 
position  of  department  chairman  is  appar¬ 
ently  rotated  on  a  yearly  basis  among  the 
senior  faculty  members.  Several  of  them 
lamented  the  fact  that  many  of  the  B.S.  and 
even  M.S.  graduates  entered  fields  not  related 
to  ocean  engineering,  such  as  banking  and 
insurance.  The  dominant  reason  seems  to 
be  substantially  higher  pay.  It  is  not  often 
that  American  students  will  go  through  a 
rigorous  engineering  curriculum  and  then 
seek  other  employment.  Also,  Japanese 
students  do  not  contact  the  companies 
directly,  but  instead  make  their  wishes  known 
to  the  department  chairman,  who  then 
“introduces”  them  to  appropriate  employers. 

The  organizational  structure  varies 
from  university  to  university.  The  senior 
faculty  members  at  the  institutes  associated 
with  Tokyo  University,  IIS  and  TNU,  enjoy 
a  degree  of  autonomy  not  apparent  else¬ 
where.  Each  senior  member  heads  a  “labo¬ 
ratory”  consisting  of  several  research  assis¬ 
tants  and  graduate  students.  At  OPU,  the 
department  is  divided  into  four  chairs,  each 
headed  by  a  full  professor.  At  RIAM,  the 
institute  is  divided  into  groups,  each  with 
one  full  professor,  one  associate  professor, 
and  several  research  associates  or  graduate 
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students.  The  organization  at  the  remaining 
universities  seems  to  be  less  structured  and 
close  to  that  at  American  schools.  The  naval 
architecture  department  at  YNU  offers  a 
special  program,  in  which  all  the  courses  are 
given  in  English,  to  foreign  students  pur¬ 
suing  doctoral  degrees. 

As  for  physical  location,  Prof.  Kinoshita 
pointed  out  that  the  Roppongi  district,  where 
IIS  is  located,  is  the  most  expensive  area  in 
Tokyo  and  hence  in  all  of  Japan.  As  a  result, 
there  are  only  two  very  small  experimental 
channels  located  in  the  basement  of  his 
building.  The  institute  does  own  the  exper¬ 
imental  facility  at  Chiba  City,  located  some 
25  miles  to  the  east,  across  the  Tokyo  Bay. 
Among  the  facilities  there  is  a  large  sea¬ 
keeping  and  maneuvering  basin  with  capa¬ 
bility  for  waves,  wind,  and  current.  In  1982, 
the  Faculty  of  Engineering  of  HIU  moved  to 
the  new  Saijo  Campus  in  Higashi-Hiroshima 
as  the  first  of  a  series  of  moves  to  have  a 
unified  campus  at  this  new  location.  RIAM 
is  located  in  a  beautiful  and  serene  ocean 
resort  area  at  the  northeastern  comer  of 
Kyushu  Island. 

The  three  industrial  research  centers, 
IHI,  MES,  and  ASMB,  represent  different 
ways  of  affiliation  with  the  parent  company. 
IHI  has  a  direct  affiliation,  MES  is  finan¬ 
cially  autonomous,  while  ASMB  represents 
a  joint  investment  by  Kawasaki  Heavy  Indus¬ 
tries,  Ltd.  and  Hitachi  Zosen  Corporation 
(HZC).  The  IHI  research  center  has  on-site 
dormitory  facilities  for  visiting  officials  from 
other  branches  of  the  company.  As  a  sign  of 
the  high  premium  on  land  in  this  country, 
some  of  the  buildings  near  the  entrance  of 
this  center  must  unfortunately  be  razed  (and 
relocated  inside  the  complex)  to  make  way 
for  a  highway.  Dr.  K.  Shimada  of  the  MES 
research  center  pointed  out  that  as  a  result 
of  its  financial  autonomy,  the  amount  of 
basic  research  work  has  decreased. 


The  government  laboratory,  SRI,  is 
composed  of  nine  technical  divisions  engaged 
in  research  on  ship  propulsion  and  dynam¬ 
ics,  structural  mechanics  and  materials  sci¬ 
ence,  power,  low  temperature  environments, 
and  nuclear  technology.  This  scope  of  work 
is  carried  out  by  a  staff  of  approximately  265. 
By  way  of  contrast,  a  similar  American  labo¬ 
ratory  would  have  a  staff  several  times  larger. 

EXPERIMENTAL  WAVEMAKING 
FACILITIES 

The  waves  are  generated  in  a  variety 
of  facilities  including  towing  basins,  wave 
tanks,  water  channels,  and  maneuvering 
basins.  The  first  three  types  are  character¬ 
ized  by  length-to-width  ratios  of  the  order  of 
ten,  while  the  last  type  is  more  squarish  with 
values  of  the  ratio  between  one  and  three. 
In  most  cases,  the  principal  wavemaking 
mechanisms  are  installed  in  the  towing  basins. 
Listed  below,  in  order  of  increasing  length, 
are  the  length,  width,  and  depth  dimensions 
in  meters  of  the  largest  towing  basin  at  each 
institute: 


1. 

OPU: 

70x3x  1.6 

2. 

RIAM: 

80  x  8  x  3 

3. 

HIU: 

80  x  8  x  3.5 

4. 

TNU: 

86  x  3.5  x  2.4 

5. 

ONU: 

100x7.8x4.4 

6. 

YNU: 

100  x  8x  3.5 

7. 

ASMB: 

200  x  13  x  6.5 

8. 

IHI: 

210  x  10  x  5 

9. 

MES: 

220  x  14  x  6.5 

10. 

SRI: 

400x18x8 

It  can  be  seen  that  the  basin  sizes  fall 
into  three  groups:  those  at  the  universities 
are  of  the  order  of  100  m,  the  industrial 
basins  are  approximately  twice  as  long,  while 
the  government  laboratory  has  the  longest 
basin.  In  all  of  these  basins  there  is  the 
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capability  of  generating  unidirectional  (or  WAVE  MEASUREMENT, 
long-crested)  waves  propagating  along  the  MODELING,  AND  MODIFICATION 
length  of  the  basin  composed  of  one  fre¬ 
quency  (regular  waves)  or  a  series  of  fre-  One  of  this  country’s  principal  inter- 

quencies  (irregular  waves).  These  waves  ests  in  the  utilization  of  ocean  space  is  the 

are  generated  by  a  single  unit  plunger  or  flap  use  of  offshore  platforms  to  serve  such  diverse 

moving  with  prescribed  motions.  In  some  of  purposes  as  factories,  storage  areas,  sea- 

these  basins  (RIAM,  YNU,  ASMB,  and  a  ports,  or  airports.  Thus,  under  the  leader- 

smaller  100-  x  5-  x  2.7-m  MES  basin),  the  ship  of  S.  Ando,  the  director  of  the  Ocean 

wavemaker  consists  of  a  series  of  snake  type  Engineering  Division,  SRI,  in  collaboration 

plungers  that  may  move  independently,  with  a  number  of  industrial  companies,  is 

resulting  in  the  generation  of  multidirec-  conducting  a  major  at-sea  experiment  involv- 

tional  (or  short-crested)  waves.  Until  recently,  ing  the  floating  platform  POSEIDON.  This 

these  waves  have  not  been  attempted  in  structure  is  34  m  long,  24  m  wide,  26  m  high, 

towing  basins  since  it  was  felt  that  the  two  submerged  5.5  m,  and  has  a  displacement  of 

sidewalls  would  limit  the  range  of  possible  531  tons.  It  has  12  columns  arranged  in 

wave  directions.  The  newest  wavemaker  of  three  rows  and  is  moored  by  six  slack  chains 

the  short-crested  type  is  the  one  at  YNU,  approximately  2  miles  north  of  Yura  Port  in 

and  it  is  composed  of  16  independent  units.  the  Japan  Sea.  The  experiment  began  in 

S.  Takezawa  and  T.  Hirayama  of  YNU  have  September  1986.  This  platform  is  meant  to 

made  detailed  studies  of  the  behavior  of  the  serve  as  the  prototype  of  an  actual  platform 

generated  waves  in  a  basin.  The  main  draw-  that  may  be  three  times  larger.  The  plat- 

backs  are  that  the  feasible  wave  directions  form  gathers  a  variety  of  data  including  the 

are  limited  to  ±  45°  measured  from  the  prin-  ocean  environment,  its  own  motions,  cable 

cipal  direction,  which  in  turn  is  now  con-  tensions,  and  structural  strains.  In  addition, 

strained  to  be  along  the  basin  length.  The  corrosion  behavior  of  the  12  columns  is 

principal  advantage  is  that  the  generated  observed,  where  1  column  is  made  of  con- 

random  wave  field  is  statistically  uniform  crete  and  the  remaining  11  are  made  of  steel 

along  most  of  the  basin.  with  different  paints  on  the  surface.  The 

On-site  maneuvering  and  seakeeping  platform  is  unmanned,  with  personnel  peri- 

basins  are  available  at  TNU,  IHI,  and  SRI  odically  going  onboard  to  collect  the  recorded 
with  respective  length  and  width  dimen-  data. 

sions  of  80  x  30  m,  70  x  30  m,  and  80  x  80  m.  The  data  are  recorded  in  two  modes. 

In  these  facilities,  irregular  long-crested  waves  The  regular  mode  is  to  record  data  for 

may  be  generated  in  two  perpendicular  direc-  34  minutes  at  intervals  of  0.5  second  every 

tions.  The  unusually  large  size  of  the  SRI  6  hours.  There  is  also  the  option  of  manually 

basin  is  due  to  the  fact  that  it  is  an  outdoor  starting  the  recording  to  collect  data  during 

pond.  M.  Kan  of  SRI  showed  films  of  capsiz-  severe  sea  conditions.  The  data  on  ocean 

ing  tests  of  a  remotely  controlled  ship  model  waves  are  obtained  by  twice  integrating 

traveling  prescribed  courses  in  the  presence  measured  vertical  accelerations  as  well  as  by 

ofthe  generated  waves.  This  project  is  being  ultrasonic  waveheight  meters  mounted  on 

done  in  cooperation  with  the  Federal  the  sea  floor,  which  is  at  a  depth  of  180  m. 

Republic  of  Germany.  The  data  from  both  methods  are  in  good 

agreement.  The  data  indicate  that  the 
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measured  energy  spectrum  is  rather  sharply 
peaked  and  is  in  better  agreement  with 
JONSWAP  (Joint  North  Sea  Wave  Project) 
than  with  that  proposed  by  ISSC  (Interna¬ 
tional  Ship  Structure  Congress).  Also,  it  is 
found  that  the  commonly  used  value  of  3.3 
for  y,  the  enhancement  factor  near  the  peak 
region,  is  correct  only  for  a  fully  developed 
sea.  The  data  are  also  analyzed  to  give 
information  on  wave  grouping,  directional 
distribution  of  the  wave  energy,  and  correla¬ 
tion  between  wave  height  and  period. 

On  a  smaller  scale,  from  the  idyllic 
seashore  town  of  Tsuyazaki,  K.  Kawatate  of 
RIAM  routinely  carries  out  at-sea  measure¬ 
ments  of  surface  as  well  as  subsurface  envi¬ 
ronmental  data  in  the  waters  surrounding 
Kyushu  Island.  First,  there  is  a  research 
tower  located  1.2  miles  offshore  that  can  be 
conveniently  reached  by  a  small  boat.  A 
larger  boat  is  available  to  deploy  or  retrieve 
data  gathering  moorings,  which  are  located 
further  offshore.  His  day-to-day  concerns 
include  frequent  communication  with  the 
local  fishermen  to  avoid  deploying  moor¬ 
ings  in  heavily  fished  areas,  cut  lines  due  to 
human  activity  or  shark  bites,  and  enlisting 
the  aid  of  the  local  authorities  in  locating 
lost  moorings.  Along  with  H.  Hongi  and 
A  Kaneko,  he  has  developed  a  self-governing 
profiling  system  that  consists  of  a  self-powered 
moving  sensor  unit  whose  vertical  motions 
are  controlled  by  a  microprocessor.  He 
calculates  the  dynamic  motions  of  his  moor¬ 
ings  by  using  a  method  that  considers  only 
the  angular  (or  transverse)  motions  of  the 
cable  system,  thus  bypassing  the  in-line  (or 
longitudinal)  motions,  which  usually  domi¬ 
nate  computer  time  requirements.  Using 
this  approach,  he  reports  a  reversal  in  the 
usual  computer  time  behavior,  i.e.,  the  cal¬ 
culation  time  for  nylon  rope  treated  as  being 
inextensible  is  only  one-sixth  of  the  extensi¬ 
ble  case. 


T.  Hirayama  of  YNU  is  investigating 
the  possibility  of  obtaining  real-time  esti¬ 
mates  of  sea  spectra  based  on  the  measured 
dynamic  motions  of  the  ship  itself.  He  feels 
that  onboard  wave  recorders  are  affected  by 
the  motions  of  the  ship,  overboard  buoys  are 
relatively  difficult  to  use,  and  visual  observa¬ 
tions  are  subjective  and  difficult  to  conduct 
at  night.  He  uses  Fast  Fourier  Transforms 
(FFT)  to  obtain  the  ship  motion  energy 
spectrum  from  the  measured  ship  motions. 
He  then  gets  the  sea  wave  spectrum  by 
dividing  the  motion  spectrum  by  the  square 
of  precalculated  motion  transfer  functions, 
which  are  stored  in  the  onboard  microcom¬ 
puter.  The  transfer  functions  are  calculated 
by  using  the  so-called  “New  Strip  Method,” 
which  better  accounts  for  hydrodynamic 
interaction  between  ship  sections  than  tra¬ 
ditional  Strip  Theory.  He  reports  that,  on 
the  whole,  there  is  reasonable  agreement  of 
his  calculated  spectra  with  those  obtained 
by  direct  wave  measurements  using  an 
onboard  Tucker  type  recorder  and  an  over¬ 
board  buoy. 

M.  Takagi  of  ASMB  and  HZC  is 
concerned  with  the  mean  water  level  of  the 
waves  generated  by  the  wavemaker  in  a 
towing  basin.  By  using  second-order  theory, 
he  shows  that  for  the  case  when  the  wave- 
maker  is  generating  two  wave  components 
of  different  frequencies,  the  mean  water 
level  is  composed  of  two  terms:  a  constant 
-at,  which  always  results  in  a  lowering  or 
“set-down”  of  the  water  level,  and  a  slowly 
varying  long-wave  sinusoidal  term  of 
amplitude  a2.  Thus,  for  those  cases  where 
at>ax,  there  locally  may  be  a  raising  or  “set¬ 
up”  of  the  water  level.  Thus,  the  traditional 
term  “set-down”  should  really  be  changed 
to  “group-induced  long-wave.” 

Traditionally,  passive  beach-type 
dampers  have  been  used  to  absorb  the  waves 
at  the  ends  of  a  wave  basin.  More  recently, 
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active  absorbers  are  being  considered.  For 
the  case  of  an  incident  wave  with  a  single 
frequency,  a  spring-dashpot-floating  body 
system,  with  natural  frequency  equal  to  the 
frequency  of  the  wave  and  damping  coeffi¬ 
cient  equal  to  the  damping  of  the  waves 
radiated  by  the  body,  is  a  highly  efficient 
wave  absorber.  S.  Naito  and  S.  Nakamura  of 
ONU  treat  the  case  of  irregular  waves 
composed  of  many  frequencies.  Here,  it  is 
necessary  to  satisfy  the  equivalent  of  the 
single-frequency  conditions  in  the  time 
domain.  The  basic  idea  is  to  continuously 
measure  the  elevation  of  the  incoming  wave 
at  a  point  in  front  of  the  body,  extrapolate 
this  elevation  to  the  actual  location  of  the 
wave  absorbing  body,  and  calculate  a  convo¬ 
lution  integral,  which  gives  the  required 
control  force  that  leads  to  motions  of  the 
body,  which  give  rise  to  maximum  energy 
absorption.  Unlike  the  single-frequency  case, 
complete  absorption  is  usually  not  possible. 

T.  Kinoshita  of  IIS  considers  a  num¬ 
ber  of  devices  to  modify  the  incoming  ocean 
waves,  to  either  mitigate  their  harmful  effects 
or  utilize  their  beneficial  effects.  In  the  case 
of  mitigation,  breakwaters  often  are  needed 
to  minimize  the  wave  effect  on  harbors  and 
offshore  structures.  Usually,  these  break¬ 
waters  are  set  with  their  lengths  perpen¬ 
dicular  to  the  incident  waves  to  maximize 
the  reflected  (and  hence  minimize  the  trans¬ 
mitted)  waves.  Unfortunately,  this  often 
places  large  loads  on  the  mooring  lines 
anchoring  the  breakwater.  He  experimen¬ 
tally  shows  that  by  using  a  “multi-body-type” 
of  floating  breakwater,  whereby  A-shaped 
flat  plates  are  placed  parallel  to  the  incom¬ 
ing  waves,  much  smaller  forces  are  experi¬ 
enced  by  the  mooring  lines  to  obtain  trans¬ 
mission  coefficients  equivalent  to  those  for 
the  conventional  type.  As  a  result  of  the 
lower  forces,  he  estimates  that  the  weight  of 
the  mooring  lines  is  reduced  by  a  factor  of  15 


compared  to  that  needed  for  the  conven¬ 
tional  type.  He  also  considers  hypothetical 
bodies,  as  represented  by  specified  singular¬ 
ity  distributions,  which  give  no  reflected  waves 
upstream  and  “focus”  the  incoming  single¬ 
frequency  wave  at  the  downstream  point 
x=L,  y=0  where  L  may  be  considered  to  be 
the  focal  length.  Typically,  wave  amplitudes 
at  the  focal  point  are  five  times  larger  than 
that  of  the  incoming  wave.  One  possible 
application  is  for  recreational  purposes. 

Along  with  H.  Maeda  and  K.  Masuda, 
T.  Kinoshita  also  investigates  the  design  of 
devices  to  extract  useful  energy  from  the 
ocean  waves.  They  estimate  that  perfect 
conversion  of  the  wave  power  along  the 
3,000-mile  coastline  of  Japan  would  supply 
about  one-half  of  its  electric  power  needs. 
There  are  presently  no  wave  energy  extrac¬ 
tion  devices  in  full-scale  operation.  The 
essential  idea  is  to  design  devices  that  min¬ 
imize  the  transmitted  and  reflected  wave 
field,  thus  maximizing  the  energy  that  is 
extracted.  Three  types  are  currently  under 
world-wide  consideration:  the  Salter  Duck, 
the  Articulated  Raft,  and  the  Oscillating 
Water  Column  (OWC).  The  first  device 
gives  the  desired  wave  field  through  design 
of  body  shape,  the  second  device  achieves 
the  result  by  means  of  body  deformation, 
while  the  last  device  concentrates  on  maxi¬ 
mizing  the  efficiency  of  the  energy  conver¬ 
sion  process.  The  authors  consider  in  detail 
the  OWC  type  where  the  ocean  wave  is 
trapped  in  a  column  and  its  energy  is  trans¬ 
mitted  by  an  air  column  to  a  turbine.  The 
authors  report  total  efficiencies  of  the  order 
of  10  to  20  percent. 

UNEARIZED  ANALYSES  OF  WAVE- 
STRUCTURE  INTERACTIONS 

The  increased  computing  speed  and 
memory  capabilities  of  modem  mainframe 
computers,  novel  uses  of  offshore  structures, 
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and  the  desire  to  minimize  needless  over- 
design  have  led  to  increasing  use  of  the 
nonlinear  analysis  approaches  described  in 
the  following  section.  Nevertheless,  linear¬ 
ized  analyses  continue  to  be  widely  used  due 
to  their  relative  simplicity  of  calculation, 
ease  of  interpretation,  as  well  as  generality 
of  application  of  the  calculated  results.  For 
example,  most  of  the  investigations  discussed 
in  the  preceding  section  make  explicit  use  of 
linearized  techniques  or  assumptions.  These 
analyses  are  usually  carried  out  in  the  fre¬ 
quency  domain. 

The  basic  building  block  of  most  linear¬ 
ized  analyses  is  the  Green  function,  which 
represents  a  singularity  that  satisfies  all  the 
boundary  conditions  of  the  problem  (includ¬ 
ing  the  troublesome  free  surface  condition) 
except  for  the  kinematic  condition  of  no 
flow  through  the  structure  surface,  which 
then  serves  to  determine  the  strengths  of  the 
singularities  distributed  on  this  surface.  The 
three  most  widely  used  Green  functions  are 
the  zero  speed  case  for  a  pulsating  stationary 
source,  the  wave  resistance  case  for  a  non¬ 
pulsating  translating  source,  and  the  gen¬ 
eral  radiation  case  for  a  pulsating  translat¬ 
ing  source.  The  current  status  is  that  robust 
and  efficient  computation  methods  exist  for 
the  zero  speed  case,  and  tractable  (though 
not  necessarily  both  efficient  and  robust) 
schemes  have  recently  been  developed  for 
the  wave  resistance  case.  Thus,  the  current 
major  effort  is  on  developing  computation 
methods  for  the  radiation  case. 

At  MES,  M.  Kobayashi  and  K.  Shimada 
have  incorporated  in  their  structural  dynamic 
response  program  “DREAMS”  a  calcula¬ 
tion  method  that  makes  a  direct  evaluation 
of  the  usual  double  integral  representation 
of  this  function.  A  series  of  transformations 
is  performed  to  convert  the  complex  inte¬ 
grals  to  a  more  tractable  real  form. 


H.  Iwashita  and  M.  Ohkusu  of  RLAM 
start  instead  from  a  unique  single  integral 
representation  of  the  function  and  evaluate 
the  resulting  integral  along  numerically  cal¬ 
culated  paths  of  steepest  descent  on  which 
the  oscillatory  behavior  of  the  integrand  is 
minimized.  They  report  that  accurate  evalu¬ 
ations  (for  given  values  of  the  frequency, 
and  locations  of  the  source  and  field  points) 
only  require  computer  times  of  the  order  of 
0.01  second  on  a  FACOM  M  780/20  main¬ 
frame  computer.  Their  work  offers  the 
hope  that  problems  involving  oscillating 
structures  with  forward  speed  or  in  the  pres¬ 
ence  of  a  current  may  be  solved  satisfac¬ 
torily  in  the  near  future.  In  fact,  Y.  Higo  and 
K.  Sato  of  HIU  use  a  preliminary  version  of 
this  method  to  calculate  the  added  resis¬ 
tance  of  a  surface  piercing  cylinder  advanc¬ 
ing  in  waves. 

S.  Yamashita  of  IHI  analyzes  the 
behavior  of  a  variety  of  structures  with  pro¬ 
truding  underwater  shapes.  It  has  been  well 
known  for  some  time  that  by  placing  a  suit¬ 
able  bulbous  bow  at  a  certain  distance  below 
the  free  surface,  the  resulting  body  experi¬ 
ences  no  exciting  forces  in  the  vertical  (heave) 
direction  at  a  certain  frequency.  While 
experimental  measurements  do  confirm  the 
existence  of  this  “wave-free”  frequency,  the 
region  is  sharply  tuned;  that  is,  the  exciting 
forces  increase  rapidly  on  either  side  of  this 
frequency.  By  changing  the  shape  of  the 
underwater  bulb  and  moving  it  further  below 
the  free  surface,  he  shows  that  the  resulting 
body  has  two  wave-free  frequencies  for  the 
heave  forces.  In  addition,  these  forces  are 
now  relatively  small  for  all  frequencies  larger 
than  the  lowest  wave-free  frequency.  More 
recently,  he  extends  this  idea  to  obtain  totally 
submerged  and  floating  bodies  that  are  wave- 
free  in  the  sway  direction,  which  is  the  hori¬ 
zontal  direction  perpendicular  to  the  inci¬ 
dent  wave.  The  basic  approach  is  to  trace 
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the  streamlines  of  velocity  potentials  that 
satisfy  the  wave-free  condition  of  no  radi¬ 
ating  waves  when  the  body  oscillates  in  sway. 
The  floating  body  cases  are  characterized  by 
a  sharp  narrowing  of  the  cross  section  as  the 
free  surface  is  approached.  Experimental 
measurements  of  the  sway  force  tend  to 
have  a  minimum  at  the  calculated  wave-free 
point  but  do  not  go  completely  to  zero.  As 
an  example  of  the  ease  and  generality  of 
application  of  the  results  obtained  from 
linearized  analyses,  he  gives  a  series  of  charts 
showing  the  wave  frequencies,  as  functions 
of  semi-submersible  parameters,  at  which 
the  heave  response  is  large.  These  widely 
used  structures  are  characterized  by  an  above 
water  deck  supported  by  vertical  columns 
that  terminate  in  underwater  pontoons. 

In  linearized  analyses,  the  amplitudes 
of  both  the  exciting  waves  and  body  motions 
are  assumed  to  be  small.  As  an  example  of 
the  relaxation  of  the  requirement  for  small 
body  motions  while  still  retaining  the  sim¬ 
plicity  of  the  frequency  domain  analysis, 
S.  Yamashita  calculates  up  to  third  order 
(or  harmonic  of  the  basic  frequency)  the 
forces  on  a  thin  body  heaving  with  large 
amplitude.  The  time  variation  of  the  body 
wetted  surface  as  a  function  of  heaving 
amplitude  is  now  considered  in  calculating 
the  exciting  forces.  As  may  be  expected,  th? 
agreement  between  calculated  and  mea¬ 
sured  forces  is  excellent  for  first  order,  rea¬ 
sonable  for  second  order,  and  only  fair  for 
third  order. 

NONUNEAR  ANALYSES  OF  WAVE- 
STRUCTURE  INTERACTIONS 

Nonlinear  analyses  are  needed  when 
the  motion  or  wave  amplitudes  are  large  or 
nonlinear  forces  such  as  the  viscous  drag  are 
significant.  These  analyses  are  usually  car¬ 
ried  out  in  the  time  domain.  Since  the  added 


mass  and  wave  damping  coefficients  are  not 
constants  but  functions  of  the  frequency  of 
motion,  the  equations  take  the  form  of 
integro-differential  equations.  The  integral 
involves  a  convolution  integration  from  the 
initial  time  t=0  to  the  present  time  t  of  the 
product  of  the  memory  (or  impulse  response) 
function  and  the  body  motion.  K.  Saito  of 
ONU  and  M.  Takagi  of  ASMB  and  HZC 
have  made  detailed  investigations  of  the 
analytic  evaluation  of  this  function,  its  deter¬ 
mination  from  experimental  tests,  and 
methods  of  solving  the  integro-differential 
equations  of  motion.  In  the  case  of  analytic 
evaluation,  the  memory  function  is  essen¬ 
tially  the  Fourier  transform  over  all  fre¬ 
quencies  of  the  wave  damping  coefficient. 
To  facilitate  this  evaluation,  they  derive 
asymptotic  formulas  for  this  coefficient  at 
high  frequencies.  In  the  case  of  solution  of 
the  equations  of  motion,  they  find  that  the 
use  of  the  convolution  integral  gives  calcu¬ 
lated  motions  that  are  consistently  in  good 
agreement  with  measured  results  for  the 
case  of  a  body  moored  by  cables  with  various 
load-deflection  curves.  When  an  attempt  is 
made  to  simplify  the  equations  by  assigning 
constant  values  to  the  added  mass  and 
damping  coefficients  corresponding  to  either 
the  natural  frequency  of  the  system  a>n  or  the 
peak  frequency  of  the  exciting  waves  oop 
(thus  avoiding  the  need  to  evaluate  the 
troublesome  convolution  integral),  the  results 
are  usually  poor  for  the  to  case.  They  are 
reasonable  for  the  con  case  when  the  mooring 
nonlinearities  are  relatively  weak,  i.e.,  when 
the  system  is  not  subject  to  sudden  impulse- 
type  loads  such  as  those  that  occur  when  a 
slack  rope  suddenly  becomes  taut. 

A  number  of  choices  are  available  for 
the  experimental  determination  of  the 
memory  function.  The  tests  may  be  con¬ 
ducted  as  forced  or  free  oscillations  of  the 
body,  the  integral  may  be  evaluated  in  the 
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time  or  frequency  domain,  and  the  integral 
may  be  evaluated  directly  or  in  correlation 
form  (Blackman-Tukey  method).  After 
systematically  carrying  out  the  above 
approaches,  they  conclude  that  the 
Blackman-Tukey  method  gives  the  most 
consistently  accurate  estimates  of  the  memory 
function  in  both  the  frequency  and  time 
domains. 

Y.  Higo  and  M.  Takaki  of  HIU  con¬ 
sider  the  forces  acting  on  two-dimensional 
and  axisymmetric  bodies  undergoing  large 
heave  motions.  Their  work  may  be  con¬ 
sidered  to  be  the  time  domain  equivalent  of 
S.  Yamashita’s  frequency  domain  approach. 
In  this  case,  the  memory  function  must  be 
updated  for  a  new  underwater  shape  at  each 
instant  of  time.  Their  approach  is  as  follows. 
The  flow  field  around  the  body  is  divided 
into  two  parts:  a  wave  field  caused  by  the 
body  motions  during  the  previous  time  steps 
(thus  accounting  for  the  past  motions)  and 
an  impulsive  (infinite  frequency)  motion  of 
the  instantaneous  underwater  body.  The 
wave  field  and  underwater  body  shape  then 
are  updated  before  proceeding  to  the  next 
time  step.  As  in  the  case  for  S.  Yamashita, 
the  agreement  of  their  calculated  results 
with  experimental  values  is  good  for  first 
order  and  progressively  worsens  with  increas¬ 
ing  order. 

In  most  of  the  preceding  analyses, 
potential  flow  methods  have  been  used  in 
that  the  effect  of  viscosity  has  been  either 
completely  neglected  or  simply  accounted 
by  adding  a  term  for  the  viscous  drag  force, 
which  is  a  nonlinear  function  of  the  fluid 
velocity  relative  to  the  body.  This  is  usually 
permissible  for  floating  bodies  near  the  free 
surface  since  the  wave  damping  is  larger 
than  the  viscous  drag.  However,  for  cases 
where  substantial  parts  of  the  body  are  located 
below  the  free  surface,  for  certain  body 


shapes,  or  modes  of  motion,  the  effect  of 
viscosity  can  no  longer  be  neglected.  In 
these  cases,  there  is  a  complex  interaction 
between  the  viscous  and  potential  flow 
components.  Thus,  the  viscosity-induced 
shedding  of  vortices  changes  the  pressure 
distribution  and  thus  the  added  mass  coeffi¬ 
cients  computed  from  potential  flow.  Simi¬ 
larly,  the  unsteady  flow  due  to  the  ocean 
waves  affects  the  development  of  the  bound¬ 
ary  layer  and  vortex  generation  and  thus  the 
viscous  drag  on  the  body.  A  key  dimension¬ 
less  parameter  governing  the  extent  of  this 
interaction  is  the  well  known  Keulegan- 
Carpenter  (KC)  number,  which  may  be 
viewed  as  the  ratio  of  the  diameter  of  the 
wave  particle  orbit  to  a  representative  body 
length.  The  effect  of  viscosity  increases  with 
increasing  KC  number.  Y.  Ikeda,  K.  Otsuka, 
and  N.  Tanaka  of  ONU  experimentally  and 
analytically  investigate  the  effect  of  KC 
number  on  the  added  mass  and  drag  coeffi¬ 
cients  of  horizontally  placed  cylinders  in  the 
presence  of  regular  and  irregular  waves. 
Experimentally,  they  find  that  at  the  rela¬ 
tively  low  KC  number  of  2,  the  measured 
added  mass  coefficients  are  only  one-half  of 
the  potential  flow  values.  The  KC  number 
has  less  of  an  effect  on  the  drag  coefficient. 
These  results  are  qualitatively  confirmed  by 
their  theoretical  model,  which  uses  discrete 
vortices.  The  vortices  are  placed  such  that 
the  circular  cylinder  shape  is  obtained  by 
using  the  well-known  potential  flow  “circle 
theorem”  and  the  free  surface  is  approxi¬ 
mated  as  a  rigid  wall  by  simple  reflection  of 
the  underwater  vortices.  The  location  and 
strength  of  the  shed  vortices  are  based  on  a 
simple  model  for  the  development  of  an 
oscillating  boundary  layer.  This  model  is 
very  useful  for  observing  the  effect  of  previ¬ 
ously  shed  vortices  or  “memory  effect”  on 
the  development  of  flow  around  the  cylinder. 
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The  principal  focus  of  the  hydrody¬ 
namic  analysis  of  floating  platforms  such  as 
the  POSEIDON  is  on  a  specific  and  seem¬ 
ingly  unimportant  aspect:  the  surge  and 
sway  motions  (the  horizontal  directions, 
respectively,  parallel  and  perpendicular  to 
the  incident  waves)  at  very  low  frequencies 
that  are  below  the  frequency  range  for  the 
principal  energy  content  of  the  incoming 
waves.  Thus,  the  forces  here  arise  due  to 
nonlinear  interactions  and  hence  are  small 
to  higher  order.  The  problem  arises  because 
of  the  low  natural  frequency  of  the  moored 
platform  due  to  the  lack  of  a  hydrodynamic 
restoring  force  in  the  horizontal  directions 
and  the  low  restoring  force  of  the  slack 
cables  mooring  the  platform.  This  fact, 
coupled  with  low  damping  for  these  modes 
of  motion,  may  lead  to  large  second-order, 
low-frequency  resonance  type  of  motions. 
As  a  result,  the  structure  and  mooring  cables 
must  be  carefully  designed  to  avoid  this 
phenomenon.  This  means  that  the  force 
coefficients  must  be  carefully  determined, 
not  only  as  a  function  of  wave  frequency  but 
also  of  wave  amplitude,  i.e.,  of  KC  number. 
The  nonlinearity  of  the  response  to  the 
exciting  random  waves  also  makes  more 
difficult  the  extrapolation  of  calculated  or 
measured  model  motions  (usually  of  the 
order  of  minutes  of  real  time)  to  extreme 
values  of  the  response  over  a  span  of  hours 
or  perhaps  days,  the  typical  duration  of  a 
storm.  The  combination  of  the  usefulness 
and  industrial  need  for  these  platforms  and 
the  complexity  of  the  analysis  of  the  low- 
frequency  horizontal  motions  has  made  this 
an  area  of  active  world-wide  research  for 
the  past  20  years. 

S.  Kato  and  S.  Ando  of  SRI  and 
T.  Kinoshita,  S.  Takase,  and  K.  Takaiwa  of 
IIS  are  actively  engaged  in  research  in  this 
area  with,  of  course,  specific  application  to 
the  POSEIDON.  There  are  basically  two 


steps  that  are  required  to  obtain  the  extreme 
response  statistics  to  expected  ocean  storm 
waves.  First,  the  short-term  response  behav¬ 
ior  may  be  obtained  in  one  of  two  ways.  One 
is  to  painstakingly  conduct  forced  oscilla¬ 
tion  tests  to  measure  the  force  coefficients 
as  a  function  of  wave  frequency  and  ampli¬ 
tude,  and  motion  frequency  and  amplitude, 
and  then  use  them  to  solve  the  integro- 
differential  equations  of  motion  to  obtain 
the  response  behavior  in  various  sea  states. 
The  second  alternative  is  to  directly  mea¬ 
sure,  in  a  model  basin,  the  platform  response 
to  these  sea  states.  When  the  response  is 
linear,  the  extrapolation  to  obtain  extreme 
responses  over  longer  times  is  relatively 
straightforward.  In  fact,  the  calculation 
procedure  is  very  similar  to  that  used  for  the 
exciting  waves  themselves.  When  the 
response  is  nonlinear,  this  is  no  longer  the 
case  and  the  calculation  procedure  is  con¬ 
siderably  more  complex.  Now,  both  the  First 
(linear)  and  second  order  (quadratic) 
response  spectra  must  be  obtained  from  the 
measured  or  calculated  responses.  The 
resulting  estimate  of  the  extreme  responses 
is  not  precise  and  involves  a  number  of 
approximations  concerning  the  interactions 
between  the  motions,  the  bandwidth  of  the 
response,  and  finite  series  representations 
of  the  probability  density  functions.  The 
present  researchers  show  that  inclusion  of 
the  statistical  interference  between  the  first 
and  second  order  responses  leads  to  extreme 
values  that  are  typically  10  percent  higher 
than  those  obtained  by  neglecting  this  inter¬ 
action. 

CONCLUDING  REMARKS 

Japanese  research  on  ocean  waves 
and  their  interactions  with  offshore  struc¬ 
tures  ranges  over  a  broad  spectrum  of  activ¬ 
ities  and  applications.  The  POSEIDON 
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floating  platform  and  Oscillating  Water 
Column  wave  energy  absorption  projects 
offer  the  possibility  of  better  future  utiliza¬ 
tion  of  the  offshore  waters  as  resources  for 
expanded  working  space  and  additional 
electric  energy,  both  of  which  are  in  short 
supply  in  this  country.  Many  of  the  research 
activities  aim  at  introducing  new  techniques 
for  wave  measurement,  generation,  and 
modification.  These  include  the  measure¬ 
ment  of  the  motions  of  the  ship  itself  to 
obtain  real  time  estimates  of  the  ambient 
sea  waves,  the  use  of  active  devices  to  dampen 
waves,  the  generation  of  short-crested  waves 
in  a  long  towing  basin,  and  the  use  of  break¬ 
waters  aligned  parallel  to  the  incoming  waves. 
While  the  use  of  the  more  accurate  time 
domain  nonlinear  analyses  is,  of  course, 
increasing,  traditional  linearized  analyses  in 
the  frequency  domain  continue  to  be  widely 
used  due  to  their  simplicity  of  calculation 
and  ease  of  interpretation.  Here,  the  recent 
development  of  a  promising  method  for  the 
calculation  of  the  radiation  Green  function 
for  a  translating  pulsating  source  offers  the 
hope  of  more  accurate  analyses  for  offshore 
structures  moving  with  forward  speed  or  in 
the  presence  of  a  current.  In  the  case  of  the 
more  difficult  and  complex  nonlinear  anal¬ 
yses,  research  activities  are  proceeding  along 
a  number  of  directions.  One  approach  is  to 
develop  more  efficient  analytical  and  numeri¬ 
cal  evaluations  of  the  memory  or  impulse 
response  function  as  well  as  to  identify  the 
optimum  experimental  techniques  for  obtain¬ 
ing  this  function.  Another  approach  is  to 


conduct  experimental  and  simplified  ana¬ 
lytical  calculations  on  standard  bodies  such 
as  circular  cylinders  to  obtain  better  physi¬ 
cal  understanding  of  the  complex  interac¬ 
tions  between  the  wave  field,  body  motions, 
and  vortex  generation  at  the  body  surface. 
For  more  complex  actual  structures  such  as 
the  POSEIDON,  the  approach  is  to  pains¬ 
takingly  measure  the  force  coefficients  over 
a  matrix  of  wave  and  motion  frequencies 
and  amplitudes  and  obtain  the  motion 
response  in  the  time  domain  by  solving 
integro-differential  equations  of  motion.  Or, 
the  motions  may  be  directly  measured  in 
model  tests  for  different  expected  sea  states. 
Finally,  the  techniques  of  probability  theory 
for  nonlinear  responses  are  being  refined  to 
obtain  better  estimates  of  expected  extreme 
responses  in  storms,  based  on  extrapolation 
of  the  above  calculated  or  measured  short¬ 
term  motions. 
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THE  PAX  COMPUTER  AND 
QCDPAX:  HISTORY  STATUS, 
AND  EVALUATION 


David  K.  Kahaner 


The  PAX  parallel  processing  project  developed  at  the  University  of  Tsukuba  by 
T.  Hoshino  is  summarized.  PAX  is  a  two-dimensional  array  of  processors  that 
operate  in  MIMD  form.  QCDPAX  has  a  fast  floating  point  unit  specialty 
programmed  for  QCD  calculations.  Peak  performance  is  14  GFLOPS. 


INTRODUCTION 

Modern  computing  for  the  solution  of 
scientific  and  engineering  problems  began 
in  the  1940s.  It  was  recognized  by  Von 
Neumann  and  others  that  simulation  or 
modelling  of  real  world  phenomena  could 
be  used  to  reduce  the  amount  of  physical 
experimentation  in  certain  situations.  When 
no  experimentation  was  possible  it  was  hoped 
that  models  would  allow  extrapolation  and 
prediction  of  new  results.  Weather  model¬ 
ling  provides  an  excellent  example.  Solving 
the  differential  equations  describing  the 
temporal  evolution  of  pressure,  tempera¬ 
ture,  and  associated  variables  was  one  of  the 
earliest  motivations  for  computer  model¬ 
ling.  Even  today  the  problem  is  largely 
unsolved  although  tremendous  progress  has 
been  made.  This  growing  desire  to  compute 
has  fueled  research  in  computers  and  com¬ 
putational  science.  There  are  many  papers 
detailing  the  need  for  ever  faster  hardware 
and  software.  While  some  of  the  original 
problems  can  now  be  solved,  science’s  appe¬ 
tite  for  computing  has  grown  at  least  as 
rapidly  as  the  technology  has  been  able  to 
provide  solutions.  There  is  no  reason  to 
think  that  this  trend  will  abate. 


Computers  were  originally  sequential, 
each  operation  occurring  in  its  own  turn. 
This  is  essentially  the  “Von  Neumann  model,” 
although  early  scientists  appreciated  that 
there  is  often  simultaneity  in  nature  that 
could  be  exploited.  A  substantial  amount  of 
simultaneous  processing  takes  place  in  all 
computers  today.  For  example,  input  and 
output  are  performed  asynchronously  with 
arithmetic  processing.  But  it  is  only  within 
the  last  10  to  15  years  that  parallel  or  nonse¬ 
quential  numerical  processing  has  been 
available. 

Here,  some  aspects  of  a  long-term 
parallel  processing  research  project  (PACS, 
PAX,  and  QCDPAX)  begun  in  1977  at  Kyoto 
University  and  Hitachi  Corporation’s  nuclear 
power  division  are  summarized.  The  initial 
name,  Processor  Array  for  Continuum 
Simulation  (PACS),  was  soon  changed  to 
Processor  Array  experiment  (PAX). 
QCDPAX  is  the  current  running  computer 
(see  below).  Unless  it  is  necessary  to  distin¬ 
guish  between  them,  we  refer  to  the  family 
as  PAX  computers.  Research  is  now  cen¬ 
tered  at  the  Institutes  of  Engineering 
Mechanics  and  Information  Sciences,  Uni¬ 
versity  of  Tsukuba,  Ibaraki  305  Japan.  The 
principal  investigator  is  Professor  Tsuiomu 
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Hoshino,  although  he  has  many  long-term  called  spatial  replication  (Ref  2),  is  today’s 
collaborators,  including  Professors  best  bet  for  substantial  performance  improve- 
Shirakawa,  Iwasaki,  Oyanagi,  and  others.  ments  and  is  a  major  part  of  current  research 

A  number  of  papers  (Ref  1-9)  detail-  in  high  performance  computing.  However, 

ing  this  work  are  listed  in  the  references.  very  few  scientific  computations  can  be 

This  report  is  the  result  of  analyzing  many  of  performed  completely  asynchronously.  That 

these,  several  visits  between  December  1989  is  generally  limited  to  searches  of  tables  for 
and  January  1990  to  Professor  Hoshino’s  state  variables,  calculation  of  material  prop- 

Tsukuba  laboratory,  and  an  examination  of  erties,  portions  of  ray  tracing  algorithms, 
some  programs  that  are  now  running  on  etc.  that  are  data  driven  and  hence  will  be 
QCDPAX.  different  for  each  process.  Thus  a  program 

The  PAX  idea  was  to  capitalize  on  the  may  need  to  determine  information  about  a 

obvious  parallelism  in  many  continuum  function  at  a  large  number  of  points  on  the 
problems,  the  rapidly  decreasing  costs  of  boundary  of  a  region  before  further  pro¬ 
computer  hardware,  along  with  its  associated  cessing.  Each  of  these  computations  may  be 
miniaturization.  These  ideas  have  been  done  asynchronously,  but  it  will  also  be 

pursued  by  other  researchers,  and  PAX  was  necessary  to  globally  synchronize  before 

in  many  ways  inspired  by  research  on  the  proceeding  from  these  results.  Similarly  in 

ILLIAC  IV.  The  latter  was  a  SIMD  com-  any  computation  involving  iteration  or  time 

puter  (single  instruction,  multiple  data)  in  integration,  a  global  synchronization  of  all 

which  a  control  unit  forced  a  large  number  parallel  tasks  must  be  performed.  However, 

of  identical  processors  to  perform  identical  between  these  points  the  calculations  are 

calculations  on  different  data.  ILLIAC  IV  independent  and  may  be  performed  asyn- 

was  seminal  for  the  subsequent  develop-  chronously.  Thus  a  good  model  for  real 
ment  of  parallel  computers,  but  it  was  never  problems  is  “quasi  MIMD”  (Ref  3). 

commercialized.  Even  today  it  is  extremely  Supercomputers  also  incorporate 

difficult  to  design  a  SIMD  computer  that  “pipelining”  within  their  arithmetic  units, 
can  cope  with  data  dependent  branching,  This  means  that  an  operation  such  as  multi- 

and  this  has  restricted  SIMD’s  applicability  ply  is  broken  down  into  a  number  of  steps, 

to  special  kinds  of  problems  such  as  image  where  each  can  proceed  concurrently  on 

processing.  Current  SIMD  machines,  for  distinct  operands.  This  provides  a  different 
example  the  ICL  DAP,  are  often  criticized  type  of  parallelism  and  has  successfully 
for  this  restriction.  generated  performance  improvements  of 

PAX  also  uses  several  hundred  iden-  perhaps  one  order  of  magnitude,  roughly 

tical  processors,  but  each  of  these  can  oper-  equal  to  the  number  of  sections  or  steps  in 

ate  independently.  This  is  usually  referred  arithmetic  operations.  Pipeline  computers 
to  as  MIMD  processing-multiple  instruc-  are  often  referred  to  as  “vector”  computers, 

tion,  multiple  data-although  the  idea  was  Sometimes  two  or  more  operations,  such  as 

known  before  PAX.  Supercomputers  such  +  and  *,  can  be  “chained”  together  to  pro¬ 

as  the  Cray  Y/MP  or  NEC  SX-3  utilize  very  vide  further  speedup.  If  there  are  indepen- 
powerful  asynchronous  processors,  but  typi-  dent  arithmetic  units  several  operations  can 
cally  fewer  than  16.  Careful  programming  be  done  simultaneously.  Breaking  arithmetic 
can  sometimes  lead  to  speedups  nearly  equal  operations  into  sections  requires  some 
to  the  number  of  processors.  Synchronous  additional  hardware  and  results  in  its  own 
or  asynchronous  multiprocessing,  sometimes  increased  overhead;  it  is  necessary  to  “latch” 


between  sections  to  synchronize.  New  hard¬ 
ware  devices  such  as  Josephson  junctions 
might  provide  a  means  of  increasing  the 
number  of  steps  without  introducing  too 
much  additional  latching  overhead.  There  is 
also  the  possibility  of  pipelining  from  mem¬ 
ory  so  that  operands  can  be  provided  to  the 
arithmetic  units  in  greater  number  per  unit 
time.  Currently  this  is  not  being  done, 
although  the  primary  restriction  seems  to  be 
economic.  Pipelining  is  an  essential  ingredi¬ 
ent  in  improving  computer  performance, 
but  additional  gains  to  be  obtained  from  it 
appear  to  be  limited.  Most  computer  scien¬ 
tists  think  of  spatial  replication  as  true  par¬ 
allel  processing  and  feel  that  it  has  the  most 
potential  for  high  performance. 

PAX 

The  Introduction  above  provides  a 
setting  for  describing  the  work  by  Hoshino 
and  collaborators.  The  physical  questions 
that  originally  motivated  Hoshino  were 
related  to  nuclear  reactor  models,  and  from 
these  he  became  interested  in  “action  through 
a  medium”  problems.  For  example,  fluid 
flow  equations  are  normally  derived  by 
describing  various  conservation  laws  on  a 
small  volume  of  material.  Thus  the  state  at 
point  P  is  mostly  determined  by  local  data 
around  P.  Similarly,  in  statics  problems,  such 
as  solving  the  Laplace  or  Poisson  equation, 
the  usual  approach  is  to  discretize  the  region 
and  compute  approximations  to  the  spatial 
derivatives  at  P  by  suitable  averages  of  data 
adjacent  to  P.  Finite  differences  in  two 
dimensions  take  five-point  or  nine-point 
averages;  finite  element  approximations  can 
be  described  in  the  same  way.  On  the  other 
hand  many  particle  problems  do  not  satisfy 
this  proximity  condition,  nor  do  many  impor¬ 
tant  numerical  methods,  such  as  Gaussian 
elimination  with  partial  pivoting. 


All  PAX  computers  have  shared  the 
same  basic  design  philosophy: 

•  A  high  performance  parallel  computer  to 
solve  specific  real  engineering  problems 
rather  than  general  problems. 

•  Multiple  processing  units  (PU)  arranged 
in  a  torus  or  two-dimensional  rectangular 
lattice  (nearest  neighbor  mesh,  NNM). 

•  Fast  connections  for  data  between  neigh¬ 
boring  PUs. 

•  Synchronization  and  broadcasting 
between  many,  or  all,  PUs  whenever 
necessary,  but  otherwise  independent 
MIMD  processing. 

Various  PAX  computers  have  been 
constructed;  they  are  named  according  to 
the  number  of  PUs:  PAX-9,  PAX-32,  PAX- 
128,  QCDP AX-240,  and  QCDPAX-488.  The 
major  differences  between  them  are  the 
number  of  processors,  the  specific  hard¬ 
ware  technology  used  in  each  PU,  the  con¬ 
nection  network,  the  memory  and,  in  the 
case  of  the  last  two,  the  application  software 
that  is  being  developed.  QCDP  AX-488  is 
still  under  construction  as  this  is  written.  It  is 
scheduled  for  completion  in  March  1990 
and  will  have  a  peak  performance  of  about 
14  GFLOPS. 

Hardware 

The  PAX  computer  system  consists  of 
an  array  of  PUs,  a  host  computer,  and  an 
interface  unit  (HP I)  to  communicate  between 
the  host  and  the  PUs.  In  QCDP  AX  each  PU 
contains  a  Motorola  MC68020  (CPU)  with 
a  clock  speed  of  25  MHz.  The  host  is  a  Sun 
3-260  workstation  with  a  color  graphics  dis¬ 
play.  In  QCDPAX-488  the  PUs  are  arranged 
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in  a  24  by  16  array;  QCDPAX-288  has  them 
arranged  in  a  24  by  12  array.  The  PUs  are 
arranged  toroidally.  By  that  we  mean  that 
the  eastern  end  of  each  row  is  also  adjacent 
to  the  western  end,  similarly  for  north  and 
south. 

One  particularly  interesting  design 
feature  of  all  PAX  computers  is  the  provi¬ 
sion  for  fast  floating  point  computation.  This 
is  an  integral  part  of  most  scientific  comput¬ 
ing  but  is  not  always  included  in  multipro¬ 
cessing  computers  of  a  size  comparable  to 
PAX.  For  example,  the  ICL  DAP  lacks 
hardware  for  floating  point.  (All  supercom¬ 
puters  provide  floating  point  hardware.) 

In  QCDPAX  each  PU  contains  a  vector 
processing  unit  (VPU)  with  floating  point 
hardware,  memory,  and  control  circuits.  For 
fast  floating  point  computation  Hoshino  and 
his  associates  decided  that  the  usual  floating 
point  processor,  Motorola  68881,  that  is 
often  provided  with  68020  chips  was  not  fast 
enough.  Instead,  they  selected  as  a  floating 
point  processor  (FPU)  the  LSI  L64133  that 
has  a  peak  arithmetic  performance  of 
33  MFLOPS.  This  chip  can  perform  both  a 
32-bit  floating  point  add  and  multiply  in  a 
single  clock  period.  They  have  incorporated 
the  L64133  into  the  VPU,  which  also  has 
2  MB  of  high  speed  (25-ns)  SRAM  for  data 
memory  (dm),  another  8  KB  of  similar 
memory  called  writable  control  storage  (wcs) 
for  the  microcode  (software)  to  support  the 
floating  point  operations,  and  1  KB  ROM 
for  a  look-up  table.  There  is  a  three-part 
pipeline  associated  with  instruction  fetch- 
decode-execute.  This  pipeline,  along  with 
the  two  simultaneous  arithmetic  operations, 
accounts  for  the  name  “vector  processor.” 
Floating  point  operations  occur  after  the 
CPU  sends  instructions  and  data  informa¬ 
tion  to  the  VPU,  which  then  operates  inde¬ 
pendently.  Computations  such  as  evalua¬ 
tion  of  the  elementary  functions,  random 
numbers,  and  some  matrix  computations 


have  been  programmed  via  microcode  into 
the  wcs.  There  are  no  wait  states  when  the 
CPU  addresses  the  wcs.  The  VPU  also 
supports  vector  operations  that  are  set  up  by 
the  CPU  but  run  by  the  wcs  and  the  FPU. 
Compared  to  the  68881,  the  FPU  is  difficult 
to  program  and  programming  problems  have 
been  a  significant  time  sink  for  the  small 
PAX  research  group.  Nevertheless,  it 
accounts  for  the  high  peak  performance  of 
QCDPAX. 

Communication  and  Memory 

In  addition  to  the  68020  CPU  and 
VPU,  each  PU  contains  communication 
hardware  and  several  different  kinds  of 
memory.  The  PU  communicates  rapidly  to 
four  PUs  on  the  east,  west,  north,  and  south 
by  communication  memories  (cm).  Each  of 
these  is  an  8-KB  two-port  RAM  and  thus 
can  be  accessed  by  the  PUs  that  share  it. 
Adjacent  PUs  thus  communicate  by  reading 
and  writing  into  the  same  memory  location. 
Synchronization  is  by  a  flag  register.  The 
PUs  also  communicate  rapidly  with  the  host 
by  another  8-KB  two-port  RAM  that  is  asso¬ 
ciated  with  the  HPI  and  is  mapped  onto  the 
address  space  of  the  host.  As  the  host  writes 
to  its  memory  it  also  sets  a  switch  indicating 
which  PU  is  to  get  the  data,  and  that  data 
item  then  becomes  available  at  the  appro¬ 
priate  PU.  Communication  between  a  PU 
and  the  host  works  in  the  same  way.  Also,  by 
overlapping  the  addresses,  i.e.,  having  all 
the  PUs  associated  with  the  same  memory 
locations  on  the  host,  the  latter  can  transmit 
identical  data  to  all  the  PUs.  This  mecha¬ 
nism  also  transmits  data  from  one  PU  to  all 
the  others.  In  PAX  descriptions  this  type  of 
communication  is  known  as  a  “ferry.”  Each 
PU  also  has  4  MB  of  DRAM  local  memory 
(lm);  this  is  used  to  store  the  PU  program 
and  whatever  local  data  it  requires. 
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Host 

The  host  computer  is  used  to  compile 
and  assemble  source  programs,  load  them 
into  the  array  of  PUs,  initiate  parallel  tasks, 
send  and  receive  data  from  the  PUs,  and 
output  the  results  to  a  user  via  files  and  the 
color  display.  Each  local  PU,  or  node,  runs 
whatever  program  is  loaded  into  it.  All  PUs 
are  normally  loaded  with  the  same  program, 
but  this  is  not  a  requirement.  (This  is  similar 
to  other  asynchronous  multiprocessors  having 
many  PUs.)  Different  data  at  each  node 
means  that  the  program  execution  path,  or 
thread,  is  likely  to  be  different  for  each  PU. 

Parallelism 

To  summarize,  PAX  has  four  levels  of 
parallelism.  First,  each  PU  can  compute 
independently.  Second,  the  CPU  can  com¬ 
pute  independently  of  the  VPU.  Third,  in 
the  VPU,  instruction  fetch,  decode,  and 
execute  can  occur  simultaneously.  Fourth, 
within  the  FPU  addition  and  multiplication 
can  occur  simultaneously.  A  user  should 
only  have  to  think  about  the  first  level.  This 
determines  the  algorithm  that  is  selected 
and  how  it  is  programmed.  System  software 
should  be  responsible  for  parallelizing  on 
the  lower  levels.  Howevet,  Hoshino  and  his 
colleagues  have  resorted  to  programming  at 
all  four  levels  in  order  to  generate  good 
performance  from  QCDPAX. 

Programming 

The  host  computer  on  QCDPAX  runs 
a  standard  version  of  UNIX.  A  user’s  job  is 
to  write  a  main  program  for  the  host  and 
programs  for  the  PUs.  Originally,  PAX 
programs  were  written  in  Fortran,  with  some 
extensions.  Now  C  is  used  for  QCDPAX. 
This  has  been  extended  to  allow  program¬ 
mers  to  communicate  between  nodes,  to  the 


host,  etc.  Complex  arithmetic  is  supported. 
Assembly  language  is  also  available  for  tak¬ 
ing  advantage  of  the  additional  levels  of 
parallelism  mentioned  above,  and  using  it 
seems  to  be  an  essential  part  of  getting 
optimal  performance. 

The  structure  of  programs  written  to 
run  on  a  PU  will  be  almost  entirely  under¬ 
standable  to  any  experienced  C  program¬ 
mer.  There  are  some  differences  due  to  the 
communication.  Thus  variables  can  be 
declared  as  “fast,”  meaning  that  they  are  to 
be  stored  in  the  dm  of  the  VPU;  “slow,” 
meaning  that  they  are  to  be  stored  in  1m  of 
the  PU;  “north,”  “east,”  “south,”  or  “west,” 
meaning  that  they  are  to  be  stored  in  the 
communication  memory  and  hence  avail¬ 
able  to  the  adjacent  processor;  or  “ferry,” 
meaning  that  they  are  to  be  stored  in  the 
memory  space  that  is  mapped  onto  the  hosts 
memory  and  hence  are  available  to  other 
PUs.  The  “vfor”  loop  forces  computations 
to  occur  in  the  FPU;  a  more  traditional  “for” 
loop  generates  code  that  executes  on  the 
CPU. 

vfor  (k-0;  k<NMAXAL;  k+-l)  { 
u[k]  -  a[k]*la[k] ; 
v[k]  -  O.O5 
> 

for  (j-0;  j<10;  j+-l) 
l_ferry  -  l[j]; 

The  “vfor”  loop  microcode  instructions  and 
data  are  sent  to  the  VPU,  where  they  begin 
executing.  In  this  example  both  loops  exe¬ 
cute  simultaneously.  Usually,  a  “vfor”  loop 
is  inside  a  “for”  loop.  In  that  case  the  same 
microcode  is  sent  to  the  VPU  during  each 
iteration  of  the  “for,”  resulting  in  substantial 
overhead.  Hoshino  is  working  to  have  this 
transmission  only  occur  once. 

Hoshino  always  has  felt  that  it  was  a 
mistake  to  try  and  develop  very  general 
parallel  computers.  At  every  opportunity  he 


ONRFE  SCI  INFO  BUL  15  (2)  90 


61 


emphasizes  “simple  is  better.”  A  good  exam¬ 
ple  is  the  two-dimensional  lattice  arrange¬ 
ment  (NNM)  of  the  PUs  in  PAX.  In  going 
from  one  to  two  to  three  dimensions  the 
mean  distance  between  PUs  goes  from  P/2 
to  sqrt(P)  to  1.5*cbrt(P)  for  P  total  pro¬ 
cessors.  If  P = 40%  these  values  are  2048, 64, 
and  24.  Thus,  the  incremental  improvement 
in  three  dimensions  is  modest.  Another  well 
known  parallel  computer  architecture  is  a 
“hypercube,”  in  which  2**n  PUs  are  arranged 
at  the  corners  of  an  n-dimensional  lattice. 
Rapid  communication  is  possible  with  PUs 
sharing  an  edge.  Communication  with  PUs 
that  are  further  away  takes  longer  because, 
on  average,  n/2  edges  must  be  traversed. 
This  is  also  an  example  of  a  nearest  neighbor 
architecture.  PAX  is  a  two-dimensional 
specialization,  but  Hoshino  believes  that 
the  more  complicated  hardware  design  of  a 
hypercube  reduces  its  efficiency.  For  exam¬ 
ple,  the  ratio  of  maximum  communication 
overhead  for  a  4096- node  hypercube  (2**12) 
to  a  40%-node  PAX  computer  (2**6  by 
2**6)  is  128  to  12,  about  a  factor  of  10.  But 
if  overall  communication  overhead  is 
15  percent,  reducing  this,  even  to  zero,  does 
not  contribute  much  to  program  speedup.  It 
can  also  be  easier  to  design,  write,  and  debug 
programs  for  a  “two-dimensional”  computer 
compared  to  a  higher  dimensional  one,  but 
this  depends  upon  the  problem  and  the 
tools  that  are  available  for  the  user.  On  the 
other  hand  communication  time  has  been  a 
serious  problem  for  many  multiprocessing 
computers,  and  in  those  cases  reducing 
communication  costs  by  a  factor  of  10  would 
be  extremely  important.  Further,  some 
problems  map  better  onto  a  higher  dimen¬ 
sional  arrangement.  But  Hoshino  and  other 
PAX  researchers  have  shown  in  a  number  of 
papers  that  often  the  improvement  is  asymp¬ 
totic  and  it  may  be  difficult  to  achieve  the 
best  rate  in  practice. 


Performance 

Hoshino  has  repeatedly  stated  that 
any  parallel  processing  research  must  be 
verified  by  actual  implementation.  Although 
some  of  his  papers  estimate  the  efficiency  of 
PAX  by  mathematical  analysis,  he  has  run 
and  measured  performance  on  calculations 
including 

•  Poisson  equation  in  two  and  three  dimen¬ 
sions  via  successive  overrelaxation  (SOR), 
alternating  direction  implicit  (ADI),  and 
Fourier  analysis  cyclic  reduction  (FACR) 
methods 

•  Navier-Stokes  equations  via  Beam- 
Warming,  MacCormack,  and  SOLA 
methods 

•  Linear  equations  via  Gauss-Jordan  and 
conjugate  gradient  methods 

•  Molecular  dynamics 

•  Plasma  simulation  via  particle  in  cell  (PIC) 
method 

•  Two-dimensional  Heisenberg  spin  model 

•  Unsteady  flow  in  flooding  river 

•  Logic  circuit  analysis 

•  Quantum  chromodynamics  (QCD) 

Peak  measured  performance  is  8.25  GFLOPS 
on  the  288-node  QCDPAX.  Expected  peak 
performance  will  be  13.75  GFLOPS  on  the 
488-node  QCDPAX  in  March  1990.  The 
table  below,  from  Reference  6,  indicates  the 
measured  performance  of  QCDP AX-240 
on  the  solution  of  a  three-dimensional  Poisson 
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equation  using  a  seven-point  discretization 
stencil  and  red/black  SOR  iteration.  The 
mesh  is  360  x  370  x  360. 


Update  time  per  iteration,  fis . 199,294 

Communication  time,  /is . 265,214 

Efficiency,  % . 42.9 

GFLOPS . 0.98 

Time  to  converge,  s  . 226.7 


The  figure  of  42.9  percent  efficiency  given 
by  Hoshino  is  computed  as  computation 
time  divided  by  total  time.  The  latter  includes 
all  the  overhead  associated  with  noncompu- 
tational  tasks,  primarily  communication.  The 
usual  measure  is  the  ratio  of  measured 
speedup  relative  to  maximum  theoretical 
speedup.  Speedup  is  the  ratio  of  time  for 
serial  processing  to  solve  the  problem  to 
time  for  parallel  processing. 

QCDPAX 

Computers  with  parallel  processing 
capabilities  that  also  require  the  user  to 
formulate  an  algorithm  in  order  to  take 
advantage  of  parallelism  are  inherently 
special  purpose.  Some,  such  as  the  NASA 
Finite  Element  Machine  (FEM)  or  Massively 
Parallel  Processor  (MPP),  have  been 
designed  from  the  ground  up  for  a  particular 
class  of  problems;  in  the  case  of  MPP  this  is 
satellite  image  processing.  Others,  such  as 
the  Cray  YM/P-4,  can  be  used  for  general 
computation  without  making  use  of  parallel 
constructs.  PAX,  like  a  hypercube,  is  gen¬ 
eral  purpose  in  the  sense  that  the  hardware 
and  software  are  not  restrictive.  Both  are 
special  purpose  in  that  some  problems  are 
much  better  suited  to  their  architecture  than 
others.  They  are  also  special  purpose  in  that 
their  power  is  from  the  large  number  of 
modestly  powerful  processors,  rather  than  a 
modest  number  of  very  powerful  ones.  Within 


the  context  of  this  special  purpose  environ¬ 
ment,  the  PAX  researchers  have  always  tried 
to  develop  a  general  purpose  computer. 
The  current  QCDPAX  is  unusual  because  it 
alters  this  trend  and  was  developed  in 
response  to  funding  pressures  from  the 
Ministry  of  Education,  Science,  and  Culture. 

QCDPAX  joins  together  most  of  the 
original  PAX  researchers  with  physicists 
interested  in  a  specific  problem,  quantum 
chromodynamics  (QCD),  and  its  discrete 
model,  lattice  gauge  theory.  QCD  calcula¬ 
tions  involve  massive  amounts  of  random 
number  generation,  multiplication  by  3  by  3 
matrices,  and  computation  of  exponentials- 
all  floating  point  operations.  Thousands  of 
hours  of  supercomputer  time  are  required 
to  generate  interesting  results.  The  compu¬ 
tations  are  highly  parallel  and  satisfy  the 
proximity  property,  and  thus  are  extremely 
well  suited  to  PAX  architecture.  Special 
efforts  have  been  made  to  speed  up  the 
software  that  is  necessary  for  QCD.  For 
example,  3  by  3  matrix  multiplication  has 
been  written  in  microcode  for  the  FPU,  and 
complex  floating  point  is  supported  in  C. 
Several  other  research  projects  are  attempt¬ 
ing  to  develop  special  purpose  computers 
for  QCD  calculations,  notably  at  Columbia, 
CalTech,  IBM,  Rome,  Edinburgh,  and  at 
FermiLab.  A  summary  discussion  of  the 
progress  of  these  groups  through  1987  is 
given  by  Christ  (Ref  1).  In  this  report  we 
have  been  considering  PAX  as  a  research 
project  in  parallel  processing. 

SUMMARY  AND  EVALUATION 

Versions  of  PAX  have  been  under 
development  since  1977.  They  do  not  repre¬ 
sent  any  new  architectural  breakthroughs. 
Instead  PAX  is  a  determined  attempt  to 
design  a  parallel  computer  that  is  simply 
built  to  work  reliably,  simply  organized  to 


ONRFE  SCI  INFO  BUL  15  (2)  90 


63 


reduce  the  chore  of  programming,  and  not  many  systems  programming  headaches.  The 

too  special  purpose.  The  application  areas  PAX  research  group  is  very  small.  Writing 

are  those  continuum  problems  satisfying  the  microcode  is  time  consuming  and  does  not 

proximity  property  and  for  which  substan-  produce  new  physics.  The  C  language  exten- 

tial  amounts  of  rapid  floating  point  compu-  sions  that  enable  users  access  to  PAX’s 

tation  are  required.  QCDPAX  is  the  most  parallelism  are  easy  to  understand  but  their 

special  purpose  version  of  PAX,  but  only  in  implementation  still  has  bugs.  Some  of  these 

the  application  programs  and  microcode  have  to  be  worked  around  via  assembly 

that  are  specialized  for  QCD  calculations,  language.  Documentation  at  all  levels  is  weak, 

not  in  the  hardware.  The  lessons  learned  What  exists  is  entirely  in  Japanese.  This  is 

building  and  running  QCDPAX  can  be  used  not  as  much  of  a  disadvantage  as  it  seems 

directly  in  other  versions,  if  these  are  built.  because  program  segments  are  the  best 

Although  many  computations  have  models  for  writing  new  programs,  and  all 

been  performed  on  PAX  computers  over  programs  are  in  English.  However,  without 

the  years,  in  this  author’s  opinion  some  of  one-on-one  assistance  it  would  be  difficult 

them  do  not  use  optimal  algorithms.  Fur-  for  a  user  outside  of  the  PAX  environment 

ther,  the  best  PAX  results  are  for  two-  to  utilize  the  machine.  Funding  for  the  cur- 

dimensional  problems,  despite  some  rent  project  is  questionable  after  this  year, 

extremely  clever  algorithmic  implementa-  At  that  time  the  Anritsu  Corporation  will 

tions  by  Hoshino  and  others  (see  Ref  5).  take  over  the  project  and  may  repackage 

Nevertheless,  their  examples  indicate  the  and  market  the  computer.  This  will  require 

scope  of  potential  applications.  A  number  a  major  commitment  to  refinement  and 

of  specific  types  of  computations  seem  very  maintenance  of  the  software.  I  have  no  infor- 

well  suited  to  the  PAX  architecture,  e.g.,  mation  if  this  will  occur;  Hoshino  is  not 

solving  two-dimensional  Laplace  or  Poisson  aware  of  Anritsu’s  intentions  either.  If  PAX 

equations  using  five-point  discretization  and  is  commercialized,  Anritsu  will  have  to  make 

periodic  boundary  conditions  because  the  use  of  the  expertise  in  the  PAX  group,  which 

numerical  solution  at  point  P  is  determined  could  diffuse  away  if  funding  is  not  con- 

by  data  at  immediately  adjacent  points.  tinued.  Hoshino  states  that  he  is  interested 

Similarly,  spin  systems,  because  they  are  in  adding  a  global  addressing  capability.  He 

infinitely  expandable  homogeneous  systems,  notes,  correctly,  that  in  many  problems  it  is 

are  appropriate  to  the  periodic  boundaries  convenient  to  have  a  single,  large  data  set 

of  PAX  computers.  Further,  these  systems  available  to  all  the  PUs.  The  major  super- 

are  controlled  by  the  proximity  property  computers  use  this  approach.  Hoshino’s 

and,  hence,  most  of  the  communication  is  research  will  continue  in  that  direction.  There 

between  nearest  neighbors.  seem  to  be  good  opportunities  for  scientists 

The  two-dimensional  nearest  neigh-  from  the  United  States  and  elsewhere  to 

bor  approach  of  PAX  has  somewhat  limited  visit  and  conduct  joint  research.  As  Hoshino 

its  applicability  at  the  same  time  that  it  has  has  been  developing  parallel  computers  to 

made  the  design  easier  to  build  and  to  pro-  solve  real  problems  for  many  years,  this 

gram.  Its  fast  floating  point  unit  is  a  great  would  be  very  productive, 

boost  to  performance,  but  it  has  caused 


ONRFE  SCI  INFO  BUL  15  (2) 


64 


REFERENCES 

1.  N.  Christ,  “QCD  machines,”  Nuclear  Physics  B 
(Proc.  Suppl.)  9, 549-556  (1989). 

2.  R.  Hockney  and  C.  Jesshope,  Parallel  Computers 
2  (Adam  Hilger,  Bristol  and  Philadelphia,  1988). 

3.  T.Hoshino,  “An  invitation  to  the  world  of  PAX,” 
Computer,  69-79  (May  1986). 

4.  T.  Hoshino,  R.  Hiromoto,  S.  Sekiguchi,  and 

5.  Majima,  “Mapping  schemes  of  the  particle-in -cell 
method  implemented  on  the  PAX  computer,”  Parallel 
Computing  9(1),  53-75  (December  1988). 

5.  T.  Hoshino,  Y.  Sato,  and  Y.  Asamoto,  “Parallel 
Poisson  solver  FAGECR--  Implementation  and 
performance  evaluation  on  PAX  computer,”/  Infor¬ 
mation  Processing  12(1),  20-26  (1988). 


6.  T.  Hoshino,  “Development  of  super  parallel 
computer  PAX  and  its  applications,”  paper  pre¬ 
sented  at  the  Internationa]  Conference  on  Super- 
computing  in  Nuclear  Applications,  SNA  *90. 

7.  Y.  Iwasaki,  K.  Kanaya,  T.  Yoshie,  T.  Hoshino, 
T.  Shirakawa,  Y.  Oyanagi,  S.  Ichii,  and  T.  Kawai, 
“Status  of  QCDPAX,”  paper  presented  at  the  Int 
Symp.  on  Lattice  Field  Theory,  LATTICE  89, 
September  1989,  Capri,  Italy. 

8.  Y.  Iwasaki,  T.  Hoshino,  T.  Shirakawa,  and 
Y.  Oyanagi,  “QCDPAX:  A  parallel  computer  for 
lattice  QCD  simulation,”  Computer  Physics 
Communications  49, 449-455  (1988). 

9.  T.  Shirakawa,  T.  Hoshino,  Y.  Oyanagi,  Y.  Iwasaki, 
T.  Yoshie,  K.  Kanaya,  S.  Ichii,  and  T.  Kawai, 
“QCDPAX  -  An  MIMD  array  of  vector  processors 
for  the  numerical  simulation  of  quantum  chromody¬ 
namics,"  in  Proc.  Supercomputing  ‘89, 13-17  November 
1989,  Reno,  Nev.,  ACM  Order  Number  415892. 


David  K.  Kahaner  joined  the  staff  of  the  Office  of  Naval  Research  Far  East  as  a  specialist  in  scientific 
computing  in  November  1989.  He  obtained  his  Ph.D.  in  applied  mathematics  from  Stevens  Institute  of 
Technology  in  1 968.  From  1 978  until  1 989  Dr.  Kahaner  was  a  group  leader  in  the  Center  for  Computing  and 
Applied  Mathematics  at  the  National  Institute  of  Standards  and  T echnology,  formerly  the  National  Bureau  of 
Standards.  He  was  responsible  for  scientific  software  development  on  both  large  and  small  computers. 
From  1 968  until  1 979  he  was  in  the  Computing  Division  at  Los  Alamos  National  Laboratory.  Dr.  Kahaner  is 
the  author  of  two  books  and  more  than  50  research  papers.  He  also  edits  a  column  on  scientific  applications 
of  computers  for  the  Society  of  Industrial  and  Applied  Mathematics.  His  major  research  interests  are  in  the 
development  of  algorithms  and  associated  software.  His  programs  for  solution  of  differential  equations, 
evaluation  of  integrals,  random  numbers,  and  others  are  used  worldwide  in  many  scientific  computing 
laboratories. 


90 


ONRFE  SCI  INFO  BUL  15  (2) 


65 


THIRD  NOBEYAMA  WORKSHOP 
ON  COMPUTATIONAL  FLUID 
DYNAMICS,  SUPERCOMPUTERS, 
AND  TURBULENCE 

Hideo  Yoshihara 

Background  and  description  of  the  Nobeyama  workshop  are  given  together  with  a 
summary  of  some  of  the  presentations.  The  workshopwas  composedof30invited 
papers  on  new  developments  in  turbulence  and  computational  fluid  dynamics 
(CFD). 


INTRODUCTION 

The  Third  Biennial  Nobeyama  Work¬ 
shop  (as  the  first  two)  was  held  in  the 
Yatsugatake  Kogen  Lodge,  a  resort  hotel 
outside  the  small  town  of  Nobeyama  in  the 
Japan  Alps,  a  2-hour  bus  ride  northwest  of 
Tokyo.  The  workshop  was  composed  of  30 
invited  papers.  The  theme  of  the  meeting 
was  new  concepts  of  turbulence  and  the 
interplay  of  computational  fluid  dynamics 
(CFD)  and  experiments.  Co-chairmen  of 
the  workshop  were  Professor  Kunio 
Kuwahara  of  the  Institute  for  Space  and 
Astronautical  Sciences  (ISAS)  and  Profes¬ 
sor  Steve  Orszag  of  Princeton  University. 
They  are  shown  in  Figure  1  in  front  of  the 
meeting  place,  which  normally  hosts  sum¬ 
mer  concerts.  (The  sign  translates  to 
Nobeyama  Workshop  Meeting  Place.)  There 
were  30  invitees  from  Europe  and  the  United 
States  and  49  from  Japan,  including  many 
young  researchers  from  research  institutes, 
universities,  and  companies.  The  program 
would  suggest  the  purpose  of  the  workshop 


to  be  exposing  Japanese  attendees  to  the 
expertise  of  the  invited  foreign  speakers. 
The  workshop  also  served  to  establish  per¬ 
sonal  relationships  among  the  more  experi¬ 
enced  attendees  that  would  promote  future 
exchange  of  ideas. 

Sponsorship  of  the  workshop  is  unique 
in  today’s  environment.  It  was  formally 
hosted  by  the  Keisan  Ryutai  Ryokugaku 
Kyokai  (KRRK)  (Computational  Fluid 
Dynamics  Society),  a  relatively  informal 
organization  composed  of  about  30 
researchers  from  Japanese  universities  and 
research  institutes.  The  president  is  Profes¬ 
sor  Kunio  Kuwahara  of  ISAS  and  the  vice 
president  is  Professor  K.  Nakahashi  of  Osaka 
Prefectural  University.  Membership  in  this 
organization  is  restricted.  Annual  dues  are 
¥50,000  (currently  about  $350).  The  theme 
of  the  workshop  and  the  invitees  were  pri¬ 
marily  chosen  by  Professor  Kuwahara.  The 
principal  functions  of  the  KRRK,  aside  from 
organizing  the  biennial  workshop,  are  to 
provide  travel  support  for  young  researchers 
to  attend  foreign  meetings  and  to  reprint 
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papers  from  foreign  conferences.  The  major 
source  of  funds  to  support  these  activities  as 
well  as  the  workshop  is  from  Fujitsu,  Hitachi, 
and  NEC,  the  three  supercomputer  manu¬ 
facturers  in  Japan. 


Figure  1.  Co-chairmen,  Professor  K.  Kuwahara 
(ISAS)  and  Professor  S.  Orszag 
(Princeton)  in  front  of  the  workshop 
auditorium. 

THE  LECTURES 

S.  Orszag  (Princeton  U.),  “Renormal¬ 
ization  Group  Theory  of  Turbulence.” 
Professor  Orszag  and  Dr.  Yakhot  extended 
an  earlier  renormalization  group  (RNG) 
turbulence  theory  by  evolving  turbulence 
models  at  several  levels  of  the  modeling 


hierarchy.  Line-by-line  it  is  a  difficult  theory 
to  follow.  Professor  Orszag  provided  a  more 
global,  comprehensive  explanation  of  the 
method  and  gave  impressive  examples. 
Comments  after  the  talk  criticized  gaps  in 
the  theory,  which  were  acknowledged  by  the 
speaker,  but  logic  gaps  are  not  uncommon 
in  even  the  most  respected  turbulence 
theories.  Two  examples  were  for  a  tran¬ 
sonic  airfoil  flow  with  a  shock  yielding  excel¬ 
lent  test/theory  agreement  and  a  large  eddy 
simulation  (LES)  channel  flow  solution 
matching  closely  with  direct  simulation  results 
with  l/60th  of  the  computing  time. 

H.  Sato  (Inst,  of  Flow  Research, 
Tokyo),  “Does  Turbulence  Keep  the  Memory 
of  Birth.”  Transition  was  induced  in  a  two- 
dimensional  wake  by  a  low  frequency  coher¬ 
ent  sound  wave  from  a  loud  speaker.  This 
low  frequency  coherent  component  persisted 
in  the  turbulence  spectrum  far  downstream 
(100  initial  wake  heights).  It  was  then  sug¬ 
gested  that  after  the  initial  trigger  distur¬ 
bance  was  damped,  there  was  rebirth  of  the 
coherent  component  in  the  far  downstream 
wake  from  seeds  remaining  from  the  initial 
sound  disturbance.  Others,  however,  would 
suggest  that  the  initial  disturbance  had  not 
yet  died  out  and  that  turbulence  sufficiently 
far  downstream  forgets  details  of  its  birth. 

Fazle  Hussain  (U.  of  Houston),  “Direct 
Simulation  of  Reconnection  Mechanism.” 
Reconnection  is  the  reconfiguration  of  two 
aligned  line  vortices  of  opposite  rotation 
into  two  U-shaped  hairpin  vortices.  This 
observed  phenomenon  was  reproduced  by 
Professor  Hussain  by  direct  numerical  simu¬ 
lation.  Professor  Nishioka  (Osaka  Prefec- 
tural  U.),  who  was  unable  to  attend  the 
workshop,  had  earlier  postulated  a  transi¬ 
tion  mechanism  in  terms  of  a  hierarchy  of 
hairpin  vortices.  He  expressed  satisfaction 
in  a  subsequent  meeting  that  Professor 
Hussain  was  able  to  reproduce  the  hairpin 
vortices  numerically. 
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K.  Kuwahara  (ISAS),  “Simulation  of 
Unsteady  Flows  and  Their  Visualization.” 
There  is  little  doubt  that  Professor 
Kuwahara ’s  graphics  capabilities  rank  among 
the  very  best  so  far  as  clarity  of  details  and 
scope  of  data  portrayed.  Data  displayed 
were  generated  by  direct  Navier/Stokes 
simulation.  In  the  turbulent  case,  open 
questions  remain  regarding  the  role  of 
numerical  viscosity  in  the  results. 

J.  Sethian  (U.  of  Calif.-Berkeley), 
“Modeling  of  Viscous,  Incompressible  Flow 
in  Arbitrary  Geometries  on  a  Massively 
Parallel  Supercomputer.”  Planar  incom¬ 
pressible  viscous  flow  equations  were  pro¬ 
gramed  on  the  64K  processor  Connection 
Machine  CM-2  using  a  random  vortex  method 
that  is  well  suited  for  the  massively  parallel 
computer.  Problems  of  large  vortex  struc¬ 
tures,  mixing,  and  shedding  were  calculated 
and  graphically  portrayed  on  a  video  tape. 
Relative  merits  of  the  CM-2  and,  for  exam¬ 
ple,  the  Cray  YMP8/32  depend  on  the  flow 
model  and  algorithm.  Further  experience 
with  massively  parallel  computers  is  needed 
to  gauge  their  relative  performance. 

N.  Zabusky  (Rutgers  U.),  “‘DAVID’ 
Environment  and  Application  to  Shock- 
Stratified  Flows  and  Compressible  Turbu¬ 
lence.”  DAVID  (Data  Visualization  and 
Diagnostics)  is  a  graphics  software  program 
that  “quantifies,  diagnoses,  and  validates” 
large  numerical  data  sets.  Its  use  is  illus¬ 
trated  using  the  planar  solution  of  the  inter¬ 
action  of  shock-generated  vorticity  on  inclined 
density-stratified  interfaces.  The  versatility 
of  DAVID,  permitting  composite  data  dis¬ 
plays  including  experimental  data,  was 
impressive.  Users  must  await  the  three- 
dimensional  version  presently  in  develop¬ 
ment. 

Two  papers  dealing  with  unstructured 
grid  generation  for  complex  configurations 
for  use  in  finite  volume  calculations  are 


T.  Baker  (Princeton  U.),  “Unstructured 
Mesh  Generation  and  Its  Application  to 
Complete  Aircraft  Flowfield  Computation,” 
and  K.  Nakahashi  (Osaka  Prefectural  U.), 
“An  Automatic  Flow  Solver.”  (In  addition 
to  the  grid  generation  part  a  Navier/Stokes 
flow  solver  was  also  described.)  These  two 
presentations  are  less  fundamental  fluid 
dynamically  than  the  other  papers  previ¬ 
ously  described,  but  they  are  invaluable 
contributions  for  aerospace  applications. 

Two  papers  summarizing  extensive 
aerospace  calculations  are  V.  Shankar 
(Rockwell  Int.),  “Supercomputing  in 
Aerospace-Opportunities  for  Multidisciplin¬ 
ary  Integration  of  Aerodynamics”,  and 
B.  Wagner  (Dornier  Co.),  “Some  Recent 
Industrial  Developments  and  Applications 
of  CFD.” 

In  all  workshops  of  this  series  there 
have  been  progress  reports  on  the  latest 
supercomputers  from  the  three  sponsor 
companies.  In  this  meeting  presentations 
were  given  by  Mr.  T.  Watanabe  (NEC)  and 
Mr.  K.  Uchida  (Fujitsu),  who  are  senior 
architects  in  the  respective  managements. 
Mr.  Kawabe  (Hitachi)  was  unable  to  be 
present.  Mr.  Uchida  and  Mr.  Watanabe 
gave  details  of  the  new  VP-2600  and  SX-3 
supercomputers,  respectively.  The  VP-2600 
(1  CPU,  4  G FLOPS  peak  CPU  speed)  will 
be  available  in  the  summer  of  1990,  while 
the  SX-3  (4  CPU,  22  GFLOPS  peak  speed) 
will  be  available  in  September  1990. 

Other  papers  in  the  order  of  presenta¬ 
tion  are: 

•  V.  Yakhot  (Princeton  U.),  “Probability 
Distributions  in  Turbulence” 

•  J.  Hill  (Iowa  State  U.),  “DI A  Theory  and 
Full  Turbulence  Simulations  of  Heat 
Transfer  in  Decaying,  Homogeneous 
Turbulence” 
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•  R.  Temam  (U.  of  Paris),  “The  Nonlinear 
Galerkin  Method  in  CFD” 

•  A.  Patera  [Mass.  Inst,  of  Tech.  (MIT)], 
“Parallel  Spectral  Element  Simulation  of 
Incompressible  Flow” 

•  H.  Keller  (Cal  Tech),  “Wavy  Taylor  Vortex 
Flows  Via  Multigrid-Continuation 
Methods” 

•  L.  Tuckerman  (U.  of  Texas-Austin), 
“Numerical  Simulation  of  Rayleigh- 
Benard  Convection  in  a  Cylindrical 
Geometry” 

•  E.  Krause  (RWTH-Aachen),  “Compari¬ 
son  of  Time-Dependent  Flow  Computa¬ 
tions  with  Experiments” 

•  O.  Inoue  (Tohoku  U.),  “Effects  of 
Multiple-Frequency  Forcing  on  Turbu¬ 
lent  Mixing  Layers” 

•  M.  Israeli  (Technion-Israel),  “Accurate 
Splitting  Schemes  for  Time  Dependent 
Flows” 

•  F.  Baetke  (Convex  Corp.),  “Optimizing 
the  Performance  of  Flow  Codes  on  Vec¬ 
tor  Computers” 

•  H.  Yoshihara  (Office  of  Naval  Research 
Far  East),  “Navier-Stokes  Algorithms  and 
Supercomputers  for  Aerospace  Applica¬ 
tions” 

•  T.  Aki  (National  Aeronautical  Lab- 
Mitaka),  “A  Numerical  Experiment  on 
Implosions  of  Polygonally  Interacting 
Shocks  and  Consecutive  Explosion  in  a 
Square  Box” 


•  R.  Pletcher  (Iowa  State  U.),  “Observa¬ 
tions  on  the  Numerical  Solution  of  the 
Compressible  Navier-Stokes  Equations 
at  Low  Mach  Numbers” 

•  J.  Ferziger  (Stanford  U.),  “Composite 
Grid  Methods”  (a  surprise  presentation 
from  a  turbulence  expert) 

•  C.  Chu  (Columbia  U.),  “Some  Comments 
on  the  Computation  of  Global  Ocean 
Circulations” 

•  R.  Brown  (MIT),  “Large  Scale  Compu¬ 
tation  of  Viscoelastic  Flows” 

•  K  Ghia  (U.  of  Cincinnati),  “Low-Reynolds 
Number  Nearly  Chaotic  Flows  Past  Air¬ 
foils  at  High  Angle  of  Attack” 

•  U.  Ghia  (U.  of  Cincinnati),  “Active  Control 
of  Two-Dimensional  Separated  Flow  by 
Unsteady  Forcing”  (presented  by  K  Ghia) 

Three  cancelled  presentations  that  would 

have  added  greatly  to  the  workshop  were: 

•  J.  Kim  (NASA-Ames),  “Understanding 
the  Basic  Physics  of  Turbulent  and  Tran¬ 
sitional  Flows  Using  Direct  Numerical 
Simulations” 

•  C.  Liu  and  C.  Hsu  (NASA-Langley), 
“Numerical  Studies  of  Vortex  Flow  Aero¬ 
dynamics  Using  Incompressible  Navier/ 
Stokes  Equations” 

•  S.  Osher  (U.  of  Calif.-Los  Angeles), 
“Numerical  Front  and  Shock  Capturing” 
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SUMMARY  COMMENTS 


The  subject  matter  and  speakers  at 
the  third  workshop  were  excellent.  Many  of 
the  foreign  participants  attended  all  three 
workshops  and  commented  on  the  continued 
improvements,  both  in  terms  of  the  physical 
arrangements  and  the  presentations. 
Undoubtedly  an  important  intangible  result 
of  the  workshop  has  been  the  close  friend¬ 
ships  developed,  as  witnessed  during  the 
workshop,  that  have  led  to  continuing 
international  exchange  of  ideas. 
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INTERNATIONAL  MEETINGS  IN  THE  FAR  EAST 

1990-1995 

Compiled  by  Yuko  Ushino 

The  Japan  Convention  Bureau,  the  Science  Council  of  Japan,  and  journals  of  professional  societies  are  the 
primary  sources  for  this  list.  Readers  are  asked  to  notify  us  of  any  upcoming  international  meetings  and  exhibitions  in 
the  Far  East  which  have  not  yet  been  included  in  this  report 


1990 


Date 

Tit le / At tendance * 

Site 

Contact  for  Information 

May 

16-18 

Symposium  on  Neural-Networks; 
Alliance  and  Perspective  in 
Senri  (SYNAPSE) 

Osaka, 

Japan 

Senri  International  Information  Inst. 
15F,  Bankyu-Grand  Building 

8-47  Kakuta-cho 

Klta-ku,  Osaka  530 

May 

16-19 

The  4th  Multinational 
Instrumentation  Conference 
(MICONEX  ' 90)  on  Improving 
Productivity  and  Quality  with 
Instrumentation  and  Control 
Systems 

Beijing, 
People ' a 
Republic 
of  China 

Instrument  Society  of  America 

P.0.  Box  12277 

67  Alexander  Drive 

Research  Triangle  Park,  NC  27709 

May 

17-19 

1990  Yukawa  International 
Seminar  "Common  Trends  in 
Mathematic*  and  Quantum  Field 
Theories"  (Seminar) 

Kyoto, 

Japan 

Secretariat 

1990  Research  Institute  for  Fundamental 
Physics  (RIFP) 

Kyoto  University 

Olwake-cho,  Kitashlrakawa 

Sakyo-ku,  Kyoto  606 

May 

20-25 

The  9th  International 

Symposium  on  Carotenoids 

Kyoto, 

Japan 

Profeasor  Masayoshi  Ito 

Kobe  Women's  College  of  Pharmacy 

4-19-1  Motoyamaklta -Machl 

Higashlnada-ku,  Kobe  658 

May 

20-25 

The  17th  International 
Symposium  on  Space  Technology 
and  Science 

Tokyo, 

Japan 

Ms.  Hlroko  Sakurai 

17th  ISTS  Secretariat 
c/o  Institute  of  Space  and 

Astronautical  Science 

3-1-1  Yoshlnodal 

Sagamihara,  Kanagawa  229 

May 

21-22 

Conference  and  Exhibition: 
Foundry  Asia  '90 

Hong  Kong 

FMJ  International  Publications 

May 

21-23 

4th  Symposium  on  Our 
Environment 

Singapore 

Wong  Ming  Keong 

Dept,  of  Chemistry 

National  University  of  Singapore 
Singapore  0511 


♦Note:  Data  format  was  taken  from  the  Japan  International  Congress  Calendar 
published  by  the  Japan  Convention  Bureau. 

No.  of  participating  countries 
F:  No.  of  overseas  participants 
J:  No.  of  Japanese  participants 
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Date 


Title /Attendance 


Site 


Contact  foe  Information 


June 

11-15 

1990  International  Conference; 
Metallurgical  Coatings 

Beijing, 
People's 
Republic 
of  China 

Chinese  Society  of  Metals 

46  Dongs ixi 

Dajie,  Beijing  100711 

June 

11-15 

1990  International  Conference: 
Special  Melting 

Beijing, 
People's 
Republic 
of  China 

Chinese  Society  of  Metals 

46  Dongslxl 

Dajie,  Beijing  100711 

June 

15-20 

The  2nd  International 
Conference;  Aluminum  Alloys  - 
Physical  and  Mechanical 
Properties 

Beijing, 
People's 
Republic 
of  China 

Beijing  University  of  Aeronautics 
and  Astronautics 

Beijing  100083 

June 

18-22 

C-MRS  International  '90 

Beijing, 
People's 
Republic 
of  China 

Professor  Hengde  Li 

C-MRS  International  '90 

Tsinghua  University 

Beijing  100084 

June 

19-21 

The  1990  Coal  Handling  and 
Utilization  Conference 

Sydney, 

Australia 

The  Conference  Manager 

Coal  Handling  fc  Utilisation  Conf.  1990 
The  Institution  of  Engineers,  Australia 
11  National  Circuit 

Barton,  ACT  2600 

June 

22-26 

International  Conference  on 
Dynamics,  Vibration,  and 
Control 

Beijing, 
People ' s 
Republic 
of  China 

Professor  Wei  Jlnduo 

Chinese  Society  of  Theoretical  and 
Applied  Mechanics 

No.  15  Zhong  Cuancun  Street 

Beijing 

June 

24-27 

The  5th  Japan-U.S.  Conference 
on  Composite  Materials 

Tokyo, 

Japan 

Professor  Akira  Kobayashi 

Department  of  Materials  Science 

Faculty  of  Engineering 

University  of  Tokyo 

7-3-1  Hongo 

Bunkyo-ku,  Tokyo  113 

June 

24-27 

Pacific  Offshore  Mechanics  - 
1990 

Seoul, 

Korea 

Secretariat:  PAC0MS  '90,  Korea 
Committee  for  Ocean  Resources  and 
Engineering  (KCORE) 
c/o  Dong-A  University 

840  Hadan-dong 

Saha-ku,  Pusan  604-714,  Korea 

June 

24-29 

The  3rd  International 

Symposium  and  Workshop, 

Society  of  Chinese 
Bloscientlsts  of  America 

Hong  Kong 

Professor  C.Y.  Lee 

c/o  Department  of  Biochemistry 

Chinese  University  of  Hong  Kong 

Shat in,  Hong  Kong 

June 

26-30 

International  Symposium  on 

High  Temperature  Corrosion 
and  Protection 

Shenyang, 
People's 
Republic 
of  China 

Professor  Man  Yongfa 

Institute  of  Metal  Research 

Academia  Sinlca 

2-6  Wenhua  Road 

Shenyang,  Liaoning  Province 

June 

27-29 

International  Symposium  on 
Chemistry  of  Microporous 
Crystals  (CMPC) 

Tokyo, 

Japan 

Professor  Tomoyuki  Inui 

Secretary,  CMPC 

Department  of  Hydrocarbon  Chemistry 
Faculty  of  Engineering 

Kyoto  University 

Kyoto  606 
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1990 

Data 

Title /Attendance 

Site 

Contact  for  Information 

June 

28- 

July 

3 

IAGC/IMA  Joint  Symposium  on 
Geochemistry,  Mineralogy,  and 
Origin  of  Ore  Deposits  In 
Relation  to  Future  Materials 

Beijing, 
People's 
Republic 
of  China 

IAGC 

c/o  Alberta  Research  Council 

P.0.  Box  8330 

Postal  Station  F 

Edmonton  T6H  5X2,  Canada 

July 

1-5 

The  1st  Tokyo  Conference  on 
Advanced  Catalytic  Science  and 
Technology  (TOCAT  1) 

20-F100-J200 

Tokyo, 

Japan 

Secretariat:  TOCAT  1 

c/o  Department  of  Synthetic  Chemistry 

Faculty  of  Engineering 

Tokyo  University 

7-3-1  Hongo 

Bunkyo-ku,  Tokyo  113 

July 

1-6 

The  3rd  International 
Conference  on  Technology  of 
Plasticity  (3rd  ICTP) 

10-F300- J700 

Kyoto, 

Japan 

The  Organizing  Committee  3rd  ICTP 
c/o  The  Japan  Society  for  Technology  of 
Plasticity 

Torlkatsu  Building 

5-2-5  Roppongl 

Minato-ku,  Tokyo  106 

July 

2-6 

The  26th  Yamada  Conference 
"Surface  as  a  New  Material" 

Toyonaka , 
Osaka, 

Japan 

Professor  A.  Okiji 

Secretariat  of  YCS'90 

Department  of  Applied  Physics 

Osaka  University 

2-1  Yamadaoka,  Suita  565 

July 

5-6 

IEEE  International  Workshop  on 
Intelligent  Robots  and 

Systems  '90  (IROS ' 90) -Towards 
a  New  Frontier  of  Applications 

Tsuchiura, 

Japan 

Mr.  Masakatsu  Fujie 

IROS '90 

Mechanical  Engineering  Research  Lab 
Hitachi,  Ltd. 

502  Kandatsu-machi 

Tsuchiura-shi,  Ibaraki-ken  300 

July 

6-7 

The  1st  KSME-JSME  Fracture 
and  Strength  Conference 
(Fracture  and  Strength  '90) 

Seoul, 

Korea 

Professor  Hideakl  Takahashl 

Research  Institute  for  Strength  and 
Fracture  of  Materials 

Tohoku  University 

Aoba  Tsurumakl  Aza 

Sendai  980 

July 

9-11 

Japan-U.S.A.  Symposium  on 
Flexible  Automation  -  A 

Pacific  Rim  Conference 

Kyoto, 

Japan 

Professor  Toshlhiro  Tsumura 
c/o  Institute  of  Systems,  Control 
at  Engineers 

16  Yoshida-Kawahara-cho 

Sakyo-ku,  Kyoto  606 

July 

9-12 

The  3rd  International 
Conference  on  the  Structure 
of  Surfaces 

Shanghai, 
People ' s 
Republic 
of  China 

Dr.  M.A.  Van  Hove 

MCSD,  Lawrence  Berkeley  Laboratory 
Berkeley,  CA  96720 

July 

10-13 

The  10th  UOEH  International 
Sympoaium  and  the  1st  Pan- 
Pacific  Conference  on 
Occupational  Ergonomics 

12-F50-J150 

Kitakyushu, 

Japan 

Secretary  General:  Dr.  M.  Kumashiro 
Department  of  Ergonomics 

ILES,  UOEH 

1-1  Iseigaoka 

Yahatanishi-ku,  Kitakyushu  807 

July 

11-13 

The  5th  International 

Conference  on  Manufacturing 
Engineering 

Wollongong, 

Australia 

The  Conference  Manager 

The  Institution  of  Engineers,  Australia 
11  National  Circuit 

Barton,  ACT  2600 
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1990 

Date 

Title/ Attendance 

Site 

Contact  for  Information 

July 

11-13 

The  3rd  Optoelectronics 
Conference  (OEC  '90) 

8-F20-J350 

Tokyo, 

Japan 

Katsuyoshi  Ito 

OEC  '90  Publicity  (  Registration 
Subcommittee  Chair 

c/o  Business  Center  for  Academic  Soc. 
Japan 

Conference  Department,  Crocevia 

Crocevla  Hongo  2F 

3-23-1  Hongo 

Bunkyo-ku,  Tokyo  113 

July 

12-14 

Meeting  on  Advanced  Research  on 
Computers  in  Education 

Tokyo , 

Japan 

M.  Hosaka 

Kikaishinko-kai 

3-5-8  Shiba  Koen 

Mlnato-ku,  Tokyo  105 

July 

12-14 

The  2nd  International 

Symposium  on  Magnetic  Bearings 

Tokyo , 

Japan 

Professor  Toshiro  Higuchl 

Institute  of  Industrial  Science 

Tokyo  University 

7-22-1  Roppongl 

Mlnato-ku,  Tokyo  106 

July 

12-17 

IMA,  the  15th  Quadrennial 
General  Meeting 

Beijing, 
People 1 s 
Republic 
of  China 

Professor  Huang  Yunhui 

Institute  of  Mineral  Deposits 

Chinese  Academy  of  Ceologlcal  Sciences 
Balvanzhuang  Road  26 

Fuchengmarwai ,  Beijing 

July 

15-21 

The  10th  International 

Congress  of  Nephrology 

10-F1 ,000-J4, 000 

Tokyo, 

Japan 

Japanese  Society  of  Nephrology 
c/o  2nd  Department  of  Internal  Medicine 
School  of  Medicine,  Nippon  University 
30-1  Oyaguchi-kamicho 

Itabashl-ku,  Tokyo  173 

July 

16-19 

1990  International  Micro- 
Process  Conference 
(MicroProcess  ’90) 

Makuhari, 

Japan 

Secretariat 

c/o  Business  Center  for  Academic 
Societies  Japan 

Conference  Department 

3-23-1  Hongo 

Bunkyo-ku,  Tokyo  113 

July 

16-20 

Pacific  Congress  on  Marine 
Science  and  Technology 
(PACON  90) 

10-F100-J130 

Tokyo , 

Japan 

PACON  90 

College  of  Science  and  Technology 

Nihon  University 

1-8-14  Surugadai,  Kanda 

Chlyoda-ku,  Tokyo  101 

July 

16-21 

ISEC  '90  International 

Solvent  Extraction  Conference 

31 -F150-J150 

Kyoto, 

Japan 

Conference  Secretariat  ISEC  '90 
Department  of  Chemistry 

Science  University  of  Tokyo 

Kagurazaka,  Shlnjuku-ku,  Tokyo  162 

July 

18-20 

Advanced  Research  on 

Computers  in  Education 

Tokyo, 

Japan 

Professor  Setsuko  Otsukl 

Faculty  of  Computer  Science  and  Systems 
Engineering 

Kyushu  Institute  of  Technology 

1-1  Sensui-cho,  Tobata-ku 

Kltakyushu-shi,  Fukuoka  804 
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July 

18-21 

The  3rd  International 

Symposium  on  Hunan  Factors 

In  Organizational  Design 
and  Management 

15-F100- J100 

Kyoto, 

Japan 

Dr.  Junzo  Watada 

Symposium  Secretariat 

3rd  Internat'l  Symp.  on  Human  Factors 
in  Organizational  Design  4  Management 
c/o  Faculty  of  Business  Administration 
Ryukoku  University 

Fukakusa,  Fushlmi,  Kyoto  612 

July 

19-21 

The  22nd  International 

Congress  of  Applied  Psychology 
(22nd  ICAP)  (Workshop) 

Kyoto, 

Japan 

Professor  Yasuhisa  Nagayama 

General  Secretariat  22ICAP 

P.0.  Box  38 

Suita-Senri,  Osaka  565 

July 

20-24 

International  Conference  on  Fukuoka, 

Fuzzy  Logic  and  Neural  Networks  Japan 
(IIZUKA'90) 

Professor  Takeshi  Yamakawa 

Dept,  of  Computer  Science  and  Control 
Engineering 

Kyushu  Institute  of  Technology 

680-4  Kawazu,  Ohaza 

Iizuka-shi,  Fukuoka  820 

July 

22-25 

The  5th  Asia-Pacific 
Confederation  of  Chemical 
Engineering  (APCChE  1990) 

Kuala  Lumpur, 
Malaysia 

Conference  Sec'y  APCChE  Conference  1990 
c/o  The  Institution  of  Engineers, 
Malaysia 

P.0.  Box  223  (Jin  Sultan) 

46720  Petaling  Jaya  Malaysia 

July 

22-27 

The  22nd  International 

Congress  of  Applied  Psychology 
(22  ICAP) 

Kyoto, 

Japan 

Professor  Yasuhisa  Nagayama 

General  Secretariat,  22ICAP 

P.0.  Box  38 

Suita-Senri,  Osaka  565 

July 

24-26 

The  29th  SICE  Annual  Conference  Tokyo, 

(SICE1 90)  Japan 

Dr.  Suguru  Arlmoto 

c/o  Society  of  Inatrument  and  Control 
Engineers  (SICE) 

1-35-28-303  Hongo 

Bunkyo-ku,  Tokyo  113 

July 

29- 

August 

3 

IAWPRD  BIENNIAL  '90 
(International  Association  on 
Water  Pollution  Research  and 
Control) 

Kyoto, 

Japan 

Jun-ichlro  Mataumoto,  Chairman, 

Japanese  Organizing  Committee 
c/o  Jpn  Soc  on  Water  Pollution  Research 
Yotsuya  New  Mansion  307 

12  Honshio-cho 

Sbinjuku-ku,  Tokyo  160 

July 

30- 

August 

2 

The  15th  International 
Conference  on  International 
Association  on  Water  Pollution 
Research  and  Control 

Kyoto, 

Japan 

Japan  Soc.  on  Water  Pollution  Research 
and  Control 

Yotsuya  New  Mansion 

12  Honahiocho 

Shinjuku-ku,  Tokyo  173 

Auguat 

1-4 

The  1st  International 
Conference  on  Mechanics  and 
Physics  of  Advanced  Materials 
and  Structures 

Singapore 

Professor  G.C.  Sih 

Institute  of  Fracture  i  Solid  Mechanics 
Packard  Laboratory  Building  No.  19 
Lehigh  University 

Bethlehem,  PA  18015 

Auguat 

2-8 

The  25th  International 
Conference  on  High  Energy 
Physics  1990 

Singapore 

Professor  K.K.  Phua 

South  East  Asia  Theoretical  Physics 
Association 
c/o  Dept,  of  Physics 

National  University  of  Singapore 

Kent  Ridge,  Singapore  0511 
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August 

3-6 

International  Conf.  on  the 
Environmental  Management  of 
Enclosed  Coastal  Seas  *90 
(EMECS  ’90) 

Kobe, 

Japan 

Secretariat 

Exec.  Comm,  of  the  Internet '1  Conf. 
on  the  Environmental  Management  of 
Enclosed  Coastal  Seas 

5-10-1  Shlmoyamate-dorl 

Chuo-ku,  Kobe  650 

August 

6-9 

International  Conference  on 
Advanced  Materials  Mechanical 
Properties  '90  (ICAMP’90) 

Utsunomiya, 

Japan 

Secretariat 

ICAMP' 90 

P.0.  Box  1234,  Shlbuya 

Shibuya-ku,  Tokyo  150 

August 

6-10 

The  Sth  International  Meeting 
of  the  International  Humic 
Substances  Society  (IHSS5) 

Nagoya, 

Japan 

Registration  Secretariat:  5th  IHSS 
c/o  Japan  Convention  Services,  Inc. 
Nagoya  Branch 

Nagoya  International  Bldg,  19th  Floor 
1-47-1  Nagono 

Nakamura -ku,  Nagoya  450 

August 

7-9 

International  Symposium  on 
Cereal  and  Other  Plant 
Carbohydrates 

15-F40-J250 

Kagoshima, 

Japan 

Organizing  Committee  of  International 
Symposium  on  Cereal  and  Other  Plant 
Carbohydrates 

c/o  Dept  of  Agricultural  Chemistry 
Faculty  of  Agriculture 

Kagoshima  University 

1-21-24  Korlmoto 

Kagoshima  890 

August 

7-10 

IUPAB  Satellite  Congress 

Palmerston 

North, 

New  Zealand 

J.  Tigyl 

Institute  of  Biophysics 

Medical  University 

Szigeti  ut  12 

7643  Pecs,  Hungary 

August 

7-11 

International  Symposium  on 
Analytical  Chemistry 

Changchun, 
People ' 8 
Republic 
of  China 

Professor  Qlnhan  Jin 

Dept,  of  Chemistry 

Changchun,  China 

August 

7-12 

6th  International  Symposium 
on  Biology  of  Turbellaria 

24-F95-J70 

Hirosaki, 

Japan 

Professor  Wataru  Teshlrogl 

Department  of  Biology 

Faculty  of  Science 

Hirosaki  University 

3  Bunkyo-cho 

Hirosaki,  Aomori  036 

August 

12-17 

The  15th  International 
Carbohydrate  Symposium 

25-F200-J800 

Yokohama, 

Japan 

Dr.  Ishido,  General  Secretary 

Faculty  of  Science 

Tokyo  Institute  of  Technology 

Ookayama,  Meguro-ku,  Tokyo  152 

August 

12-18 

The  14th  International 

Congress  of  Soil  Science 

86-F1 , 200-J600 

Kyoto, 

Japan 

Secretariat  of  14th  ICSS 

6-26-10-202  Hongo 

Bunkyo-ku,  Tokyo  113 

August 

13-17 

The  4th  Asia  Pacific  Physics 
Conference 

Seoul, 

Korea 

Program  Committee,  AAPC 

Department  of  Physics 

Yonsel  University 

Seoul  120-749,  Republic  of  Korea 
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August 

13-18 

Tsukuba  International 
Conference  on  Representations 
of  Algebras  and  Related  Topics 

13-F50-J80 

Tsukuba, 

Japan 

Institute  of  Mathematics 

University  of  Tsukuba 

1-1-1  Tennodai 

Tsukuba,  Ibarakl  305 

August 

15-19 

International  Conference  "Knot 
Theory  and  Related  Topics" 

20-F100-J150 

Osaka, 

Japan 

Professor  A.  Kawauchi 

Faculty  of  Science 

Osaka  City  University 

3-3-138  Suglmoto 

Sumlyoshi-ku,  Osaka  558 

August 

16-18 

SIGNAL  International  Symposium 
on  Algorithms 

Tokyo, 

Japan 

Professor  Tetsuo  Asano 

Department  of  Applied  Electronics 

Osaka  Electro-Communications  University 
18-8  Hatsu-cho 

Neyagawa-shi,  Osaka  572 

August 

18-20 

General  Assembly,  Internet' 1 
Mathematical  Union 

52-F124-J6 

Kobe  9 

Japan 

ICM  90  Secretariat 

c/o  Research  Institute  for  Mathematical 
Sciences 

Kyoto  University 

Oiwake-cho,  Kltashirakawa 

Sakyo-ku,  Kyoto  606 

August 

20-21 

The  7th  VLSI  Process /Device 
Modeling  Workshop 

Kawasaki, 

Japan 

Mr.  Hiroshi  Iwal 

ULSI  Research  Center,  Toshiba  Corp. 

1  Komukai-Toshlba-cho 

Salwai-ku,  Kawasaki  210 

August 

20-21 

International  Workshop  on 
Frontiers  in  Plant  and 
Microbial  Glycans:  Structure, 
Functions,  and  Biotechnology 

Kyoto, 

Japan 

Dr.  Naoto  Shibuya 

National  Food  Research  Institute 
Ministry  of  Agriculture 

Forestry  and  Fisheries 

2-1-2  Kannondal 

Tsukuba,  Ibarakl  305 

August 

20-24 

1990  International  Symp.  on 
Symbolic  and  Algebraic 
Computation 

15-F60-J140 

Tokyo, 

Japan 

ISSAC  '90  Conference  Office 
c/o  Scientist,  Inc. 

Yamazaki  Building 

3-2  Kanda-Surugadai 

Chlyoda-ku,  Tokyo  101 

August 

21-25 

Western  Pacific  Geophysics 
Meeting  (WEPGEM) 

15-F200-J600 

Kanazawa, 

Japan 

Organizing  Committee  of  WEPGEM 
c/o  Professor  I.  Kimura 

Dept  of  Electrical  Engineering 

Faculty  of  Engineering 

Kyoto  University 

Yoshlda-honmachl 

Sakyo-ku,  Kyoto  606 

August 

21-29 

International  Congress  of 
Mathematicians  1990 

84 -FI ,500- J2 ,000 

Kyoto, 

Japan 

ICM  90  Secretariat 

c/o  International  Relations  Office 

Research  Inst,  for  Math.  Sciences 

Kyoto  University 

Kltashirakawa  Oiwake-cho 

Sakyo-ku,  Kyoto  606 

August 

22-23 

International  Workshop 
"Conformational  Studies  of 
Polysaccharides" 

Osaka, 

Japan 

Dr.  Kozo  Ogawa 

Osaka  Prefectural  Radiation 

Research  Institute 

174-16,  Shinke-cho 

Sakai,  Osaka  593 
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August 

22-24 

1990  International  Conference 
on  Solid  State  Devices  and 
Materials 

N.A.-F100-J900 

Sendai, 

Japan 

c/o  Business  Center  for  Academic 
Societies  Japan 

Crocevia  Building  2F 

3-23-1  Hongo 

Bunkyo-ku,  Tokyo  113 

August 

23-30 

V  International  Congress 
of  Ecology 

70-F1.300-J700 

Yokohama, 

Japan 

Sec'y  General's  Office  for  INTECOL  1990 
c/o  Inst,  of  Environmental  Science  and 
Technology 

Yokohama  National  University 

156  Tokivadai 

Hodogaya-ku,  Yokohama  240 

August 

26-31 

AUSTCERAM  90 

Perth, 

Australia 

Conference  Secretariat,  AUSTCERAM  90 

P.0.  Box  515 

South  Perth,  WA  6151,  Australia 

August 

27-31 

The  2nd  Japan-Soviet  Union 
Joint  Symposium  on 
Computational  Fluid  Dynamics 

4-F36-J100 

Tsukuba , 

Japan 

Dr.  Y.  Yoshizawa 

Institute  of  Engineering  Mechanics 
University  of  Tsukuba 

1-1-1  Tennodai 

Tsukuba  305 

August 

27-31 

Step  into  the  Nineties  Joint 
Conference  on  Corrosion, 
Metals,  Materials,  Metal 
Finishing 

Gold  Coast, 
Australia 

The  Covener,  Step  into  the  90' a 

P.0.  Box  198 

Upper  Mount  Gravatt,  QL,  Australia  4122 

August 

28-31 

International  Conference  6 
Exhibition  on  Computer 
Applications  to  Materials 
Science  and  Engineering 
(CAMSE  90) 

Tokyo, 

Japan 

Professor  M.  Doyama 

CAMSE  '90 

c/o  The  Nikkan  Kogyo  Shlmbun,  Ltd. 
Business  Bureau 

1-8-10  Kudan  Kite 

Chiyoda-ku,  Tokyo  102 

August 

29- 

S«pt ember 

4 

The  11th  International 
Symposium  on  Biotelemetry 

N.A. -F120- J250 

Yokohama , 
Japan 

Professor  A.  Uchlyama 

Dept,  of  Electronics  &  Communication 
School  of  Science  and  Engineering 

Waseda  University 

3-4-1  Okubo 

Shinjuku-ku,  Tokyo  169 

August 

30- 

September 

4 

International  Conference  on 
Potential  Theory 

29-F100-J100 

Nagoya, 

Japan 

Secretariat 

International  Conf.  on  Potential  Theory 
c/o  Department  of  Mathematics 

College  of  General  Education 

Nagoya  University 

Furo-cho,  Chikusa-ku,  Nagoya  464-01 

August 

30- 

September 

4 

International  Symposium  on 
Computational  Mathematics 

10-F30- J50 

Matsuyama, 

Japan 

Professor  T.  Yamamoto 

Department  of  Mathematics 

Ehlme  University 

2-5  Bunkyo-machi,  Matsuyama  790 

August 

31- 

September 

4 

International  Symposium  on 
Functional  Analysis  and 

Related  Topics 

13-F40-J60 

Sapporo, 

Japan 

Department  of  Mathematics 

Faculty  of  Science 

Hokkaido  University 

Nlshi  8-chome,  Rita  10-jo 

Kita-ku,  Sapporo  060 
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September 

3-5 

International  Symposium  on 
Diagnostics  and  Modeling  of 
Combustion  in  Internal 
Combustion  Engines  (C0M0D1A  90) 

Kyoto, 

Japan 

Professor  Makoto  Ike garni 

Dept,  of  Mechanical  Engineering 

Kyoto  University 

Sakyo-ku,  Kyoto  606 

September 

4-7 

The  2nd  International 

Symposium  on  Chemical 

Synthesis  of  Antibiotics  and 
Related  Microbial  Products 

0i8O, 

Japan 

Faculty  of  Pharmaceutical  Sciences 

University  of  Tokyo 

7-3-1  Hongo 

Bunkyo-ku,  Tokyo  113 

23-F200-J500 

September 

10-14 

The  17th  Congress  of  the 
Collegium  International 

Neuro • Psychopharaaco log icura 

Kyoto, 

Japan 

The  17th  CINP  Congress 
c/o  Simul  International  Inc. 

Kowa  Building  No.  1 

1-8-10  Aka8aka 

Minato-ku,  Tokyo 

September 

12-14 

The  12th  International 

Symposium  on  Capillary 
Chromatography 

30-P100-J250 

Kobe, 

Japan 

Assistant  Professor  Klyokatsu  Jlnno 

School  of  Material  Engineering 

Toyohashl  University  of  Technology 

1-1  Hibarigaoka,  Aza 

Tenpaku-cho,  Toyohashi-shi  440 

September 

16-20 

The  4th  International 

Conference  on  Copepoda 

Karuizawa, 

Japan 

Professor  Shin-ichi  Uye 

Hiroshima  University,  Shltaml 

Shaijo-cho,  Hlgashi-Hiroshlma  724 

September 

16-22 

The  15th  IUMS  Congress: 
Bacteriology  6  Mycology  - 
Osaka,  Japan  -  1990 

7 1 -F2 , 000- J3 ,500 

OsaKa, 

Japan 

Secretary  General 

c/o  Department  of  Microbiology 

Faculty  of  Medicine 

Kyoto  University 

Yoshida,  Konoe-cho 

Sakyo-ku,  Kyoto  606 

September 

18-21 

The  3rd  Asia-Pacific  Microwave 
Conference  (APMC  ’90) 

30-F150-J350 

Tokyo , 

Japan 

APMC  '90  Secretariat 

c/o  Business  Ctr  for  Academic  Societies 

Japan 

3-23-1  Hongo 

Bunkyo-ku,  Tokyo  113 

September 

19-22 

The  2nd  World  Congress  on 
Particle  Technology 

H.A.-F100-J400 

Kyoto, 

Japan 

Secretariat:  2nd  World  Congress  on 

Particle  Technology 

c/o  Society  of  Powder  Technology,  Japan 
Shibunkaku-kaikan 

2-7  Tanakaseklden-cho 

Sakyo-ku,  Kyoto  606 

September 

20-23 

The  3rd  Shallow  Tethys 
International  Symposium 

15-F40-J60 

Sendai, 

Japan 

Secretariat  of  Shallow  Tethys  3 
c/o  Inst,  of  Geology  and  Paleontology 

Faculty  of  Science 

Tohoku  University 

Aobayama,  Sendai  980 

September 

23-28 

The  57th  World  Foundry 

Congress  (WFC) 

31-F400-J800 

Osaka, 

Japan 

Secretariat 

Japan  Foundrymen's  Society 

Toyokawa  Building 

8-12-13  Ginza 

Chuo-ku,  Tokyo  104 
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September 

24-27 

The  6th  International  Congress 
on  Polymers  in  Concrete 

Shanghai, 
People's 
Republic 
of  China 

ICPIC-90  Secretariat 
c/o  Associate  Professor  Tan  Muhua 
Institute  of  Materials  Science  and 
Engineering 

Tongji  University 

Shanghai 

September 

24-27 

The  3rd  International  Aerosol 
Conference 

29-F200-J300 

Kyoto, 

Japan 

Professor  Kanji  Takahashi 
c/o  Institute  of  Atomic  Energy 

Kyoto  University 

UJ1,  Kyoto  611 

September 

24-28 

The  12th  International 
Conference:  Boundary  Element 
Method  Conference  (BEM  12) 

15-F60-J160 

Sapporo, 

Japan 

Mr.  Hiroshi  Mizoguchi 

J AS CHOKE,  KKE  Inc. 

Dai-ichi  Selmei  Building  24F 

2-7-1  Nishi-Shinjuku 

Shinjuku-ku,  Tokyo  160 

September 

24-29 

The  5th  International 
Conference  on  Polymer 

Supported  Reactions  '.a 

Organic  Chemistry  (POC  90) 

Kyoto, 

Japan 

Dr.  Takeo  Shlmidzu 

Division  of  Molecular  Engineering 
Faculty  of  Engineering 

Kyoto  University 

Yoshida-honmachi 

Sakyo-ku,  Kyoto  606 

September 

24-29 

IAGC,  Geological  Society  of 
Australia,  the  7th 
International  Conference  on 
Geochronology,  Cosmochronology 
and  Isotope  Geochemistry 

Canberra, 

Australia 

IAGC 

c/o  Alberta  Research  Council 

P.0.  Box  8330 

Postal  Station  F 

Edmonton  T6H  5X2,  Canada 

September 

28- 

October 

2 

4th  International  Symposium 
on  Benthic  Foraminlfera: 

BENTHOS  '90 

Sendai, 

Japan 

Organizing  Comm,  for  4th  Int'l.  Symp. 

on  Benthic  Foraminlfera 
c/o  Inst,  of  Geology  and  Paleontology 
Faculty  of  Science 

Tohoku  University 

Aobayoma,  Sendai  980 

October 

1-5 

International  Conference  on 
Information  Technology 
Commemorating  the  30th 
Anniversary  of  the  Information 
Processing  Society  of  Japan 
(IPSJ)  -  InfoJapan  '90 

Tokyo , 

Japan 

InfoJapan  '90  Secretariat:  IPSJ 

Hoshina  Building  3F 

2-4-2  Azabudal 

Mlnato-ku,  Tokyo  106 

20-F200-J 1,000 

October 

1-5 

The  3rd  International  N^e 
Materials  Conference  (New 
Materials  90  Japan) 

12-F100-J300 

Osaka, 

Japan 

Secretariat:  Nev  Materials  90  Japan 
c/o  Inter  Group  Corp. 

Shohaku  Building 

6-23  Chayamachi 

Kita-ku,  Osaka  530 

October  4th  International  Conference  Tokyo,  Secretariat:  ICSP 

2*5  on  Shotpeenlng  Japan  Japan  Society  of  Precision  Engineering 

Ceramics  Building 
2-22*17  Ryakunln-cho 
Shinjuku-ku,  Tokyo  169 
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October 

8-9 

The  let  International 

Symposium  on  functionally 
Gradient  Materials 

Sendai, 

Japan 

The  Society  of  Non-Traditlonal 
Technology 

Toranooon  Kotohlra  Kaikan  Building  3F 
1-2-8  Toranomon 

Mlnato-ku,  Tokyo  105 

October 

8-10 

The  1st  International  Workshop 
on  Algorithmic  learning  Theory 
(ALT' 90) 

Tokyo, 

Japan 

Professor  Setsud  Arlkava 

Program  Chairman  of  ALT' 90 

Research  Inst,  of  Fundamental 
Information  Science 

Kyushu  University  33 

Fukuoka  812 

October 

8-11 

"Benibana"  International 
Symposium  on  How  to  Improve 
the  Toughness  of  Polymers 
and  Composites 

10-F40-J1 10 

Yamagata, 

Japan 

Professor  Ikuo  Narlsawa 

Polymer  Materials  Engineering 

Yamagata  University 

4-3-16  Jonan 

Yonezawa-shl,  Yamagata  992 

October 

8-12 

The  5th  Australasian  Remote 
Sensing  Conference 

Perth, 

Australia 

Gold  West  Conventions 

P.0.  Box  411 

West  Perth,  WA  6005 

October 

8-12 

The  6th  International 

Symposium  on  Biomineralization 

17-F55-J45 

Od aware, 
Japan 

Professor  S.  Suga 

Department  of  Pathology 

School  of  Dentistry 

The  Nippon  Dental  University 

1-9-20  Fujiml 

Chiyoda-ku,  Tokyo  102 

October 

9-12 

Fracture  and  Fatigue  of  High- 
Performance  and  Multi-Phase 
Polymeric  Materials 

8-F25-J60 

Undecided, 

Japan 

Faculty  of  Engineering 

Yamagata  University 

4-3-16  Jonan 

Yonezawa,  Yamagata  992 

October 

11-12 

The  2nd  KS ME -MS ME  Fluids 
Engineering  Conference 

Seoul, 

Korea 

Professor  Jinho  Lee 

Dept  of  Mechanical  Engineering 

College  of  Engineering 

Yonsel  University 

134  Shlnchon-dong,  Seodaemun-ku 

Seoul  120-749 

October 

12-14 

The  3rd  Japan-China  Joint 
Meeting  on  Optical  Fiber 

Science  and  Electromagnetic 
Theory 

Fukuoka , 
Japan 

Professor  Takashi  Hldaka 

Dept  of  Electronic  Engineering 

College  of  Science  &  Technology 

1-8-14  Kanda  Surugadal 

Chlyoda-ku,  Tokyo  101 

October 

14-19 

International  Conference  for 
Hew  Smelting  Reduction  and 

Near  Net  Shape  Casting 
Technologies  for  Steel 

Pohang, 

Korea 

Conference  Department 

Institute  of  Metals 

1  Carlton  House  Terrace 

London,  SW1Y  5  5DB,  U.K. 

October 

15-18 

The  1 st  Asian-Pacific 
International  Symposium  on 
Combustion  and  Energy 
Utilization 

Beij ing, 
People ' s 
Republic 
of  China 

Prof.  Huang,  Zhao  Xiang  and  Song  Jlalln 
Institute  of  Engineering  Therroophysics 
Chinese  Academy  of  Sciences 

P.0.  Box  2706,  Beijing 
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October 

15-19 

The  4th  International 

Symposium  on  Marine 

Engineering  (ISME  KOBE  '90) 

25-F120- J330 

Kobe, 

Japan 

ISME  Organizing  Committee 
c/o  Kobe  Shosen  Daigaku 

5-1-1  Fukae-Mlnami 

Hlgaahinada-ku,  Kobe  658 

October 

15-19 

The  8th  International 

Symposium  on  Vlt.  B6  and 
Carbonyl  Catalysis 

20-F80-J80 

Osaka, 

Japan 

Institute  of  Scientific  and  Industrial 
Research 

8-1  Mihogaoka 

Ibaraki,  Osaka  567 

October 

15-19 

The  11th  Meeting  of 
International  Collaboration  on 
Advanced  Neutron  Sources 
(ICANS  -  XI) 

Taukuba, 

Japan 

National  Laboratory  for  High  Energy 
Physics 

1-1  Oho 

Taukuba,  Ibaraki  305 

10-F40-J60 

October 

15-27 

International  Electrotechnical 
Commission,  the  54th  Annual 
General  Meeting 

BelJ  ing. 
People's 
Republic 
of  China 

Chinese  National  Committee  of  the  IEC 
P.0.  Box  820 

Beijing 

October 

21-26 

The  6th  International 

Iron  and  Steel  Congress 

50-F300-J500 

Nagoya, 

Japan 

International  Conference  Department 

Iron  and  Steel  Institute  of  Japan 

3F,  Keidanren  Kalkan 

1-9-4  Otemachi 

Chiyoda-ku,  Tokyo  100 

October 

22-25 

The  11th  International  Coal 
Preparation  Congress 

N.A. -F250-J150 

Tokyo, 

Japan 

Secretariat 

11th  Int'l  Coal  Preparation  Congress 
c/o  Slmul  International,  Inc. 

Kowa  Building,  No.  9 

1-8-10  Aka s aka 

Mlnato-ku,  Tokyo  107 

October 

22-25 

10th  International  Acoustic 
Emission  Symposium  1990 
(IAES-10) 

Sendai, 

Japan 

Professor  Hiroaki  Niltsuma 

Faculty  of  Engineering 

Tohoku  University 

Aramaki  Aza  Aoba,  Sendai  980 

October 

22-26 

International  Conference  on 
Signal  Processing 

Beij ing, 
People's 
Republic 
of  China 

Professor  Zong  Sha  Nongzhan  Guan 

Nan  Lu  12,  Room  2307 

Beijing  100026 

October 

22-26 

International  Conference  on 
Information  Technology  in 
Connection  with  30th 
Anniversary  Celebration  of 
Information  Processing  Society 
of  Japan 

Osaka, 

Japan 

Secretariat:  International  Conference 
on  Information  Technology 
c/o  Slmul  International,  Inc. 

Kowa  Building,  No.  9 

1-8-10  Akasaka 

Mlnato-ku,  Tokyo  107 

N.A.-F200-J1 ,000 

October 

25-31 

The  1st  Japanese  Knowledge  for 
Knowledge -Baaed  Systems 

Workshop  (JKAW) 

Kyoto, 

Japan 

Assoc.  Professor  Rtlchlro  Hizoguchi 
Institute  of  Scientific  and  Industrial 

8-1  Mihogaoka 

Ibaraki,  Osaka  567 
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Dace 
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Site 

Contact  for  Information 

October 

28- 

November 

2 

The  2nd  International 
Conference:  HSLA  Steels 

Beijing, 
People's 
Republic 
of  China 

Chinese  Society  of  Metals 

46  Dongslxi 

Dajie,  Beijing  100711 

October 

29-31 

The  5th  International  Workshop 
on  Atmospheric  Icing  of 
Structures 

20-F100-J150 

Tokyo , 

Japan 

Secretariat  for  5th  International 
Workshop  on  Atmospheric  Icing  of 
Structures 

c/o  Inter  Group  Corp. 

Aka 8 aka  Tamakatsu  Building 

8-5-32  Akaaaka 

Minato-ku,  Tokyo  107 

October 

29- 

Nov ember 

1 

Japan  International  Tribology 
Conference  Nagoya  -  '90 

N.A. -F200-J600 

Osaka, 

Japan 

Secretariat:  Japan  ITC  Nagoya  -  '90 
c/o  Toyota  Technological  Institute 
2-chome,  Hisakata 

Tempaku-ku,  Nagoya  468 

October 

29- 

November 

2 

International  Symposium  on 
Mineral  Exploration  -  Use  of 
Artificial  Intelligence,  1990 
(ISME  -  AI  '90) 

N.A. -F30-J70 

Tokyo, 

Japan 

Mining  and  Materials  Processing 

Institute  of  Japan  (MMIJ) 

Noglzaka  Building 

9-6-41  Akasaka 

Minato-ku,  Tokyo  107 

October 

29- 

Nov ember 

4 

The  14th  International  Pigment 
Cell  Conference 

12-F200-J200 

Kobe, 

Japan 

Organizing  Committee 

XIVth  International  Pigment  Cell 
Conference 

c/o  Department  of  Dermatology 

Kobe  University 

School  of  Medicine 

7-5-2  Kusunoki-cho 

Cbuo-ku,  Kobe  650 

November 

4-8 

International  Symposium  on 
Carbon,  1990:  "New  Processing 
and  New  Applications" 

15-F50-J200 

Tsukuba, 

Japan 

The  Carbon  Society  of  Japan 

Salto  Building  2F 

2-16-13  Yujima 

Bunkyo-ku,  Tokyo  113 

November 

5-8 

The  9th  International 

Symposium  on  Zirconium  in 
the  Nuclear  Industry 

20-F100-J100 

Kobe, 

Japan 

American  Society  for  Testing  and 
Materials  (ASTM) 

Conference  Department 

1916  Race  Street 

Philadelphia,  PA  19103 

November 

7-9 

IUPAC  International  Symposium 
on  Specialty  Polymers 

Singapore 

Professor  S.H.  Goh 

Department  of  Chemistry 

National  University  of  Singapore 

10  Kent  Ridge  Crescent 

Singapore  0511 

November 

11-13 

Metal  Bulletins  6th 
International  Aluminum 
Conference 

Singapore 

Metal  Bulletin  Conferences  Ltd. 

House  Park  Terrace,  Worcester  Park 
Surrey  KT4  7HT,  U.K. 
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Data 
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Contact  for  Information 

November 

26-29 

The  3rd  SPSJ  International 
Polymer  Conference  (IPC  90)  - 
Membranes  and  Interfacial 
Phenomena  of  Polymers 

Nagoya, 

Japan 

1PC  Secretariat 

c/o  Society  of  Polymer  Science,  Japan 
5-12-8  Ginza 

Chuo-ku,  Tokyo  104 

5-F100-J200 

November 

26-30 

The  5th  International 
Photovoltaic  Science  and 
Engineering  Conference 
(International  PVSEC-5) 

Kyoto, 

Japan 

Professor  Junjl  Sarale 

Secretariat  of  International  PVSEC-5 
c/o  Japan  Convention  Services,  Inc. 
Nippon  Press  Center  Building 

2-2-1  Uchlsaiwai-cho 

Chiyoda-ku,  Tokyo  100 

November 

28-30 

IAPR  Workshop  on  Machine 

Vision  Applications 

Tokyo , 

Japan 

Professor  Mikio  Takagl 

Institute  of  Industrial  Science 
University  of  Tokyo 

7-22-1  Roppongl 

Mlnato-ku,  Tokyo  106 

November 
(tent. ) 

The  3rd  International 

Symposium  on  Superconductivity 

17-F100-J700 

Undecided, 

Japan 

International  Superconductivity 
Technology  Center 

6F  Eishln  Kaihatsu  Building 

5-34-3  Shinbashl 

Mlnato-ku,  Tokyo  105 

December 

4-7 

ICCV3  Third  International 
Conference  on  Computer 

Vision 

Osaka, 

Japan 

Professor  Saburo  Tsujl 

Dept  of  Control  Engineering 

Faculty  of  Engineering  Science 

1-1  Machlkaneyama-cho 

Toyonaka,  Osaka  560 

December 

6-8 

The  6th  International 
Conference  on  Biomedical 
Engineering 

Singapore 

Dept  of  Orthopedic  Surgery 

National  University  Hospital 

Lower  Kent  Ridge  Road 

Singapore  0511 

1990 
(tent. ) 

Chemeca  1990  Applied 

The  rmodynamlcs 

New  Zealand 

Conference  Manager 

The  Institution  of  Engineers,  Australia 
11  National  Circuit 

Barton,  ACT  2600 

1991 

Date 

Title /Attendance 

Site 

Contact  for  Information 

February 

7-12 

The  10th  International 
Conference  on  Offshore 
Mechanics  and  Arctic 
Engineering 

Seoul, 

Korea 

Korea  Cmt  for  Ocean  Resources  6  Engrag 
Dong-A  University 

840  Sahagu 

Puaan,  Korea 

February 

10-15 

Polymer  91  IUPAC  International 
Symposium  -  Polymer  Materials 
Preparation,  Characterization, 
and  Properties 

Melbourne, 

Australia 

Polymer  91  Secretariat 

R.A.C.I.  Polymer  Division 

P.0.  Box  224 

Belmont,  VIC  3216,  Australia 

March 
11-15 
(tent. ) 

Computing  In  High  Energy 
Physics  '91 

19-F100-J100 

Tsukuba, 

Japan 

National  Laboratory  for  High  Energy 
Physics 

1-1  Oho 

Tsukuba,  Ibarakl  305 
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Contact  for  Information 

November 

13-16 

Annual  Meeting  of 

International  Lead  Zinc 
Research  Organization 

10-F40-J20 

Tokyo, 

Japan 

Process  Technology  Department 

Mitsubishi  Metal  Corp. 

P.0.  Box  227 

Sonic  City  Building 

1-441  Sakuragi-cho,  Omiya  331 

November 

14-15 

TECHNO-OCEAN  '90  International 
Symposium 

7-F50- J350 

Kobe, 

Japan 

Kobe  International  Association 

6-9-1  Mlnatojima-Nakamachi 

Chuo-ku,  Kobe  650 

November 

14-15 

Kyushu  International  Techno- 
Business  Forum 

7-F150-J450 

Kitakyushu, 

Japan 

Secretariat 

Kyushu  International  Techno-Business 
Forum 

c/o  Enterprise  Development  Division, 
Economic  Affairs  Bureau 

Kitakyushu  City  Office 

1-1  Jonai 

Kokura  Kita-ku,  Kitakyushu  803 

November 

14-16 

Pacific  Rim  International 
Conference  on  Artificial 
Intelligence  '90 

Nagoya, 

Japan 

Secretariat:  PRICAI  '90 

Inter  Group,  Inc. 

Aka 8 aka  Yamakatsu  Building 

8-5-32  Akasaka 

Minato-ku,  Tokyo  107 

November 

14-16 

Rare  Metals  '90 

15-F100-J200 

Kitakyushu, 

Japan 

Mining  and  Materials  Processing  Inst. 

of  Japan  (MMIJ) 

Noglzaka  Building 

9-6-41  Akasaka 

Mlnato-ku,  Tokyo  107 

November 

18-20 

I-SAIRAS  '90  -  International 
Symposium  on  Artificial 
Intelligence,  Robotics,  and 
Automation  for  Space 

Kobe, 

Japan 

Ren  Planning  A  Research  Office 

2nd  Klkke  Building 

2-12-14  Hamamatsu-cho 

Mlnato-ku,  Tokyo  105 

7-F50-J200 

November 

18-22 

1990  International  Conference 
on  Spoken  Language  Processing 
(ICSLP-90) 

Kobe, 

Japan 

Secretariat:  ICSLP-90 

c/o  Simul  International,  Inc. 

Kowa  Building  No.  9 

1-8-10  Akasaka 

Minato-ku,  Tokyo  107 

November 

19-22 

Asia  Pacific  Interfinish  90 
Congress 

Hong  Kong 

Asia  Pacific  Interfinish  90 
c/o  Australasian  Inst,  of  Metal 

Finishing 

191  Royal  Parade 

Parkville,  VIC  3052,  Australia 

November 

19-23 

Corrosion-Air,  Sea,  Soil: 

CASS  90 

Auckland, 

New  Zealand 

Professor  G.  Trabanelll 

Corrosion  Study  Center 

Department  of  Chemistry 

University  of  Ferrara 

Via  L.  Borsari  46 

1-4410  Ferrara,  Italy 
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March 

12-15 

The  4th  International 
Conference  on  Ceramic 

Powder  Processing  Science 

20-F250- J350 

Nagoya, 

Japan 

Secretariat 

4th  Int'l  Conf.  on  Ceramic  Powder 

Processing  Science 
c/o  Dept  of  Applied  Science 

Faculty  of  Engineering 

Nagoya  University 

Furo-cho,  Chikusa-ku,  Nagoya  464 

March 

18-19 

V1SIC0M  '91  (Packet  Video  91) 

Kyoto, 

Japan 

Mr.  Hideo  Hashimoto 

Visual  Media  Laboratories 

Human  Interface  Laboratories 

1-2356  Take 

Yokosuka -ahi,  Kanagawa  238-03 

March 

18-20 

IEEE  International  Conference 
on  Microelectronic  Teat 
Structures 

Kyoto, 

Japan 

Dr.  Yukinori  Kuroki 

Micro  Electronics  Research  Labs 

NEC  Corp. 

1120  Shimokuzawa 

Sagarolhara-shi,  Kanagawa  229 

March 

21-23 

i 

VISICOM  '91  (PCS  91) 

Tokyo , 

Japan 

Mr.  Hideo  Hashimoto 

Visual  Media  Laboratories 

Human  Interface  Laboratories 

1-2356  Take 

Yokosuka-shl,  Kanagawa  238-03 

April 

4-12 

The  7th  International 
Conference  on  Data 

Engineering 

Kobe, 

Japan 

Professor  Tadao  Ichikawa 

Dept  of  Circuit  6  Electrical 

System  Engineering 

Faculty  of  Engineering 

Shltami,  Saijo-cho 

Higashi  Hlro8hima-shl,  Hiroshima  724 

April 

8-12 

1FIP  TC5/WG  5.10  Working 
Conference  on  Modeling  in 
Computer  Graphics 

Tokyo , 

Japan 

Toshlyasu  L.  Kunil,  Program  Chairperson 

IFIP  TC5/WG  5.10  Working  Conference 

Dept  of  Information  Science 

University  of  Tokyo 

7-3-1  Hongo 

Bunkyo-ku,  Tokyo  113 

April 

11-13 

IMAC  '91  -  International 
Conference  on  Image  Management 
and  Communication  in  Patient 
Care 

Kyoto, 

Japan 

Dr.  Ellchi  Takenaka 

Department  of  Radiology 

National  Defense  Medical  College 

3-2  Namikl,  Tokorozawa  359 

May 

7-13 

Beijing  Essen  Welding  '91 

Beijing, 
People's 
Republic 
of  China 

Messe  Essen  Nobert  Street 

D-4300  Essen 

Federal  Republic  of  Germany 

May 

13-17 

The  5rh  International 

Symposium  on  Toxicity 
Assessment 

25-F100- J30 

Kurashiki, 

Japan 

Toxicity  Assessment:  Kurashiki 

Symposium  Committee 
c/o  Dr.  I.  Aoyama 

Research  Institute  for  Bioresources 

Okayama  University 

2-20-1  Chuo,  Kurashiki  710 

June 

3-6 

International  Conference  on 
Superplasticity  in  Advanced 
Materials  (ICSAM-91) 

Osaka, 

Japan 

Assistant  Professor  N.  Furushiro 

Dept  of  Materials  Science  and 

Engineering 

Osaka  University 

2-1  Yamadaoka 

Suita-shi,  Osaka  565 
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June 

10-12 

The  3rd  International 
Conference  on  Hot  Iaoatatlc 
Pressing 

Osaka, 

Japan 

Professor  Yoshinarl  Miyamoto 

Processing  Research  Center  for  High 
Performance  Materials 

Institute  of  Scientific  and  Industrial 
Research 

Osaka  University 

8-1  Mihogaoka 

Ibaraki-shi,  Osaka  567 

June 

10-13 

International  Conference  on 
Stainless  Steels 

20-F50-J100 

Tokyo, 

Japan 

Secretariat:  STAINLESS  STEELS  '91 

The  Iron  and  Steel  Institute  of  Japan 
Keidanren  Kaikan 

1-9-4  Otemachi 

Chlyoda-ku,  Tokyo  100 

June 

10-14 

The  4th  International 
Conference  on  Nucleus-Nucleus 
Collisions 

20-F200- J200 

Kanazawa, 

Japan 

The  Inst,  of  Physical  and  Chemical 
Research  (RIKEN) 

2-1  Hirosawa 

Wako,  Saltama  351-01 

June 

12-15 

The  18th  International 
Materials  Handling,  Storage 
and  Distribution  Exhibition 
and  Conference 

Tokyo, 

Japan 

Show  Management,  Yoahihisa  Shiraiahi 
Convention  Department,  Japan  Management 
Association 

3-1-122  Shiba  Koen 

Mlnato-ku,  Tokyo  105 

June 

17-21 

10th  International  Symposium 
on  the  Mathematical  Theory  of 
Networks  and  Systems  (MTNS-91) 

30-F180- J170 

Kobe, 

Japan 

Professor  Hidenori  Kimura 

Dept,  of  Mechanical  Engineering  for 
Computer-Controlled  Machinery 

Faculty  of  Engineering 

Osaka  University 

2-1  Yamadaoka 

Sulta-shl,  Osaka  565 

June 

23-29 

The  6th  Marcel  Grossmann 
Meeting  on  General  Relativity 

31-F380-J170 

Kyoto, 

Japan 

Secretariat 

6th  Marcel  Grossmann  Meeting  on 

General  Relativity 
c/o  Professor  F.  Sato 

Department  of  Physics 

Faculty  of  Science 

Kyoto  University 

Kitashirakawa  Olwake-cho 

Sakyo-ku,  Kyoto  606 

June 

30- 

July 

4 

The  6th  Southeast  Asian/ 
Western  Pacific  Regional 
Meeting  of  Pharmacologists 

Hong  Kong 

Professor  C.W.  Ogle 

Hong  Kong  Pharmacology  Society 

Department  of  Pharmacology 

Faculty  of  Medicine 

University  of  Hong  Kong 

Pokfularo  Road,  Hong  Kong 

June 
(tent. ) 

JIMIS-6:  Intermetallic 

Compound  -  Properties  and 
Applications 

Tokyo, 

Japan 

Professor  Osamu  Waizumi 

Institute  for  Materials  Research 

2-1-1  Katahlra 

Sendai  980 

July 

1-5 

International  Conference  on 

Hot  Carriers  in  Semiconductors 

Nara, 

Japan 

Department  of  Electronics 

Osaka  University 

2-1  Yamadaoka 

10-F70-J80  Suita,  Osaka  565 
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July 

7-12 

World  Congress  on  Medical 
Physics  and  Biomedical 
Engineering  (9th  ICMP  - 
16th  IFMBE) 

54 -FI , 000- J1 ,500 

Kyoto, 

Japan 

Secretariat 

World  Congress  on  Medical  Physics  and 

Biomedical  Engineering 
c/o  Japan  Convention  Services,  Inc. 

Kansai  Branch 

Sumitomo  Selmel  Midosuji  Building 

4-14-3  Kishltemma 

Kita-ku,  Osaka  530 

July 

8-13 

The  6th  International 
Conference  on  Chemistry  of 
Selenium  and  Tellurium 

12-F30-J120 

Osaka, 

Japan 

Professor  N.  Sonoda,  Chairman  of 

Organizing  Committee  of  6th  Int'l. 

Conf.  on  Chemistry  of  Selenium  and 

Tellurium 

Dept  of  Applied  Chemistry 

Faculty  of  Engineering 

Osaka  University 

2-1  Yamadaoka 

Suita,  Osaka  565 

July 

15-17 

The  7th  International 
Conference  on  Vapor  Growth 
and  Epitaxy 

15-F60- J350 

Nagoya, 

Japan 

ICVGE-7  Committee 

c/o  Professor  I.  Akasaki 

Department  of  Electronics 

Faculty  of  Engineering 

Nagoya  University 

Furo-cho,  Chikusa-ku,  Nagoya  464-01 

July 

24-26 

The  3rd  International 
Conference  on  Residual 

Stresses  (ICRS-3) 

Tokushima, 

Japan 

Society  of  Materials  Sciences,  Japan 

1-101  Yoshida  Izumidono-cho 

Sakyo-ku,  Kyoto  606 

30-F150-J200 

July 

24-30 

The  17th  International 
Conference  on  the  Physics  of 
Electronic  and  Atomic 
Collisions 

Brisbane, 

Australia 

Dr.  W.R.  Newell 

Department  of  Physics 

University  College  of  London 

Gower  Street 

London  WC1E  6BT  UK 

July 

29- 

Auguet 

The  6th  International 
Conference  on  Mechanical 
Behavior  of  Materials  (ICM-6) 

Kyoto, 

Japan 

Society  of  Materials  Sciences,  Japan 

1-101  Yoshida  Izumidono-cho 

Sakyo-ku,  Kyoto  606 

30-F300-J300 

August 

18-23 

The  11th  International 
Conference  on  Structural 
Mechanics  In  Reactor 

Technology 

45-F300-J700 

Tokyo, 

Japan 

SMirt-1 1 

c/o  Komae  Research  Laboratory 

Central  Research  Institute  of 

Electric  Power  Industry 

2-11-1  Iwato-kita 

Komae,  Tokyo  201 

August 

25-31 

IUPAC  International  Congress 
on  Analytical  Science-1991 
(ICAS  '91) 

25-F500-J1 ,000 

Chiba, 

Japan 

Japan  Society  for  Analytical  Chemistry 

Rm  304  Gotanda  Sun  Heights 

1-26-2  Nishi  Gotanda 

Shlnagawa-ku,  Tokyo  141 

August 

28-30 

1991  International  Conference 
on  Solid  State  Devices  and 
Materials 

Yokohama, 

Japan 

c/o  Business  Center  for  Academic 

Societies  Japsn 

Crocevla  Building  2F 

3-23-1  Hongo 

Bunkyo-ku,  Tokyo  113 
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September 

5-10 

The  12th  IMEKO  World  Congress 
(IMEKO  XII) 

Beijing, 
People's 
Republic 
of  China 

Instrument  Society  of  America 

P.0.  Box  12277 

67  Alexander  Drive 

Research  Triangle  Park,  NC  27709 

September 

8-13 

The  5th  International 
Conference  on  II-VI  Compounds 

Okayama, 

Japan 

Professor  Makoto  Konagai 

Electrical  and  Electronic  Engineering 
Faculty  of  Engineering 

2-12-1  Ookayana 

Meguro-ku,  Tokyo  152 

September 

12-14 

IMEKO  TC7  8th  International 
Symposium  on  Artificial 
Intelligence  Based  Measurement 
and  Control  (AlMaC  '91) 

25-F50-J100 

Kyoto, 

Japan 

Professor  Komyo  Kariya 

Dept  of  Electrical  Engineering 

Faculty  of  Science  and  Engineering 
Titsumelkan  University 

56-1  Tojin-Kita 

Klta-ku,  Kyoto  603 

September 

24-28 

The  9th  International 

Symposium  on  Affinity 
Chromatography  and  Biological 
Recognition  (9th  AC  i  BR) 

23-F150-J200 

Yokohama, 

Japan 

9th  AC  i  BR  Organizing  Committee 
c/o  Professor  K.  Kasai 

Faculty  of  Pharmaceutical  Sciences 
Teikyo  University 

Sagamlko  Tsukul,  Kanagawa  199-01 

September 

25-28 

International  Congress  on 
Flexible  Manufacturing  in 

Sheet  Metal  Working 

Hong  Kong 

Uta  Bates 

Mack-Brooks  Exhibitions,  Ltd. 

Forum  Place,  Hatfield  ALIO  ORN,  U.K. 

September 

26-29 

The  3rd  International 
Conference  on  Steel-Concrete 
Composite  Structures 

7-F30-J270 

Fukuoka, 

Japan 

Department  of  Architecture 

Faculty  of  Engineering 

Kyushu  University 

6-10-1  Hakozakl 

Higashl-ku,  Fukuoka  812 

September 

27-30 

The  2nd  Osaka  International 
Colloqulm  on  Ship  Viscous  Flow 
<OC  ’91) 

Osaka, 

Japan 

Professor  T.  Okuno 

Department  of  Naval  Architecture 

College  of  Engineering 

University  of  Osaka  Prefecture 
Nozu-Umeoachi,  Sakai,  Osaka  591 

September 

29- 

October 

4 

The  5th  Aslan  Pacific  Congress 
of  Clinical  Biochemistry 
(5th  APCCB) 

20-F300- J600 

Kobe, 

Japan 

Secretariat:  5th  APCCB 
c/o  Central  Laboratory  for  Clinical 
Investigation 

Osaka  University  Hospital 

1-1-50  Fuku8hlma 

Fukushlma-ku,  Osaka  553 

September 

30- 

October 

3 

IFAC/ISHS  Workshop  on 
Mathematical  and  Control 
Applications  in  Agriculture 
and  Horticulture 

Matsuyama, 

Japan 

Assistant  Professor  Hiroshi  Nonaml 

Dept  of  Biomechanical  Systems 

Ehlme  University 

3-5-7  Tarumi 

Matsuyama-shi,  Ehiroe  790 

October 

2-4 

International  Conference  on 

Ion  Exchange  '91 

Tokyo , 

Japan 

Secretariat 

c/o  Japan  Assoc,  of  Ion  Exchange  (JAIE) 
Department  of  Chemistry 

Faculty  of  Science 

2-12-1  Ookayama 

Meguro-ku,  Tokyo  152 
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Date 

Tit le / Attendance 

Site 

Contact  for  Information 

October 

14-16 

International  Conference  on 
Recent  Advances  in  Science 
and  Engineering  of  Light 

Metals  (RASELM  *91)  in 
Commemoration  of  40th 
Anniversary  of  the  Japan 
Institute  of  Light  Metals 

Sendai, 

Japan 

Secretariat:  RASELM  '91 
c/o  Professor  K.  Ikeda 

Dept  of  Materials  Processing 

Faculty  of  Engineering 

Tohoku  University 

Aoba,  Sendai  980 

October 

15-19 

1991  International  Symposium 
on  Lattice  Field  Theory 

20-F140- J30 

Tsukuba, 

Japan 

National  Laboratory  for  High  Energy 
Physics 

1-1  Oho 

Tsukuba,  Ibarakl  305 

October 

21-23 

Small  Engine  Technology 
Conference  (SETC) 

Yokohama , 
Japan 

Gobancho  Center  Building  5F 

Society  of  Automotive  Engineers  of 

Japan,  Inc. 

10-2  Goban-cho 

Chiyoda-ku,  Tokyo  102 

October 

27- 

Nov ember 

2 

The  6th  International  Pacific 
Conference  on  Automotive 
Engineering 

Seoul, 

Korea 

IPC-6  Conference  Secretariat 

Korea  Society  of  Automotive  Engineers, 
Inc. 

1638-3  Socho-dong,  Socho-ku 

Seoul  137 

October 

28-31 

International  Conference  on 
Fast  Reactors  and  Fuels  Cycles 

8-F150- J350 

Kyoto, 

Japan 

Power  Reactor  &  Nuclear  Fuel 

Development  Corp. 

1-9-13  Akasaka 

Mlnato-ku,  Tokyo  107 

October 

28- 

Nov ember 

1 

International  Conference  on 
Industrial  Electronics, 

Control,  and  Instrumentation 
(IECON  ’91) 

Kobe, 

Japan 

Society  of  Instrument  and  Control 
Engineers  of  Japan  (SICE) 

#303,  1-35-28  Hongo 

Bunkyo-ku,  Tokyo  113 

30-F200- J400 

November 

5-8 

The  13th  International 
Telecommunications  Energy 
Conference  (INTELEC  ’91) 

Kyoto, 

Japan 

Mr.  Katsuichi  Yotsumoto 

Applied  Electronics  Laboratories 

3-9-11  Midori-cho 

Musashlno-shl,  Tokyo  180 

November 

18-22 

International  Conference  on 
Accelerator  and  Large 
Experimental  Physics  Control 
Systems 

Tsukuba, 

Japan 

National  Laboratory  for  High  Energy 
Physics 

1-1  Oho 

Tsukuba,  Ibarakl  305 

15-F90- J80 

November 
(tent. ) 

The  1 st  Japan-China  Joint 
Seminar  of  Filtration  and 
Separation 

Shanghai  or 
Hangchou, 
People's 
Republic 
of  China 

The  Society  of  Chemical  Engineers, 

Japan 

4-6-19  Kohinata 

Bunkyo-ku,  Tokyo  112 
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1992 

Date 

Title/ Attendance 

Site 

Contact  for  Information 

February 

(tent.) 

The  19th  Australian 

Polymer  Symposium 

Perth, 

Australia 

RACI  Polymers  Division 

P.0.  Box  224 

Belmont,  VIC  3216 

May 

17-22 

NETWORKS  '92:  The  5th 
International  Network 

Planning  Symposium 

Kobe, 

Japan 

NTT  Telecommunication  Networks  Labs 
3-9-11  Midorl-cho 

Musashlno-shl,  Tokyo  180 

30-F200-J150 

August 

3-8 

The  8th  International  Congress 
of  Biorheology 

20-F300-J200 

Yokohama , 
Japan 

8th  ICBR  (Int'l  Cong.  Blorheology) 

Local  Organizing  Committee 
c/o  Department  of  Physics 

Keio  University 

832  Hlyoshi-cho 

Kohoku-ku,  Yokohama  223 

Augusc 

24- 

September 

3 

The  29th  International 
Geological  Congress 

30-Fl,000-J2,000 

Kyoto, 

Japan 

IGC  Secretariat 

P.0.  Box  65 

Tsukuba,  Ibaraki  305 

Augusc 

30- 

September 

4 

The  9th  International  Congress 
on  Photosynthesis 

20-F500- J500 

Nagoya, 

Japan 

Professor  Norio  Murata 

Okazaki  National  Research  Institute 
National  Research  for  Basic  Biology 

38  Salgou-Naka,  Mlyoudaigi-cho-aza 
Okazaki,  Aichl 

September 

22-25 

1992  International  Symposium 
on  Antennas  and  Propagation 
(ISAP-92) 

Sapporo, 

Japan 

Dr.  Ken-lchi  Kagoshima 

RfcD  Laboratories 

NTT 

1-2356  Take 

Yokosuka-shi,  Kanagawa  238-03 

October 

26-30 

The  14th  International 
Switching  Symposium  (ISS  ’92) 

60-F1 , 200- J800 

Yokohama , 
Japan 

NTT  Communication  Switching  Labs 

3-9-11  Mldori-cho 

Musashlno-shl,  Tokyo  180 

November 

9-12 

The  8th  International  Congress 
on  Heat  Treatment  of  Materials 
(HEAT  i  SURFACE  '92) 

N.A. -F-J500 

Kyoto, 

Japan 

(tentative) 

Secretariat  of  8th  Int'l  Congress  on 
Heat  Treatment  of  Materials 
c/o  Research  Inst,  for  Applied  Science 

49  Tanaka  Ohl-cho 

Sakyo-ku,  Kyoto  606 

1993 

Date 

Title/ Attendance 

Site 

Contact  for  Information 

February 

16-20 

International  Federation  of 
Clinical  Chemistry,  the  15th 
Triennial  World  Congress 

Melbourne , 
Australia 

Dr.  David  Hay 

Royal  Womens  Hospital 

132  Gratten  Street 

Carlton,  VIC  3053,  Australia 

May 

23-28 

The  18th  International  Mineral 
Processing  Congress 

Sydney, 

Australia 

AUSIMM,  Conference  Department 

P.0.  Box  122 

Parkville,  VIC  3052 
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1993 


Dace 

Tit le /Attendance 

Site 

Contact  for  Information 

July 

12-15 

1993  International  Powder 
Metallurgy  Conference  & 
Exhibition 

40-P200-J400 

Kyoto, 

Japan 

Japan  Powder  Metallurgy  Association 
Tamagawa  Building 

2-2-16  Ivamoto-cho 

Chiyoda-ku,  Tokyo  101 

August 

8-13 

The  8th  International  Congress 
of  Virology 

Osaka, 

Japan 

Prof.  Stuart  W.  Glover,  Sec'y  General 
Int'l  Union  of  Microbiological  Soc. 

Dept  of  Genetics,  Catherine  Cookson 

Bldg 

Framllngton  Place 

Nerocastle  upon  Tyne  NE2  4 HU,  UK 

Late 

August 

The  25th  Congress  of  IAHR 
(International  Assoc,  of 
Hydraulic  Research) 

40-F400-J600 

Tokyo , 

Japan 

Local  Organizing  Committee  of  25th 

IAHR  Cong. 

c/o  Professor  Y.  Iwasa 

School  of  Civil  Engineering 

Kyoto  University 

Yoshida-honmachl 

Sakyo-ku,  Kyoto  606 

1993 

(tent.) 

International  Federation  of 
Automatic  Control  Congress 

Sydney, 

Australia 

Conference  Manager 

The  Institution  of  Engineers,  Australia 
11  National  Circuit 

Barton,  ACT  2600 

1994 

Date 

Title / Attendance 

Site 

Contact  for  Information 

Tentative 

XXX  International  Conference 
on  Coordination  Chemistry 

Kyoto, 

Japan 

Professor  Hltoshl  Ohtakl 

Coordination  Chemistry  Laboratories 
Institute  for  Molecular  Science 
Myodaijl-cho,  Okazaki  444 

Tentative 

The  10th  International 
Conference  on  the  Strength 
of  Metals  and  Alloys 
(ICSMA-10) 

Undecided, 

Japan 

Professor  Hiroshi  Oikava 

Faculty  of  Engineering 

Tohoku  University 

Aoba,  Aramakl  Aza 

Sendai  980 

1995 

Date 

Title /Attendance 

Site 

Contact  for  Information 

Tentative 

The  13th  International  Vacuum 
Congress  (IVC-13) 

7th  International  Conference 
on  Solid  Surfaces  (ICSS-7) 

Undecided, 

Japan 

Vacuum  Society  of  Japan 

302  Klkal  Shinko  Kalkan  Annex 

3-5-22  Shiba-koen 

Mlnato-ku,  Tokyo  105 

Undecided 

World  Energy  Council  '95 

90-F-J5.000 

Undecided, 

Japan 

Japan  Energy  Association 

Uchisalval  Building 

1-4-2  Uchisalval -cho 

Chiyoda-ku,  Tokyo  100 

Yuko  Ushino  is  a  technical  information  specialist  for  ONR  Far  East.  She  received  a  B.S.  degree  from 
Brigham  Young  University  at  Provo,  Utah. 
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